Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,341,128 T→A 6.5% intergenic (‑165/+78) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,341,1280TA 6.5% 482.0 / 3.2 163intergenic (‑165/+78)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base T (77/70);  new base A (11/2);  total (89/74)
Fisher's exact test for biased strand distribution p-value = 3.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCC  >  NC_002947/3340980‑3341275
                                                                                                                                                    |                                                                                                                                                   
 ggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATAc                                                                                                                                                   <  2:615527/150‑1 (MQ=255)
   cggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCt                                                                                                                                                 <  1:27991/150‑1 (MQ=255)
      ttGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACGAGCGTCCAATACCtgtg                                                                                                                                              <  1:2739250/150‑1 (MQ=255)
        gggCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCg                                                                                                                                            <  1:1733521/150‑1 (MQ=255)
            gacttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCCCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTcc                                                                                                                                        >  2:1811360‑M1/91‑150 (MQ=255)
             tGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTccc                                                                                                                                       <  2:4302240/150‑1 (MQ=255)
              cttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTC‑‑‑‑‑‑‑‑‑‑CCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                                                            >  1:690250‑M1/89‑150 (MQ=255)
                tCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgt                                                                                                                                    <  2:2956061/150‑1 (MQ=255)
                tCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCCCTCACCCAAAGACCCACAACCGGCCAAAACCAGTTCGGTCCCagt                                                                                                                                    >  1:2274043/1‑150 (MQ=255)
                  ccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTgg                                                                                                                                  <  2:1560961/150‑1 (MQ=255)
                  ccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCAGTgg                                                                                                                                  >  2:645754/1‑150 (MQ=255)
                    atcatGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGGTTTGCCGCTGTGCCAAGGCCCCCCAAGCGTGAAATAACTTTTCGGTGGCTGTGGgt                                                                                                                                <  1:3549361‑M2/150‑3 (MQ=255)
                     tcatGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTGGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgtcgatc                                                                                                                               <  2:3047006/150‑6 (MQ=255)
                       atGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAAAACCAGTGCGGTCCCGGTGGgagcg                                                                                                                             <  1:4057802‑M2/150‑6 (MQ=255)
                        tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgg                                                                                                                            >  2:1483627‑M2/1‑144 (MQ=255)
                        tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgg                                                                                                                            >  1:4002699‑M2/1‑144 (MQ=255)
                           tacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCa                                                                                                                         >  1:2498630‑M1/76‑150 (MQ=255)
                           ccGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCCCGCTCCCCACGCCCCACCAGCGGCCACTACCAGTGCGGTCCCGGTGGgagcgggcc                                                                                                                         <  2:3431450‑M2/150‑10 (MQ=255)
                           ccGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGAGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggca                                                                                                                         >  2:2582425‑M2/1‑141 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaa                                                                                                                        >  1:731757‑M2/1‑140 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaa                                                                                                                        >  1:1084543‑M2/1‑140 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaa                                                                                                                        >  1:164817‑M2/1‑140 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaa                                                                                                                        >  1:2961661‑M2/1‑140 (MQ=255)
                            acgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                                                        >  1:3024314‑M1/75‑150 (MQ=255)
                                  gcgAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCGCTCTCCCCAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgc                                                                                                                  <  2:2280031‑M2/150‑17 (MQ=255)
                                     aaGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaa                                                                                                               >  2:3713921‑M2/1‑131 (MQ=255)
                                      aGGACAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGTTCGGTGGGATGGGGATGTACTCTTCGAGGGCTTGCCCTCTCTCCCAAGGCCACACTAGCGTCCAATACCTGTGCGGTCCCTGTGTgagcgggcaagcccgcgaag                                                                                                              >  1:3733422‑M2/1‑130 (MQ=255)
                                      aGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCGGACTCCCAAGGCCCCACAAGCGTCCAATACCTTGGCGGTCCCTGTGGgagcgggcacgcccgcgacg                                                                                                              <  2:1914275‑M2/150‑21 (MQ=255)
                                      aGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAACCGTCCCATAACTGTGCGGTCCCTGTGGgagcgggcacgcccgcggag                                                                                                              <  2:1914259‑M2/150‑21 (MQ=255)
                                       tgaAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGTTGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaaga                                                                                                             >  1:3374360‑M2/2‑129 (MQ=255)
                                       ggAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaaga                                                                                                             >  2:2604543‑M2/1‑129 (MQ=255)
                                       ggAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGGCCCCTCACCCCAGGCCCCACAAACCTCCAAAAACAGTTCGGTGCCAGTTGGCGCGGGCCAGCCCCCGAaga                                                                                                             >  1:3068690/1‑150 (MQ=255)
                                        gAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaagag                                                                                                            <  2:3852134‑M2/150‑23 (MQ=255)
                                        gAAAAGACAGAAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCGGTGGgagcgggcaagcccgcgaagag                                                                                                            <  2:1368744‑M2/150‑23 (MQ=255)
                                          aaaGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGGTTGGCCCCGCACCCAAGGCCCCCCAAACCGCCAAAAACAGTGCGGTGCCAGGGGGTGCGGGGAAGCCCCCCCAGAGGc                                                                                                          >  1:1204634/1‑150 (MQ=255)
                                                      cTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                              <  2:4207617/150‑1 (MQ=255)
                                                      cTAACGAATTATTCGTCGGGCCGTGCCGGTCGATGGGATGGGGATGTACTCTTCGCGGGCTGGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                              <  2:2634451/150‑1 (MQ=255)
                                                       tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCCCACTCCCCACGCCCCACAAACGTCCAATAACATTTAGGTCCCAGTTGGAGAGGGCAAGCCCGCGCAGAGGCCCCGcccaaccag                                                                                             >  1:3117684/1‑143 (MQ=255)
                                                       tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCCATACCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCgag                                                                                             >  2:1755113/1‑150 (MQ=255)
                                                       tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgag                                                                                             >  1:4133116/1‑150 (MQ=255)
                                                       tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTACCAAGGCCCCACAAGCGTCCAATAACAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCCCGAAGAGGCCCGGCCCTCCgag                                                                                             >  2:533251/1‑150 (MQ=255)
                                                         aCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                           <  2:2841933/150‑1 (MQ=255)
                                                          cGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATg                                                                                          >  2:373268/1‑150 (MQ=255)
                                                              ttattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCCCGCTCCCAAGGCCCCCCAACCGTCCAAAACCAGTGCCGGCCCAGTGGGGGCGGGCCACCCCCCCAAGGGGCCCCGCCCACCCAGATTGccc                                                                                      >  1:3485698/1‑149 (MQ=255)
                                                               tattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAAAACCAGTGCGGTCCCTGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAg                                                                                     >  2:2401230/1‑150 (MQ=255)
                                                               tattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCGCTCTCCCAAGGCCCCCCAAGCGTCCAAAACCAGTGCGGTCCCGGTGGGAGCGGGCAACCCCGCGAAGAAGCCCCGCCATCCCAGATGGCCAg                                                                                     >  2:2401244/1‑150 (MQ=255)
                                                                attCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAgg                                                                                    <  2:717384/150‑1 (MQ=255)
                                                                 ttCGTCGGGCCGTGCCGGTCGGTGGGTTGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACGTGTGCGGTCCCTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                   <  1:4363407/150‑1 (MQ=255)
                                                                 ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                   <  2:4474982/150‑1 (MQ=39)
                                                                 ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                   <  2:1482957/150‑1 (MQ=39)
                                                                 tgcgtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                   <  1:1816400/148‑1 (MQ=255)
                                                                    gtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTCCCAAGGCCCCACAAGCGGCCAATACCAGTGCGGTCCCGGTGGGCGCGGGCCAGCCCGCGAAGAGGCCCGGCCAACCGAGATGGCCAGGCCgg                                                                                >  1:1394820/1‑149 (MQ=255)
                                                                     tcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATACGACATGGCCAGGCCCTc                                                                               >  2:1720048/1‑150 (MQ=255)
                                                                        ggCCGTGCCGGTCGGTGGGAGGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCa                                                                            <  2:2680517/150‑1 (MQ=255)
                                                                         gCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCGGCTCTCCCAAGGCCCCACAAGGGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCat                                                                           <  1:2931164/150‑1 (MQ=255)
                                                                            gTGCCGGTCGGTGGGCTGGGGAGGTACTCTGCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGGCCAAGACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaa                                                                        <  2:3035596/150‑1 (MQ=39)
                                                                                 ggtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                   <  2:3099549/150‑1 (MQ=255)
                                                                                              gggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACACGCTTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTg                                                      <  1:3001512/150‑1 (MQ=255)
                                                                                              gggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTg                                                      <  2:1047245/150‑1 (MQ=255)
                                                                                               ggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGCAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTgg                                                     >  2:2166826/1‑150 (MQ=255)
                                                                                                    tACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCAGGCTACg                                                >  2:2937927/1‑150 (MQ=255)
                                                                                                     aCTCTTCGCGGGCTTGCCCGCTCCCCCAAGGCCCCACAAGCGTCCAATACGTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGc                                               <  1:3659181/150‑1 (MQ=255)
                                                                                                      CTCTTCGCGGGCGTGCCCGCTCTCCCAAGGCCCCACACGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCa                                              >  2:3218815‑M1/1‑150 (MQ=255)
                                                                                                       gctggGCGGGCTGGCCCGCACCCCCACGGCCCCACACGCGCCCATTCCCTGGGCGGTCCGGGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                             >  1:4198306‑M1/6‑150 (MQ=255)
                                                                                                       cCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCGGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                             >  2:882787‑M1/2‑150 (MQ=255)
                                                                                                       TCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGTGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                             >  1:2537211‑M1/1‑150 (MQ=255)
                                                                                                         TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTAAGCAGGt                                           <  2:3283297‑M1/150‑1 (MQ=255)
                                                                                                         TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGAGATGGCCAGGCCGTCCCAGCAATCAGACCAGCACACCCTGGCGACGCAGGt                                           <  2:238320‑M1/150‑1 (MQ=255)
                                                                                                          gcgggggCTGGCCCGCTCCCCCACGGCCCAACACGCGTCCAATACCTGTGGGGTCCCGGGGGGTGGGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                          <  1:117288/146‑1 (MQ=255)
                                                                                                          TCGCGGGCTAGCCCGCTCTCCCAACCCCCCACAAGCGTCCAATACCTCTGCGGTCCCAGTGCGAGCGGGCAAGCCCCCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                          >  2:3549361‑M1/1‑150 (MQ=255)
                                                                                                             CGGGCTGGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTc                                       >  2:2274043‑M1/1‑150 (MQ=255)
                                                                                                                  ttGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaagagtacatccccatcccaccgaccggcacggcccgacgaataattcgttagcctgctgtcttttccttcgcgcaacg                                  <  1:3315477‑M2/150‑97 (MQ=255)
                                                                                                                   tGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaagagtacatccccatcccaccgaccggcacggcccgacgaataattcgttagcctgctgtcttttccttcgcgcaacgg                                 <  2:3774847‑M2/150‑98 (MQ=255)
                                                                                                                   gGCCCGCTCCCCCACGGCCCCACACGCGTCCAACCCGTGTGCGGTCCCTGGGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                 >  1:4205979‑M1/2‑150 (MQ=255)
                                                                                                                   TGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                 >  2:1339601‑M1/1‑150 (MQ=255)
                                                                                                                    GCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                >  1:353186‑M1/1‑150 (MQ=255)
                                                                                                                     cccGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaagagtacatccccatcccaccgaccggcacggcccgacgaataattcgttagcctgctgtcttttccttcgcgcaacggcc                               <  2:580389‑M2/150‑100 (MQ=255)
                                                                                                                      ccGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGgagcgggcaagcccgcgaagagtacatccccatcccaccgaccggcacggcccgacgaataattcgttagcctgctgtcttttccttcgcgcaacggcca                              <  1:1796021‑M2/150‑101 (MQ=255)
                                                                                                                      CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                              <  2:2234473‑M1/150‑1 (MQ=255)
                                                                                                                            tcCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCt                        >  2:1055388/1‑150 (MQ=255)
                                                                                                                              ccACGGCCCCACATGCGTCCAATCCCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGa                      <  1:843556/150‑1 (MQ=255)
                                                                                                                               catggCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                     <  2:428437/147‑1 (MQ=255)
                                                                                                                                aaGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                    >  2:4408274/1‑150 (MQ=255)
                                                                                                                                   gCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcc                 >  1:3161862/1‑150 (MQ=255)
                                                                                                                                      ccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacc              >  2:41577/1‑150 (MQ=255)
                                                                                                                                       cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccaccc             >  2:2091741/1‑149 (MQ=255)
                                                                                                                                       cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacca             >  2:2467933/1‑150 (MQ=255)
                                                                                                                                       cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacca             >  2:849735/1‑150 (MQ=255)
                                                                                                                                        acaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccac            >  1:3925101/1‑150 (MQ=255)
                                                                                                                                          aaGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccacac          >  2:3113901/1‑150 (MQ=255)
                                                                                                                                          aaGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccacac          >  1:2263107/1‑150 (MQ=255)
                                                                                                                                              gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc      <  2:3956835/150‑1 (MQ=255)
                                                                                                                                              gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc      <  1:1073188/150‑1 (MQ=255)
                                                                                                                                              gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc      <  2:3068690/150‑1 (MQ=255)
                                                                                                                                               tCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCg     >  2:3832306/1‑150 (MQ=255)
                                                                                                                                                 cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGc   >  1:3142605/1‑150 (MQ=255)
                                                                                                                                                  aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc  <  1:4408274/150‑1 (MQ=255)
                                                                                                                                                  aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc  >  2:1240924/1‑150 (MQ=255)
                                                                                                                                                  aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc  >  1:1917911/1‑150 (MQ=255)
                                                                                                                                                  aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCAGCGGcc  >  1:473924/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCC  >  NC_002947/3340980‑3341275

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