Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,402,466 T→G 5.2% M264L (ATG→CTG)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,402,4660TG 5.2% 1664.9 / 35.2 512M264L (ATG→CTG) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base T (231/247);  new base G (2/32);  total (233/279)
Fisher's exact test for biased strand distribution p-value = 3.24e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 527 total aligned reads displayed.

GGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAACGTTA  >  NC_000913/1402317‑1402615
                                                                                                                                                     |                                                                                                                                                     
ggCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAt                                                                                                                                                       >  1:8147721/1‑150 (MQ=255)
  ccTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAttt                                                                                                                                                     >  1:6472840/1‑150 (MQ=255)
  ccTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAttt                                                                                                                                                     >  1:4240407/1‑150 (MQ=255)
    tGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCg                                                                                                                                                   >  1:6162652/1‑150 (MQ=255)
     ggTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCgg                                                                                                                                                  <  2:7542774/150‑1 (MQ=255)
      gTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCggg                                                                                                                                                 <  2:4817189/150‑1 (MQ=255)
      gTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCggg                                                                                                                                                 >  1:7620282/1‑150 (MQ=255)
             gCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCATTGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAGATATGCTGCACGTCGGTCATTTCGGGGGCGCgg                                                                                                                                          >  2:849646/1‑150 (MQ=255)
             gCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCgg                                                                                                                                          >  1:5384860/1‑150 (MQ=255)
             gCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCgg                                                                                                                                          <  1:3094946/150‑1 (MQ=255)
               aTTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgat                                                                                                                                        <  1:4005736/150‑1 (MQ=255)
                ttttCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgata                                                                                                                                       <  1:5061215/150‑1 (MQ=255)
                   tCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgataaga                                                                                                                                    <  1:7337781/150‑1 (MQ=255)
                   tCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgataaga                                                                                                                                    >  2:5023856/1‑150 (MQ=255)
                      aaGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgataagataa                                                                                                                                 <  1:2846807/150‑1 (MQ=255)
                      aaGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGgataagataa                                                                                                                                 <  2:7914585/150‑1 (MQ=255)
                       aGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATaaa                                                                                                                                <  1:3859362/150‑1 (MQ=255)
                       aGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATaaa                                                                                                                                <  1:8736476/150‑1 (MQ=255)
                        ggTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAg                                                                                                                               <  1:1913282/150‑1 (MQ=255)
                           gCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAAAATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCAc                                                                                                                            >  2:172539/1‑150 (MQ=255)
                             gATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACtt                                                                                                                          >  2:5213722/1‑150 (MQ=255)
                             gATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGGCAGGGCGGGGGCGCGGATAAGATAAAGCACtt                                                                                                                          <  1:6814393/150‑1 (MQ=255)
                               ttCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTct                                                                                                                        >  1:1971442/1‑150 (MQ=255)
                                 cGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTctgc                                                                                                                      <  2:7636311/150‑1 (MQ=255)
                                 cGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTctgc                                                                                                                      <  2:4091320/150‑1 (MQ=255)
                                  gggAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTctgcc                                                                                                                     <  1:7156388/150‑1 (MQ=255)
                                    gagtTAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTAGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAGTTCGGGGGCGCGGATAAGATAAAGCACTTctgcctg                                                                                                                   <  1:6006874/150‑1 (MQ=255)
                                          aCTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGAGTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTACTGCCTGCGCCTg                                                                                                             >  2:9368794/1‑150 (MQ=255)
                                           cTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTgg                                                                                                            >  2:2743227/1‑150 (MQ=255)
                                           cTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTgg                                                                                                            >  1:1069119/1‑150 (MQ=255)
                                            tGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGa                                                                                                           <  1:9285319/150‑1 (MQ=255)
                                            tGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAGTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGa                                                                                                           <  1:4866484/150‑1 (MQ=255)
                                            tGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAGTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGa                                                                                                           >  2:5503452/1‑150 (MQ=255)
                                             ggAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGAc                                                                                                          <  2:9516604/150‑1 (MQ=255)
                                              gAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGAcc                                                                                                         >  1:2295120/1‑150 (MQ=255)
                                               aaCAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCa                                                                                                        >  2:9885824/1‑150 (MQ=255)
                                               aaCAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGGCAGGGCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCa                                                                                                        <  1:10240559/150‑1 (MQ=255)
                                                 cAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACg                                                                                                      <  2:2434425/150‑1 (MQ=255)
                                                  aGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGt                                                                                                     >  2:1087663/1‑150 (MQ=255)
                                                    gCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTg                                                                                                   >  1:8371540/1‑150 (MQ=255)
                                                     cTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTgg                                                                                                  >  1:5011299/1‑150 (MQ=255)
                                                         gTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCga                                                                                              <  1:10135156/150‑1 (MQ=255)
                                                          tCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCgag                                                                                             >  2:6295286/1‑150 (MQ=255)
                                                            gAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATCAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGAt                                                                                           >  2:1341066/1‑150 (MQ=255)
                                                                 ggCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCccgc                                                                                      >  1:9783485/1‑150 (MQ=255)
                                                                   cATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCccgccg                                                                                    <  2:2541029/150‑1 (MQ=255)
                                                                    aTTCTACCGTTGTTTCGGTCATCTATGCCGCACCTTCGTCGATTTTGGCGACCCGATCATAAATATGCTTCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCcgccgc                                                                                   <  1:5492860/150‑1 (MQ=255)
                                                                    aTTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCcgccgc                                                                                   >  1:2536598/1‑150 (MQ=255)
                                                                    aTTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCcgccgc                                                                                   >  2:78289/1‑150 (MQ=255)
                                                                     ttCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCt                                                                                  >  2:6902060/1‑150 (MQ=255)
                                                                        aaCCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTAtt                                                                               >  2:3326939/1‑150 (MQ=255)
                                                                         aCCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTAttt                                                                              <  2:7501439/150‑1 (MQ=255)
                                                                          ccGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTtg                                                                             >  2:7964172/1‑150 (MQ=255)
                                                                           cGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTtgt                                                                            <  2:5675326/150‑1 (MQ=255)
                                                                           cGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTtgt                                                                            >  2:378357/1‑150 (MQ=255)
                                                                             tGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGa                                                                          <  1:9124302/150‑1 (MQ=255)
                                                                             tGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGa                                                                          <  2:2252611/150‑1 (MQ=255)
                                                                               gTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATg                                                                        <  1:6922229/150‑1 (MQ=255)
                                                                                tttCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATgg                                                                       >  2:6062846/1‑150 (MQ=255)
                                                                                 ttCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGc                                                                      <  2:9075027/150‑1 (MQ=255)
                                                                                    ggTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATa                                                                   >  1:2143700/1‑150 (MQ=255)
                                                                                         tcaATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCCTAGTGTa                                                              >  1:4470543/1‑150 (MQ=255)
                                                                                           aaTGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTAcg                                                            >  2:975650/1‑150 (MQ=255)
                                                                                           aaTGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTAcg                                                            >  1:1422296/1‑150 (MQ=255)
                                                                                               ccGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTAcgcgcg                                                        <  1:7066983/150‑1 (MQ=255)
                                                                                                    ccTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCtttt                                                   <  2:9600352/150‑1 (MQ=255)
                                                                                                     cTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCttttt                                                  <  2:1881872/150‑1 (MQ=255)
                                                                                                     cTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCAGTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGTATAGTGTACGCGCGCttttt                                                  <  1:2708480/150‑1 (MQ=255)
                                                                                                        cGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCaa                                               <  2:6472840/150‑1 (MQ=255)
                                                                                                          gCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCaata                                             >  2:8636626/1‑150 (MQ=255)
                                                                                                               tttGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATg                                        <  2:4531610/150‑1 (MQ=255)
                                                                                                                   gCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCg                                    >  1:1678436/1‑150 (MQ=255)
                                                                                                                    cGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGt                                   >  2:1022434/1‑150 (MQ=255)
                                                                                                                      aCCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTg                                 >  2:126373/1‑150 (MQ=255)
                                                                                                                       cccGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGa                                <  1:3744538/150‑1 (MQ=255)
                                                                                                                       cccGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGa                                <  1:3415909/150‑1 (MQ=255)
                                                                                                                       cccGATCATAAATATGCTGCACGTCGGTCAGGTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGa                                <  1:7832826/150‑1 (MQ=255)
                                                                                                                        ccGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATACGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTCCGCGCGCTTTTTCAATAATATGTTCGTTGAg                               >  1:5252480/1‑150 (MQ=255)
                                                                                                                        ccGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAg                               >  1:5774569/1‑150 (MQ=255)
                                                                                                                         cGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAgg                              <  2:8371540/150‑1 (MQ=255)
                                                                                                                           aTCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGaa                            >  2:2863998/1‑150 (MQ=255)
                                                                                                                            tCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAg                           <  2:3661882/150‑1 (MQ=255)
                                                                                                                              aTAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGtt                         <  1:3232004/150‑1 (MQ=255)
                                                                                                                              aTAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGtt                         <  2:8860254/150‑1 (MQ=255)
                                                                                                                              aTAAATATGCTGCACGTCGGTCAGTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGtt                         <  1:5691732/150‑1 (MQ=255)
                                                                                                                               tAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTg                        <  1:940144/150‑1 (MQ=255)
                                                                                                                               tAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTg                        <  2:10038034/150‑1 (MQ=255)
                                                                                                                                 aaTATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGt                      <  2:5384860/150‑1 (MQ=255)
                                                                                                                                 aaTATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGt                      <  1:9883386/150‑1 (MQ=255)
                                                                                                                                 aaTATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGt                      <  2:4240407/150‑1 (MQ=255)
                                                                                                                                   gatGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGc                    >  2:7069120/2‑150 (MQ=255)
                                                                                                                                        tgcACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTg               >  1:8663323/1‑150 (MQ=255)
                                                                                                                                        tgcACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTg               <  1:9240369/150‑1 (MQ=255)
                                                                                                                                             gTCGGTCAGGTCGGGGGCGCTGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGtcat          <  1:4624520/150‑1 (MQ=255)
                                                                                                                                             gTCGGTCAGGTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGtcat          <  1:1194659/150‑1 (MQ=255)
                                                                                                                                                ggTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCaa       <  2:3396570/150‑1 (MQ=255)
                                                                                                                                                 gTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAAc      >  1:1788722/1‑150 (MQ=255)
                                                                                                                                                 gTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAAc      <  2:9025369/150‑1 (MQ=255)
                                                                                                                                                 gTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAAc      <  2:2085935/150‑1 (MQ=255)
                                                                                                                                                     tttCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAACGTTa  >  2:5835115/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAATGCCGCACCTTCGGCGATTTTGGCGACCCGATCATAAATATGCTGCACGTCGGTCATTTCGGGGGCGCGGATAAGATAAAGCACTTCTGCCTGCGCCTGGACCACGTTGGGCGAGATCCCGCCGCTATTTGTGATGGCATAGTGTACGCGCGCTTTTTCAATAATATGTTCGTTGAGGAAGTTGGTGCCAGTGGTCATCAACGTTA  >  NC_000913/1402317‑1402615

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