Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,480,810 G→C 10.3% intergenic (+98/+347) PP_2172 → / ← PP_2173 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,480,8100GC10.3% 521.9 / 39.0 166intergenic (+98/+347)PP_2172/PP_2173hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (72/76);  new base C (12/5);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 170 total aligned reads displayed.

TTCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGATCA  >  NC_002947/2480661‑2480959
                                                                                                                                                     |                                                                                                                                                     
ttCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTg                                                                                                                                                       >  2:255607/1‑150 (MQ=255)
ttCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTg                                                                                                                                                       >  2:2085405/1‑150 (MQ=255)
   gAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGc                                                                                                                                                    >  1:4875233/1‑150 (MQ=255)
         ccTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGc                                                                                                                                              >  2:5031957/1‑150 (MQ=255)
           tCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAg                                                                                                                                            >  1:1073844/1‑150 (MQ=255)
           tCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAg                                                                                                                                            >  1:3814408/1‑150 (MQ=255)
              gCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGgt                                                                                                                                         <  1:3630944/150‑1 (MQ=255)
              gCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt                                                                                                                                         <  1:1636940/150‑1 (MQ=255)
               cAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGgtg                                                                                                                                        >  1:4111159/1‑150 (MQ=255)
                aaCGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGgtgt                                                                                                                                       >  1:3785759/1‑150 (MQ=255)
                 aCGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGtgtt                                                                                                                                      <  1:3945315/150‑1 (MQ=255)
                  cGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGtgttt                                                                                                                                     >  1:4759417/1‑150 (MQ=255)
                     tGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgtgtggggg                                                                                                                                  <  2:3634581/150‑6 (MQ=39)
                         gggTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGttttgttttg                                                                                                                              <  2:319818/150‑1 (MQ=255)
                         gggTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGttttgttttg                                                                                                                              <  2:4149526/150‑1 (MQ=255)
                          ggTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTgg                                                                                                                             <  1:2081341/150‑1 (MQ=255)
                           gTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTggg                                                                                                                            <  1:255607/150‑2 (MQ=255)
                           gTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga                                                                                                                            >  1:316459/1‑150 (MQ=255)
                           gTCTTCGTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga                                                                                                                            <  2:2183586/150‑1 (MQ=255)
                               tcttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGAcc                                                                                                                        <  1:1123032/150‑1 (MQ=255)
                                 ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcat                                                                                                                      >  1:4124434/1‑150 (MQ=255)
                                 ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcag                                                                                                                      >  2:3233594/1‑149 (MQ=255)
                                  tcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgtgtgggg                                                                                                                                   <  1:302112/136‑5 (MQ=38)
                                    aaaaaaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGgtgtggggg                                                                                                                                  <  1:2005666/135‑6 (MQ=255)
                                    aaaaaaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGttttgttttg                                                                                                                              <  1:1642853/139‑1 (MQ=255)
                                    aaaaaaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcatcat                                                                                                                   >  2:3692240/1‑150 (MQ=255)
                                         aTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcga                                                                                                              <  2:2996241/150‑1 (MQ=255)
                                           gATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                            <  2:2070970/150‑1 (MQ=255)
                                             tCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagcga                                                                                                          <  1:2877975/150‑1 (MQ=255)
                                              cGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATAAcgagcgag                                                                                                         >  1:5225102/1‑150 (MQ=255)
                                                  gatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAg                                                                                                     <  1:897887/150‑1 (MQ=255)
                                                   atTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGc                                                                                                    >  2:3761014/1‑150 (MQ=255)
                                                        aCCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTg                                                                                               <  2:777781/150‑1 (MQ=255)
                                                         ccGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTgg                                                                                              >  2:2002912/1‑150 (MQ=255)
                                                            cATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCgg                                                                                           <  1:2473920/150‑1 (MQ=255)
                                                               ccATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgat                                                                                        <  2:1311194/150‑1 (MQ=255)
                                                                cATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgctc                                                                                       >  1:3212131/1‑150 (MQ=255)
                                                                cATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgatc                                                                                       <  2:317575/150‑1 (MQ=255)
                                                                cATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgatc                                                                                       <  1:669813/150‑1 (MQ=255)
                                                                 aTGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgatcg                                                                                      >  2:2826838/1‑150 (MQ=255)
                                                                      gTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGAtctc                                                                                 >  2:3234816/1‑150 (MQ=255)
                                                                      gTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGAtctc                                                                                 >  1:1523853/1‑150 (MQ=255)
                                                                       tACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCg                                                                                <  2:483766/150‑1 (MQ=255)
                                                                         ccGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGt                                                                              >  1:2267964/1‑150 (MQ=255)
                                                                           ggCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTaa                                                                            <  1:1432997/150‑1 (MQ=255)
                                                                           ggCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTaa                                                                            >  1:1931766/1‑150 (MQ=255)
                                                                            gCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAAc                                                                           >  2:1465492/1‑150 (MQ=255)
                                                                            gCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAAc                                                                           <  2:3086515/150‑1 (MQ=255)
                                                                              tGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACca                                                                         <  1:4290999/150‑1 (MQ=255)
                                                                              tGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTATGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACca                                                                         >  2:5046815/1‑150 (MQ=255)
                                                                                tAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagcga                                                                                                          >  2:68406/1‑115 (MQ=255)
                                                                                tAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagcga                                                                                                          <  1:68406/115‑1 (MQ=255)
                                                                                        aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTc                                                               >  2:5424558/1‑150 (MQ=255)
                                                                                        aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTc                                                               <  2:815246/150‑1 (MQ=255)
                                                                                         cGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTcc                                                              <  1:3397893/150‑1 (MQ=255)
                                                                                          gagaGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTc                                                               <  2:2034881/148‑1 (MQ=255)
                                                                                          gagaGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGGTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCg                                                             >  1:5012487/1‑150 (MQ=255)
                                                                                            gagCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGtt                                                                    <  1:4056691/141‑1 (MQ=255)
                                                                                            gagCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGtt                                                                    >  2:4056691/1‑141 (MQ=255)
                                                                                                        ggCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTa                                               <  2:1561943/150‑1 (MQ=255)
                                                                                                         gCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTAc                                              >  1:4066220/1‑150 (MQ=255)
                                                                                                          ccGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTAGCAGGGTGTTTTGTTTTGTAGACAATAATAAAGAGCGAGTCAGCACGTGGCGGGATCGATCTCTGTAACCAAAGATAACAACGGCAAGGCCAAGGTACg                                             >  1:3998232/1‑150 (MQ=255)
                                                                                                           cGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCAAGTGGCGGGATCGATCTCGGTAACCTCAGTTACCTCCGGCATGGCCAAGGTACGt                                            <  1:274570/150‑1 (MQ=255)
                                                                                                             tGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTg                                          <  1:3122111/150‑1 (MQ=255)
                                                                                                                tAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTc                                       <  1:4660223/150‑1 (MQ=255)
                                                                                                                   ttGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCCGGGGGGTTTGTTTTGGAGACCCTCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCAAGGCCAAGGTACGTTGGTCCaa                                    >  1:5112803/1‑150 (MQ=255)
                                                                                                                   ttGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCaa                                    <  1:1516071/150‑1 (MQ=255)
                                                                                                                    tGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAg                                   >  1:4904695/1‑150 (MQ=255)
                                                                                                                     gAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGc                                  >  2:3515311/1‑150 (MQ=255)
                                                                                                                     gAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGc                                  >  1:1699315/1‑150 (MQ=255)
                                                                                                                     gAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGc                                  >  1:5180068/1‑150 (MQ=255)
                                                                                                                       aCGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCgg                                <  1:805809/150‑1 (MQ=255)
                                                                                                                        cGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGgt                               <  1:2092276/150‑1 (MQ=255)
                                                                                                                         gggTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGgtg                              >  1:2819743/1‑150 (MQ=255)
                                                                                                                           gTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGtt                            >  2:1240598/1‑150 (MQ=255)
                                                                                                                            ttACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGttt                           <  1:510194/150‑1 (MQ=255)
                                                                                                                              aaacaaaacacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGttt                           >  1:4313452/8‑148 (MQ=38)
                                                                                                                              aaacaaaacacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGttt                           <  2:4313452/141‑1 (MQ=38)
                                                                                                                              aCCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCg                         <  2:1899259/150‑1 (MQ=255)
                                                                                                                               ccccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGc                        <  2:1802269/150‑1 (MQ=255)
                                                                                                                                cccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCa                       >  2:1089760/1‑150 (MQ=255)
                                                                                                                                 ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTc                          >  2:4658418/1‑146 (MQ=255)
                                                                                                                                 ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTc                          <  1:4658418/146‑1 (MQ=255)
                                                                                                                                  cccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag                     <  1:83621/150‑1 (MQ=255)
                                                                                                                                  aaaacacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag                     >  1:3109398/4‑150 (MQ=39)
                                                                                                                                    cacaCCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaagaa                   <  2:3769532/150‑1 (MQ=255)
                                                                                                                                    cacaCCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaagaa                   <  1:1849817/150‑1 (MQ=255)
                                                                                                                                    cacaCCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaagaa                   <  2:1699315/150‑1 (MQ=255)
                                                                                                                                      cacCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAgg                 >  2:892895/1‑150 (MQ=255)
                                                                                                                                       acCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGt                >  2:3811871/1‑150 (MQ=255)
                                                                                                                                       acCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGt                >  2:4314306/1‑150 (MQ=255)
                                                                                                                                         ccTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCCTCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTCCGTTGGTCCACGCGGTGTTTCGCAAGAAGGTAg              >  2:530647/1‑150 (MQ=255)
                                                                                                                                         ccTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAg              >  2:236346/1‑150 (MQ=255)
                                                                                                                                         ccTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAg              >  1:2893728/1‑150 (MQ=255)
                                                                                                                                          cTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGa             >  1:2475678/1‑150 (MQ=255)
                                                                                                                                          cTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGa             >  1:3850823/1‑150 (MQ=255)
                                                                                                                                            gCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACg           <  1:1813586/150‑1 (MQ=255)
                                                                                                                                                gcgcTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCAc       <  2:5111884/150‑1 (MQ=255)
                                                                                                                                                    tGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGGGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGATc   >  2:1951301/1‑150 (MQ=255)
                                                                                                                                                     gAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGATCa  <  2:5185444/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGATCA  >  NC_002947/2480661‑2480959

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: