Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,153,214 T→A 7.5% intergenic (+34/‑86) fabF → / → pabC 3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,2140TA 7.5% 814.1 / 42.1 240intergenic (+34/‑86)fabF/pabC3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (114/108);  new base A (10/8);  total (124/116)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 242 total aligned reads displayed.

GATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCAGTT  >  NC_000913/1153067‑1153361
                                                                                                                                                   |                                                                                                                                                   
 aTCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGc                                                                                                                                                   >  2:224927/1‑149 (MQ=255)
          tCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtttt                                                                                                                                          <  4:213218/149‑1 (MQ=255)
           cGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCttttt                                                                                                                                         <  3:205167/149‑1 (MQ=255)
           cGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCttttt                                                                                                                                         >  1:90900/1‑149 (MQ=255)
            gTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttttat                                                                                                                                        <  3:835159/149‑3 (MQ=255)
                cGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTa                                                                                                                                    >  1:728559/1‑149 (MQ=255)
                cGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTa                                                                                                                                    <  3:1576660/149‑1 (MQ=255)
                    cGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCtt                                                                                                                                <  4:61786/149‑1 (MQ=255)
                     gAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTc                                                                                                                                       >  1:11633/1‑141 (MQ=255)
                     gAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTc                                                                                                                                       <  2:11633/141‑1 (MQ=255)
                     gAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCttt                                                                                                                               <  1:1732804/149‑1 (MQ=255)
                         cgcgTCAGGGTAGCGGAATGGAATACACACTGTCTAACTCCTTCGGCTTCCGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTAtt                                                                                                                           <  2:1432427/149‑1 (MQ=255)
                             tCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGa                                                                                                                       <  4:1443817/149‑1 (MQ=255)
                             tCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGa                                                                                                                       <  4:687461/149‑1 (MQ=255)
                              cAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGAc                                                                                                                      <  1:1254270/149‑1 (MQ=255)
                                           tggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTg                                                                                                         >  3:105511/1‑149 (MQ=255)
                                            ggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGa                                                                                                        >  2:1666133/1‑149 (MQ=255)
                                             gaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttt                                                                                                                                           >  2:1391946/1‑113 (MQ=255)
                                             gaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttt                                                                                                                                           <  1:1391946/113‑1 (MQ=255)
                                                 acacTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTGCCGTAGTGAACAtg                                                                                                   <  2:675646/149‑1 (MQ=255)
                                                   acTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgct                                                                                                 >  3:443956/1‑149 (MQ=255)
                                                     tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCgga                                                                                                                                                >  3:1526380/1‑99 (MQ=255)
                                                     tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTATTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctgc                                                                                               <  2:274814/149‑1 (MQ=255)
                                                          gtAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCacac                                                                                          <  1:560165/149‑1 (MQ=255)
                                                                 cttcggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGtt                                                                                                                 <  1:1216879/119‑1 (MQ=255)
                                                                 cttcggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTa                                                                                                            >  2:74626/1‑124 (MQ=255)
                                                                 cttcggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAAt                                                                                   >  1:1544700/1‑149 (MQ=255)
                                                                  ttcggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgc                                                                                                  >  4:560232/1‑133 (MQ=255)
                                                                     ggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtctc                                                                               <  3:685902/149‑1 (MQ=255)
                                                                        ttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtct                                                                                <  4:1597382/145‑1 (MQ=255)
                                                                        ttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtct                                                                                >  3:1597382/1‑145 (MQ=255)
                                                                        ttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                            <  1:1098025/149‑1 (MQ=255)
                                                                          cggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATa                                                                          <  1:10117/149‑1 (MQ=255)
                                                                           ggTGGCACTAATGGTTCTTTGATCTTTAAACAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGGCTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATat                                                                         <  4:610462/149‑2 (MQ=255)
                                                                           ggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATaa                                                                         <  2:1090050/149‑1 (MQ=255)
                                                                           ggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATaa                                                                         <  3:310636/149‑1 (MQ=255)
                                                                                 aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                                   <  2:1092817/149‑1 (MQ=255)
                                                                                 aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                                   <  3:1406468/149‑1 (MQ=255)
                                                                                 aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                                   >  3:286212/1‑149 (MQ=255)
                                                                                 aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                                   >  3:267929/1‑149 (MQ=255)
                                                                                  cTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCg                                                                  >  1:1749120/1‑149 (MQ=255)
                                                                                  cTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCg                                                                  >  1:1051924/1‑149 (MQ=255)
                                                                                   tAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACt                                                                                                                     >  3:756883/1‑97 (MQ=255)
                                                                                      tGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTAt                                                              <  1:1713429/149‑1 (MQ=255)
                                                                                      tGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCAGGCCGTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTAt                                                              <  3:1627500/149‑1 (MQ=255)
                                                                                        caacttagATCTTTTTAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGt                                                            >  2:1687178/8‑149 (MQ=255)
                                                                                           cTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGAt                                                                           <  1:1456267/131‑1 (MQ=255)
                                                                                            tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                        >  3:560571/1‑149 (MQ=255)
                                                                                            tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                        <  4:1556592/149‑1 (MQ=255)
                                                                                            tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                        >  4:663211/1‑149 (MQ=255)
                                                                                            tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                        >  1:1654485/1‑149 (MQ=255)
                                                                                               gATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaat                                                     <  2:1051924/149‑1 (MQ=255)
                                                                                                   tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAAc                                                 >  2:552518/1‑149 (MQ=255)
                                                                                                   tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAAc                                                 <  1:1418714/149‑1 (MQ=255)
                                                                                                   tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACg                                                >  1:336812/1‑149 (MQ=255)
                                                                                                    ttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACg                                                <  3:1174877/149‑1 (MQ=255)
                                                                                                     tAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACgg                                               <  4:1083468/149‑1 (MQ=255)
                                                                                                      aaaaaGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGt                                              <  3:990742/149‑1 (MQ=255)
                                                                                                           gATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATaa                                         >  1:466902/1‑149 (MQ=255)
                                                                                                            aTCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAg                                        <  3:630869/149‑1 (MQ=255)
                                                                                                             tCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGc                                       >  4:1673181/1‑149 (MQ=255)
                                                                                                              cTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGCCTTGTTCCGTAGTTAACATGCTGCCACACTAACAAGTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACTGTCATAAGCa                                      <  3:1032965/149‑1 (MQ=255)
                                                                                                              cTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCa                                      <  1:851790/149‑1 (MQ=255)
                                                                                                               tAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                     >  1:286074/1‑149 (MQ=255)
                                                                                                               tAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                     <  1:1337497/149‑1 (MQ=255)
                                                                                                                 aGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTc                                             >  4:1195315/1‑139 (MQ=255)
                                                                                                                 aGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTc                                             <  3:1195315/139‑1 (MQ=255)
                                                                                                                  attGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTCTCCTTTATTCCGACTTGTTCGGCAGTGAACATGCTCTCACAGTACCAATTTTCTTCTAAGGAGCCGGTCTGTTCTTAATTATCGTTGATAAGCAGGaa                                  >  4:1512866/2‑149 (MQ=255)
                                                                                                                   ttGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTAt                                                              >  3:1076408/1‑120 (MQ=255)
                                                                                                                   ttGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAAt                                 >  1:1345177/1‑149 (MQ=255)
                                                                                                                      tCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGc                              <  2:1544700/149‑1 (MQ=255)
                                                                                                                      tCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGc                              >  3:773612/1‑149 (MQ=255)
                                                                                                                       cATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCt                             >  2:1244140/1‑149 (MQ=255)
                                                                                                                       cATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCt                             <  4:268382/149‑1 (MQ=255)
                                                                                                                       cATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCCGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCt                             <  2:1638092/149‑1 (MQ=255)
                                                                                                                        aTTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTg                            >  4:887775/1‑149 (MQ=255)
                                                                                                                        aTTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTg                            >  3:1138902/1‑149 (MQ=255)
                                                                                                                         ttttCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTgg                           >  1:969489/1‑149 (MQ=255)
                                                                                                                         cattttccccTTATAAAAGGTCCGCCTGCGGGCCTTATTTCTTAGCTTTTATTCCGACTAGTTACGTAGGGAACATGCTGTCACACTAACAATACTCTGATAAGGAGCCGGTATGTTCTTAATATACTGTCAGAAGCAGGAATCGCAgg                           <  4:1072655/142‑1 (MQ=255)
                                                                                                                            tCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAg                        >  4:803035/1‑149 (MQ=255)
                                                                                                                              cACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTGATTAACGGTCATAAGCAGGAATCGCTGGCAGTa                      >  4:1140130/1‑149 (MQ=255)
                                                                                                                              cACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTa                      >  1:1248102/1‑149 (MQ=255)
                                                                                                                              cACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTa                      >  2:102366/1‑149 (MQ=255)
                                                                                                                               aCCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTgg                           >  4:416441/1‑143 (MQ=255)
                                                                                                                               aCCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                     >  3:373943/1‑149 (MQ=255)
                                                                                                                               aCCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                     <  2:1170170/149‑1 (MQ=255)
                                                                                                                                cccTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAg                    >  2:1269134/1‑149 (MQ=255)
                                                                                                                                     ataAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATc               <  3:929055/149‑1 (MQ=255)
                                                                                                                                     ataAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATc               <  1:1292225/149‑1 (MQ=255)
                                                                                                                                     ataAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATc               >  2:386456/1‑149 (MQ=255)
                                                                                                                                      tagaaGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATGCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCg              <  4:1170218/146‑1 (MQ=255)
                                                                                                                                      taAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCg              <  1:769058/149‑1 (MQ=255)
                                                                                                                                         aaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTa                                                                                                            <  2:1478349/52‑1 (MQ=21)
                                                                                                                                         aaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTCACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGc           >  1:1279583/1‑149 (MQ=255)
                                                                                                                                          aGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCa          >  4:386549/1‑149 (MQ=255)
                                                                                                                                            gTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACg                                                <  1:432257/109‑1 (MQ=255)
                                                                                                                                            gTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACg                                                >  2:432257/1‑109 (MQ=255)
                                                                                                                                               cGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATc                                <  2:878298/122‑1 (MQ=255)
                                                                                                                                               cGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCa     <  2:1343764/149‑1 (MQ=255)
                                                                                                                                                  ttGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCAGtt  <  1:103553/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGCAGTT  >  NC_000913/1153067‑1153361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: