Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,135 G→A 16.9% intergenic (+420/+38) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947322,1350GA16.9% 477.5 / 78.1 201intergenic (+420/+38)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (60/107);  new base A (15/19);  total (75/126)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

AAACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321986‑322281
                                                                                                                                                     |                                                                                                                                                  
aaaCACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCgg                                                                                                                                                    <  2:3170096/150‑1 (MQ=255)
aaaCACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCgg                                                                                                                                                    <  2:2432153/150‑1 (MQ=255)
aaaCACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCgg                                                                                                                                                    <  1:5208507/150‑1 (MQ=255)
   cacCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCa                                                                                                                                                 <  1:3122981/150‑1 (MQ=255)
    acCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCAt                                                                                                                                                >  1:3910735/1‑150 (MQ=255)
          gAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                          <  2:2484118/150‑1 (MQ=255)
          gAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                          <  1:2049621/150‑1 (MQ=255)
          gAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                          <  1:937348/150‑1 (MQ=255)
              cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATCCCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAag                                                                                                                                      <  1:4782142/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  2:3647691/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  2:4551991/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  1:4556282/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  1:5137073/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  2:157807/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  1:4016088/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  1:1982867/150‑1 (MQ=255)
                   gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                                 <  1:4839836/150‑1 (MQ=255)
                   gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                                 <  1:4984331/150‑1 (MQ=255)
                   gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                                 >  2:2813177/1‑150 (MQ=255)
                    ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                >  1:1371551/1‑150 (MQ=255)
                    ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAg                                                                                                                                <  2:2960882/150‑1 (MQ=255)
                      tCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGt                                                                                                                              >  1:1530628/1‑150 (MQ=255)
                            cgaccgacCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAggt                                                                                                                        <  1:848290/150‑2 (MQ=255)
                               ccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCTATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAAc                                                                                                                     >  2:3552654/1‑150 (MQ=255)
                                cgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCGCCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAggtgcaa                                                                                                                    <  1:3213519/150‑6 (MQ=255)
                                 gaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAggtgcaat                                                                                                                   <  2:4439663/150‑7 (MQ=255)
                                    ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaa                                                                                                                <  1:2635216/150‑1 (MQ=255)
                                    ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaa                                                                                                                >  2:5181821/1‑150 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGTGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                               <  1:3588968/150‑1 (MQ=255)
                                        gacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCa                                                                                                            <  2:599634/150‑1 (MQ=255)
                                        gacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCa                                                                                                            <  2:3182867/150‑1 (MQ=255)
                                         acgTCTGCTTGCCGGGCGATACCCCCACGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           <  2:141244/150‑1 (MQ=255)
                                         acgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           >  1:5548945/1‑150 (MQ=255)
                                             cTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAAc                                                                                                       <  2:2182657/150‑1 (MQ=255)
                                              tGCTTGCCGGGCGATACCCCCAGGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAgg                                                                                                                         <  1:762375/131‑1 (MQ=255)
                                              tGCTTGCCGGGCGATACCCCCAGGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAgg                                                                                                                         >  2:762375/1‑131 (MQ=255)
                                              tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACt                                                                                                      >  1:385674/1‑150 (MQ=255)
                                              tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACt                                                                                                      <  1:1175174/150‑1 (MQ=255)
                                                cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGt                                                                                                    <  2:5258015/150‑1 (MQ=255)
                                                 ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCTGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                   <  1:4242384/150‑1 (MQ=255)
                                                 ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                   <  2:5301800/150‑1 (MQ=255)
                                                 ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                   <  2:1498415/150‑1 (MQ=255)
                                                 ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                   <  1:2952988/150‑1 (MQ=255)
                                                   gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGc                                                                                                 >  2:1820602/1‑150 (MQ=255)
                                                        gCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGTTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                            <  1:1042888/150‑1 (MQ=255)
                                                        gCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                            <  2:4198697/150‑1 (MQ=255)
                                                         cGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCg                                                                                           >  2:4943462/1‑150 (MQ=255)
                                                            tACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTca                                                                                        <  2:4045021/150‑1 (MQ=255)
                                                            tACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGCGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTca                                                                                        <  2:3403212/150‑1 (MQ=255)
                                                               ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCGCCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                       <  2:3824176/100‑1 (MQ=255)
                                                                  caagttttGATCAGCGAAGCAGAGATGCAGGCGCTGTTGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACcacac                                                                                  <  1:1447879/144‑2 (MQ=255)
                                                                          gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          <  1:4325166/150‑1 (MQ=255)
                                                                          gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          <  1:2369072/150‑1 (MQ=255)
                                                                          gATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          <  2:4135381/150‑1 (MQ=255)
                                                                             cAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGa                                                                       <  1:503850/150‑1 (MQ=255)
                                                                              aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGAt                                                                      >  2:5559909/1‑150 (MQ=255)
                                                                              aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGAt                                                                      <  1:3814571/150‑1 (MQ=255)
                                                                               gCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATc                                                                     <  1:5434263/150‑1 (MQ=255)
                                                                                      agagATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCCGATCGCCCAcc                                                              >  2:3034291/1‑150 (MQ=255)
                                                                                       gagaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTCCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCt                                                             <  2:2112420/150‑1 (MQ=255)
                                                                                        agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTg                                                            >  2:3094655/1‑150 (MQ=255)
                                                                                               ggCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgc                                                     >  2:2493707/1‑150 (MQ=255)
                                                                                                gcgcTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgcg                                                    <  2:1224958/150‑1 (MQ=255)
                                                                                                     gATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGACACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCa                                               <  1:1257375/150‑1 (MQ=255)
                                                                                                        ggCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTg                                            >  2:1728543/1‑150 (MQ=255)
                                                                                                        ggCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTg                                            >  1:2707734/1‑150 (MQ=255)
                                                                                                        ggCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTg                                            <  2:2727460/150‑1 (MQ=255)
                                                                                                           ggTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTc                                         >  1:1804576/1‑150 (MQ=255)
                                                                                                             tgtgGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGc                                       <  1:2408706/150‑1 (MQ=255)
                                                                                                                ggCACTGATCGATTGCACCTGTACCTGGCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGa                                    >  1:464782/1‑150 (MQ=255)
                                                                                                                 gCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCTACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAg                                   >  2:1345930/1‑150 (MQ=255)
                                                                                                                 gCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACaa                                                                                  <  2:1522105/103‑1 (MQ=255)
                                                                                                                 gCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAg                                   <  2:4680887/150‑1 (MQ=255)
                                                                                                                   aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCt                                 <  1:3449386/150‑1 (MQ=255)
                                                                                                                   aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCt                                 >  2:5351810/1‑150 (MQ=255)
                                                                                                                      gATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                              <  1:2472848/150‑1 (MQ=255)
                                                                                                                               cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                     <  2:851282/150‑1 (MQ=255)
                                                                                                                               cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                     <  2:3910735/150‑1 (MQ=255)
                                                                                                                               cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                     <  1:1329190/150‑1 (MQ=255)
                                                                                                                                 ccTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGt                   <  2:3709714/150‑1 (MQ=255)
                                                                                                                                       ccTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa             >  1:1099160/1‑150 (MQ=255)
                                                                                                                                        cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt            >  1:1737081/1‑150 (MQ=255)
                                                                                                                                         tGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg           >  2:3989675/1‑150 (MQ=255)
                                                                                                                                           ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt         <  2:3604948/150‑1 (MQ=255)
                                                                                                                                           ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt         <  2:2058969/150‑1 (MQ=255)
                                                                                                                                           ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt         >  2:1575887/1‑150 (MQ=255)
                                                                                                                                            ctctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        >  1:2427689/1‑150 (MQ=255)
                                                                                                                                             tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCACTGTTGc       >  2:3621267/1‑150 (MQ=255)
                                                                                                                                             tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc       >  2:5530474/1‑150 (MQ=255)
                                                                                                                                             tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           >  1:1516811/1‑50 (MQ=17)
                                                                                                                                             tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc       >  1:619847/1‑150 (MQ=255)
                                                                                                                                             tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc       >  2:3217443/1‑150 (MQ=255)
                                                                                                                                             tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACAACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCACTGTTGc       >  2:3304319/1‑150 (MQ=255)
                                                                                                                                              ctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt      <  1:1728543/150‑1 (MQ=255)
                                                                                                                                              ctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCa              >  2:1090645/1‑142 (MQ=255)
                                                                                                                                               ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc     >  2:1679006/1‑150 (MQ=255)
                                                                                                                                               ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                              >  1:5387356/1‑125 (MQ=255)
                                                                                                                                               ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                              <  2:5387356/125‑1 (MQ=255)
                                                                                                                                               ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc     >  1:1591871/1‑150 (MQ=255)
                                                                                                                                                  gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt  <  1:1904660/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321986‑322281

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: