Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,142 T→C 17.0% intergenic (+427/+31) 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,1420TC17.0% 234.5 / 37.1 106intergenic (+427/+31)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base C (6/12);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 111 total aligned reads displayed.

ATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGCA  >  NC_002947/321993‑322291
                                                                                                                                                       |                                                                                                                                                      
     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:127652/150‑1 (MQ=255)
        aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                               <  1:1602805/150‑1 (MQ=255)
        aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                               <  2:754404/150‑1 (MQ=255)
        aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                               <  2:421450/150‑1 (MQ=255)
         tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAagag                                                                                                                                              <  1:1286055/150‑1 (MQ=255)
         tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGAGCATGCCCGCGAagag                                                                                                                                              <  1:1315488/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)
                  ggACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                                     >  1:1711590/1‑150 (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                                                                                                                                                 <  2:2181398/115‑1 (MQ=255)
                                         cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                 >  1:2181398/1‑115 (MQ=255)
                                          ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAATTCAGAAACTGTa                                                                                                             <  1:859022/150‑1 (MQ=255)
                                              cGGGCTATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                                    <  2:2354409/123‑1 (MQ=255)
                                              cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                               <  1:33646/112‑1 (MQ=255)
                                              cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                                    >  1:2354409/1‑123 (MQ=255)
                                              cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                            >  2:1039603/1‑147 (MQ=255)
                                              cGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                            <  1:1039603/147‑1 (MQ=255)
                                                     tACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTca                                                                                                  <  2:841746/150‑1 (MQ=255)
                                                      aCCCCCAAGGGTTGATCAGCGAAGCAGAGATGCCGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑GTGTACCTTCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGAAGCTACCCGTcac                                                                                                 <  2:471348/150‑1 (MQ=255)
                                                       cccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacc                                                                                                >  2:826681/1‑150 (MQ=255)
                                                        ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTTTTCGCGGGCATGCTCGCTAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                               >  1:282438/1‑150 (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)
                                                                                      cAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCt                                                                 <  2:485860/150‑1 (MQ=255)
                                                                                      cAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCt                                                                 <  2:1700428/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                                                                                                                                     <  1:474405/72‑1 (MQ=255)
                                                                                                tGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                                     >  2:474405/1‑72 (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)
                                                                                                              tGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgt                                         >  2:1561371/1‑150 (MQ=255)
                                                                                                               gATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                                      <  1:813969/88‑1 (MQ=255)
                                                                                                               gATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAccc                                                                                                      >  2:813969/1‑88 (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)
                                                                                                                                        tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc                 >  2:2147408/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                                    <  2:1356034/83‑1 (MQ=255)
                                                                                                                                        tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGa                                              >  1:1154624/1‑121 (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                                                                                                                     >  1:482952/1‑48 (MQ=14)
                                                                                                                                          ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                     <  2:482952/48‑1 (MQ=14)
                                                                                                                                          ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                     <  2:488156/48‑1 (MQ=14)
                                                                                                                                          ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                     >  1:488156/1‑48 (MQ=14)
                                                                                                                                          ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTc               >  1:1227797/1‑150 (MQ=255)
                                                                                                                                            cgcgAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCaa             <  1:528006/150‑1 (MQ=255)
                                                                                                                                             gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                                        <  1:590921/58‑1 (MQ=255)
                                                                                                                                             gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                                        >  2:590921/1‑58 (MQ=255)
                                                                                                                                             gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGATGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt            <  2:666140/150‑1 (MQ=255)
                                                                                                                                                   aTGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccgc                                                          >  2:526535/1‑98 (MQ=255)
                                                                                                                                                   aTGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccgc                                                          <  1:526535/98‑1 (MQ=255)
                                                                                                                                                      cTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt                <  2:612281/137‑1 (MQ=255)
                                                                                                                                                      cTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt                >  1:612281/1‑137 (MQ=255)
                                                                                                                                                       tCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc                 >  1:1834626/1‑135 (MQ=255)
                                                                                                                                                       tCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc                 <  2:1834626/135‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
ATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCAC‑‑CTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGCA  >  NC_002947/321993‑322291

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