Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,500 T→C 6.5% intergenic (‑113/+54) otsB ← / ← araH trehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,5000TC 6.5% 165.1 / 3.2 65intergenic (‑113/+54)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

GGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGA  >  NC_000913/1982352‑1982646
                                                                                                                                                    |                                                                                                                                                  
ggATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGt                                                                                                                                                    >  2:77531/1‑149 (MQ=255)
 gATAGTTAAGGGTTGTCGGTTAACGGTGCTGGCACCCGTTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGtt                                                                                                                                                   <  1:684113/149‑1 (MQ=255)
      ttCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGCACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGc                                                                                                                                              >  2:484172/1‑149 (MQ=255)
          ggggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                          >  2:395007/1‑149 (MQ=255)
          ggggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGtt                                                                                                                                          <  1:289744/149‑1 (MQ=255)
           gtgTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGAGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTAAGGGCGGTTg                                                                                                                                         >  2:972787/3‑149 (MQ=255)
           gggTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATAGTGTTGCGAAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGCCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTg                                                                                                                                         >  1:8674/1‑149 (MQ=255)
             ggttCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACGGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                       <  1:415446/147‑1 (MQ=255)
             gTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGCCATAACGGTTGGc                                                                                                                                       <  2:287243/149‑1 (MQ=255)
             gTTTCGGTTAACGGTTCTGAAACCCGGTTCTCCTACTTATCTTTTCTTATTGTGTTGCGCACAAAGCGCTGTACAGACATTTAAGTATAGACAGTGTGACGGGGGTCGCCATGTGCAACGCAGAACAGCCAAACGATCTGGCGGTGGGc                                                                                                                                       >  2:926743/1‑149 (MQ=255)
                 ctgTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATTGACAGTGTGACGGGTGTCGCCATTTGCAACGAAGATCAGCCAACCGTTATGGCGGTTGGCTAtg                                                                                                                                   <  1:747633/147‑1 (MQ=255)
                 cGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTCGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtg                                                                                                                                   <  1:289694/149‑1 (MQ=255)
                  ggTTCACGGTGCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtgg                                                                                                                                  <  1:43501/149‑1 (MQ=255)
                   gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTTTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCACAGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                 >  1:844045/1‑149 (MQ=255)
                   gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                 >  2:246771/1‑149 (MQ=255)
                   gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTATCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAtggt                                                                                                                                 >  2:91526/1‑149 (MQ=255)
                   gTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTGTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCGTTTGCAAGGAAGAACAGCCACCCGTTATGGCGGTAGGCTAgggt                                                                                                                                 >  2:580288/1‑149 (MQ=255)
                          gTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaa                                                                                                                          >  2:433411/1‑149 (MQ=255)
                           ttCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGTAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaaa                                                                                                                         >  2:683545/1‑149 (MQ=255)
                                   cccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTaaa                                                                                                                 >  1:94061/1‑149 (MQ=255)
                                         tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGCCATGTAAGTATAGACAGTGGGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAAttgttg                                                                                                           <  2:577097/149‑1 (MQ=255)
                                         tctcCTCCTTCTCTTTTCTTATTGTCTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGGCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAAttgttg                                                                                                           <  2:220483/149‑1 (MQ=255)
                                           tcctccTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGAAACGAAGAACAGCCATCCGTTATGGCGCTTGGCTATGGTGGGCAAAAACGCTAAATTGTTCCa                                                                                                         >  2:807959/1‑149 (MQ=255)
                                                 ttCTGTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCAAACCGTTATGGCGTTTGGCTATGGTGGGAAAACACGCTAAATTGTAGCAGaaaaa                                                                                                   >  1:574512/1‑149 (MQ=255)
                                                          cTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGa                                                                                          <  2:770576/149‑1 (MQ=255)
                                                             aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCTACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAg                                                                                       >  1:284483/1‑149 (MQ=255)
                                                             aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAg                                                                                       >  1:288433/1‑149 (MQ=255)
                                                                   ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt                                                                                 <  1:520891/149‑1 (MQ=255)
                                                                   ttGCGCAGAACGCGCGTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGtt                                                                                 >  1:448442/1‑149 (MQ=255)
                                                                     gcgcAGAAGGCGCTTTACAGAAATGTAAGTATCGACAGTGTGACGGGGTTCGCGATTTGCAACGAAGAACAGCCAACCGTTTTGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTc                                                                               <  1:588360/149‑1 (MQ=255)
                                                                      cgcAGCTCGCGCGATACAGCAAGGCAAGTATAGACAGTTGGAGCGGGTGCGCCTTTTGCAACGAAGCACAGCCAACCGTTATGGCGGTTTGCTATGGGGGGAAAACTCGCTAAATTGTTGCAGAAAAAAGCAGCCGACAGTGCGtttgg                                                                              <  2:673394/149‑3 (MQ=255)
                                                                      cgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCg                                                                              <  2:201512/149‑1 (MQ=255)
                                                                         acaaCGCGCTTTACAGACATGTAAGTCTAGACAGGGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAATAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtt                                                                           <  2:597816/147‑1 (MQ=255)
                                                                          gAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGTAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttt                                                                          <  1:683545/149‑1 (MQ=255)
                                                                           aaCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtttt                                                                         >  1:678023/1‑149 (MQ=255)
                                                                           aaCGCGCTTTACAGACATGAAAGTATAGACAGTCAGACCGGCCTCGCCATATGCAACGAAGAACAGACAACAGTTATGGCGGTTGGCTATGGTGTGAATAACCGCTAAATTGTTGCATAAAAAAGCATCAGACAGTACCATTCGCtttt                                                                         >  2:465706/1‑149 (MQ=255)
                                                                            aCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttttt                                                                        <  1:555546/149‑1 (MQ=255)
                                                                                gcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCtt                                                                    <  1:895083/149‑1 (MQ=255)
                                                                                  tttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGt                                                                  <  1:246771/149‑1 (MQ=255)
                                                                                       aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACgg                                                             >  1:243737/1‑149 (MQ=255)
                                                                                          cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGGTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCTGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCg                                                          <  1:77531/149‑1 (MQ=255)
                                                                                           aTGTAAGAATAGACAGTGTGACGGGGGTCCACATATGCAACGAAGAACAGCCAACCGATATGGAGGTTGGCTATGGAGGGAAAAAACGCTAAATTGATGCAGAAAAAAGCATCAGACAGAGCGTTTCCCTATTAGCCTGTAACGGACGa                                                         >  1:234481/1‑149 (MQ=255)
                                                                                            tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGaa                                                        <  2:390171/149‑1 (MQ=255)
                                                                                              tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGa                                                      >  1:169433/1‑149 (MQ=255)
                                                                                               aattataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGAt                                                     <  2:844045/146‑1 (MQ=255)
                                                                                               aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACCAAGAACAGCCAACCGTAATGGCGGGTGGCTATGGTGGGAAAAAACGCTAAATTGTTTCAGAAAAAAGCATCAGACAGTGCGTTGCGCTTTTTGCTTGTAACGGACGAAGAt                                                     <  2:97321/149‑1 (MQ=255)
                                                                                                                cGGGCGTCGCCATTTGCAACGAAGAACAGCCCAACGTTATGGCGGGTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAACGCCTCAGACAGCGCGGTTCGCTTTTTGCTTGTAACGGGCGAAGAGCACCGCTGCCAGCAGGa                                    <  2:151314/149‑1 (MQ=255)
                                                                                                                   gggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                 >  2:432597/1‑149 (MQ=255)
                                                                                                                     gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAg                               <  1:484172/149‑1 (MQ=255)
                                                                                                                        gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcg                            <  1:844545/149‑1 (MQ=255)
                                                                                                                         ccATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc                           >  1:882350/1‑149 (MQ=255)
                                                                                                                               gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAATCCTCGAACCa                     >  1:1024998/1‑149 (MQ=255)
                                                                                                                                cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                    >  2:423347/1‑149 (MQ=255)
                                                                                                                                cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAc                    <  1:91526/149‑1 (MQ=255)
                                                                                                                                cAACGAAGAACAGCCAACCGTAATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGAAGGATTAAGCCCCGaaaaac                    >  2:425507/1‑145 (MQ=255)
                                                                                                                                 aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCCCg                   >  2:472469/1‑149 (MQ=255)
                                                                                                                                   cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGCCAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTa                 >  2:483004/1‑149 (MQ=255)
                                                                                                                                   cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTa                 >  2:1024174/1‑149 (MQ=255)
                                                                                                                                     aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAATACGCGAAATTGTTGCAGAAAAAAACATCAGACAGTGCGTTTCGGTGTTTGCTTGTAACGATCGAACTTCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt               >  1:984769/1‑149 (MQ=255)
                                                                                                                                     aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACt               <  2:956960/149‑1 (MQ=255)
                                                                                                                                      accataGCCAACCGCCATAACGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTg              >  2:126825/6‑149 (MQ=255)
                                                                                                                                       ccataGCCAACCGCCATAACGGTGGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTGTTTGCTTGTAACGGGCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCCCGTACTgc             >  1:429259/5‑149 (MQ=255)
                                                                                                                                        ctgagtcAAGCGTTATTGCGGTTGGATTTGGTGGGAAAAAACGCTATATTGTTGCAGAATAAAGCATCAGTCAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCAGGTACTgcg            <  2:1024998/143‑1 (MQ=255)
                                                                                                                                           aGCCACCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGCACCCGTTACTGCGCGa         >  1:159364/1‑149 (MQ=255)
                                                                                                                                             ccAACAGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa       <  1:545175/149‑1 (MQ=255)
                                                                                                                                              cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg      <  2:94061/149‑1 (MQ=255)
                                                                                                                                                  cGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTCTAACGGTACAACAACACCTCTCCCAGCAGGATTAAGCCGAGAACCACGTACTTCGCGAAAGgaga  >  1:154692/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGATAGTTCAGGGGTTTCGGTTAACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGA  >  NC_000913/1982352‑1982646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: