Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,135 G→A 18.3% 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,1350GA18.3% 216.2 / 41.1 104intergenic (+420/+38)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base A (6/13);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 107 total aligned reads displayed.

AAACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321986‑322281
                                                                                                                                                       |                                                                                                                                                   
aaaCACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCgg                                                                                                                                                     >  2:1634278/1‑150 (MQ=255)
   cacCATGGAAGCATCACGGTCAGGACGACCGACCCACCCGTTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCa                                                                                                                                                  <  1:723504/150‑1 (MQ=255)
    acCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCAt                                                                                                                                                 <  1:1658959/150‑1 (MQ=255)
       aTGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCt                                                                                                                                              >  2:2227237/1‑150 (MQ=255)
            aGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGa                                                                                                                                         >  1:485860/1‑150 (MQ=255)
            aGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGa                                                                                                                                         >  1:1700428/1‑150 (MQ=255)
             gCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                       <  2:497778/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  1:438403/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  1:1602805/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  1:127652/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  2:421450/150‑1 (MQ=255)
               aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  2:754404/150‑1 (MQ=255)
                tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAagag                                                                                                                                    <  1:1315488/150‑1 (MQ=255)
                tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAagag                                                                                                                                    <  1:1286055/150‑1 (MQ=255)
                  aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGc                                                                                                                                  <  2:1575156/150‑1 (MQ=255)
                    ggTCAGGACGCCCGACCCACCACGTCTGCTTGCCGGGCGATGCCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                <  1:21401/150‑1 (MQ=255)
                    ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                <  1:2227237/150‑1 (MQ=255)
                    ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                <  2:28662/150‑1 (MQ=255)
                      tCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                           >  2:2289294/1‑138 (MQ=255)
                      tCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                           <  1:2289294/138‑1 (MQ=255)
                              accgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGcca                                                                                                                      >  1:2181408/1‑148 (MQ=255)
                              accgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCAAAGTGGCACTGATCCAATGCGTTGTTGTGCCTGAGTC‑‑CTCGGGCATGCTCGCGAAGAGGCAGGTACAGGCaa                                                                                                                      <  1:415526/150‑1 (MQ=255)
                                cgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACt                                                                                                                    >  1:1724934/1‑150 (MQ=255)
                                cgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACt                                                                                                                    >  2:2233933/1‑150 (MQ=255)
                                  accCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTa                                                                                                                  <  2:610927/150‑1 (MQ=255)
                                    ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaa                                                                                                                >  1:330588/1‑150 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                                   <  1:2052424/130‑1 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                                   >  2:2052424/1‑130 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                               <  2:326926/150‑1 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                               >  1:326926/1‑150 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTAGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                               <  1:1245549/150‑1 (MQ=255)
                                     cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCAAGATGAGATATAGGCAGGTACAGGCAACTTaaaa                                                                                                               <  1:1881572/150‑1 (MQ=255)
                                      acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAAt                                                                                                              >  2:1773333/1‑150 (MQ=255)
                                        gacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCa                                                                                                            <  2:1711590/150‑1 (MQ=255)
                                        gacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCa                                                                                                            <  2:1850794/150‑1 (MQ=255)
                                                cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                       >  1:2181398/1‑115 (MQ=255)
                                                cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                       <  2:2181398/115‑1 (MQ=255)
                                                 ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAATTCAGAAACTGTa                                                                                                   <  1:859022/150‑1 (MQ=255)
                                                     cGGGCTATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                          <  2:2354409/123‑1 (MQ=255)
                                                     cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     >  2:33646/1‑112 (MQ=255)
                                                     cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  1:33646/112‑1 (MQ=255)
                                                     cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                          >  1:2354409/1‑123 (MQ=255)
                                                     cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                  <  1:1039603/147‑1 (MQ=255)
                                                     cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                  >  2:1039603/1‑147 (MQ=255)
                                                            tACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTca                                                                                        <  2:841746/150‑1 (MQ=255)
                                                              cccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacc                                                                                      >  2:826681/1‑150 (MQ=255)
                                                               ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTTTTCGCGGGCATGCTCGCTAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                     >  1:282438/1‑150 (MQ=255)
                                                                cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCgg                                                                                                                                                     <  1:2122590/86‑1 (MQ=255)
                                                                cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCgg                                                                                                                                                     >  2:2122590/1‑86 (MQ=255)
                                                                cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACGGGCAACTTAAAATCAGAAACTGTAGCTACCCGGcaccac                                                                                    >  1:1319250/1‑150 (MQ=255)
                                                                cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGACACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                >  1:1849779/1‑106 (MQ=255)
                                                                cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGACACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                <  2:1849779/106‑1 (MQ=255)
                                                                 ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACcaca                                                                                   >  1:1649012/1‑150 (MQ=255)
                                                                  cAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGATACTGTAGCTACCCGTCACCACaa                                                                                  <  2:1185638/150‑1 (MQ=255)
                                                                   aaGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAg                                                                                 <  2:242978/150‑1 (MQ=255)
                                                                      ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCCCTTGCCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCt                                                                              <  1:1967758/150‑1 (MQ=255)
                                                                      ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAAc                                                                                                                     <  2:2296803/111‑1 (MQ=255)
                                                                               gCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGCTCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                <  2:1331427/91‑1 (MQ=255)
                                                                               gCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                                >  1:1331427/1‑91 (MQ=255)
                                                                                  aaGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGcc                                                                  >  1:1141238/1‑150 (MQ=255)
                                                                                     cAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                               >  1:66940/1‑150 (MQ=255)
                                                                                        agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCt                                                                              <  1:173894/132‑1 (MQ=255)
                                                                                        agaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCt                                                                              >  2:173894/1‑132 (MQ=255)
                                                                                             cAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCt                                                       <  2:1700428/150‑1 (MQ=255)
                                                                                             cAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCt                                                       <  2:485860/150‑1 (MQ=255)
                                                                                               ggCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgc                                                     <  2:1057065/150‑1 (MQ=255)
                                                                                                 cgcTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgcgc                                                   <  2:1724934/150‑1 (MQ=255)
                                                                                                  gcTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                             >  2:850228/1‑60 (MQ=255)
                                                                                                  gcTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                             <  1:850228/60‑1 (MQ=255)
                                                                                                       tGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                           >  2:474405/1‑72 (MQ=255)
                                                                                                       tGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                           <  1:474405/72‑1 (MQ=255)
                                                                                                              gtgGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGcc                                      <  2:2062159/150‑1 (MQ=255)
                                                                                                                   aCTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCt                                 >  1:72885/1‑150 (MQ=255)
                                                                                                                     tGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGCTCCGCGCCGCACTGGTCACCGAGCTgt                               >  1:411247/1‑150 (MQ=255)
                                                                                                                      gATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                            <  1:813969/88‑1 (MQ=255)
                                                                                                                      gATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                            >  2:813969/1‑88 (MQ=255)
                                                                                                                          gATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAAATGTAGGTACCCGTCACCAAAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                          >  2:566245/1‑150 (MQ=255)
                                                                                                                                    cTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc                  >  2:1746326/1‑150 (MQ=255)
                                                                                                                                          cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATTAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt            >  1:768793/1‑150 (MQ=255)
                                                                                                                                          cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt            >  2:1643820/1‑150 (MQ=255)
                                                                                                                                           tGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAtt           >  2:2401234/1‑149 (MQ=255)
                                                                                                                                             ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                 >  1:2346335/1‑142 (MQ=255)
                                                                                                                                             ccTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                 <  2:2346335/142‑1 (MQ=255)
                                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          >  1:1356034/1‑83 (MQ=255)
                                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          <  2:1356034/83‑1 (MQ=255)
                                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGa                                    >  1:1154624/1‑121 (MQ=255)
                                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGa                                    <  2:1154624/121‑1 (MQ=255)
                                                                                                                                                ctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt      <  1:878214/150‑1 (MQ=255)
                                                                                                                                                ctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGCCGCCTAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt      <  1:1773333/150‑1 (MQ=255)
                                                                                                                                                 ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        <  2:376400/147‑1 (MQ=255)
                                                                                                                                                 ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        >  1:376400/1‑147 (MQ=255)
                                                                                                                                                 ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc     <  1:2147408/150‑1 (MQ=255)
                                                                                                                                                 ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           <  2:488156/48‑1 (MQ=14)
                                                                                                                                                 ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           >  1:488156/1‑48 (MQ=14)
                                                                                                                                                 ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           <  2:482952/48‑1 (MQ=14)
                                                                                                                                                 ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                           >  1:482952/1‑48 (MQ=14)
                                                                                                                                                 ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc     >  1:1227797/1‑150 (MQ=255)
                                                                                                                                                    gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                              <  1:590921/58‑1 (MQ=255)
                                                                                                                                                    gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                              >  2:590921/1‑58 (MQ=255)
                                                                                                                                                    gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGATGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt  <  2:666140/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
AAACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321986‑322281

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