Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,452 G→C 6.2% intergenic (‑263/+121) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,084,4520GC 6.2% 918.6 / 7.1 300intergenic (‑263/+121)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (140/135);  new base C (20/3);  total (161/139)
Fisher's exact test for biased strand distribution p-value = 7.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGG  >  NC_002947/4084303‑4084601
                                                                                                                                                     |                                                                                                                                                     
  gcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCa                                                                                                                                                     <  1:2631502‑M1/48‑1 (MQ=255)
  gAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCa                                                                                                                                                     <  1:6894165/150‑1 (MQ=255)
   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGCaa                                                                                                                                                    >  2:7787884/1‑150 (MQ=255)
   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCaa                                                                                                                                                    <  2:3753651/150‑1 (MQ=255)
      cccTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCCCCAAGCGGCCCCAGccaatt                                                                                                                                                 >  2:4765199/1‑145 (MQ=255)
        cTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCa                                                                                                                                               <  1:3878406/150‑1 (MQ=255)
          tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                                                             <  1:5782160/150‑1 (MQ=255)
          tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                                                             <  2:1820731/150‑1 (MQ=255)
             ccGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGc                                                                                                                                          <  2:3966438/150‑1 (MQ=255)
              cGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGGGTGGGGGGCGCTAGGGGCCCAAAGCGGCCCCCCCCAGGTCAGCGGcc                                                                                                                                         >  1:3362064/1‑149 (MQ=255)
                        ttGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTAAGCGGCCAAGCTGCaaa                                                                                                                               >  2:3115015/1‑149 (MQ=255)
                           cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCGCAAAATCAGCGGCCAAGCGGCACTCGt                                                                                                                            >  2:6786966/1‑150 (MQ=38)
                           cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGAGAGGGGCGCCAAGCGGCCCCAGCCATTTCAGCCGCACAGCGGCAATCgg                                                                                                                            >  2:7672650/1‑149 (MQ=38)
                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                                           <  1:2715490/150‑1 (MQ=255)
                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                                           <  1:3741201/150‑1 (MQ=255)
                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                                           <  2:7931019/150‑1 (MQ=255)
                             cTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAg                                                                                                                          <  1:5167548/150‑1 (MQ=255)
                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                                         <  2:6005050/150‑1 (MQ=255)
                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                                         <  2:5976040/150‑1 (MQ=255)
                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                                         <  2:1592420/150‑1 (MQ=255)
                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                                         <  1:4404840/150‑1 (MQ=255)
                               tGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAggg                                                                                                                        <  1:5846418/150‑1 (MQ=255)
                               tGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAggg                                                                                                                        <  2:6609797/150‑1 (MQ=255)
                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgggg                                                                                                                       <  2:3207383/150‑1 (MQ=255)
                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgggg                                                                                                                       <  1:8542870/150‑1 (MQ=255)
                                 tCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCACAGCCGCCCCCGCAATTTAAGCGGCAAAGCGGCAAACGTAGGGgc                                                                                                                      >  1:2092957/1‑150 (MQ=255)
                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTAAGCGGCAAAGCTGCAATCGTAGGGgcgc                                                                                                                    >  2:1367917/1‑150 (MQ=255)
                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGGGCGGGGGGAGAGGGCCGCCCAACGGCCCCCCCCAATTAAGAGGCAAAGCGGCAAAAGAAGGGgcgc                                                                                                                    >  1:8097946/1‑150 (MQ=255)
                                      aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTcc                                                                                                                 <  1:5327235/150‑1 (MQ=255)
                                        aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGGGTCGCGGGTGGGCCGCAAAGCGGCCCCCCCAATTTAACCGGCCAAGCTGCCAACGGAGGGGGGcgccgc                                                                                                               >  1:7754055/1‑150 (MQ=255)
                                         aGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCa                                                                                                              <  1:3072310/150‑1 (MQ=255)
                                          gCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAc                                                                                                             <  2:3730059/150‑1 (MQ=255)
                                             gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGGTATGGCCGCAAAACGGCCCCACCAAAATAAGCGGCAAAGCGGCAATAGTAGGGGCGCTCCCCAccca                                                                                                          >  2:7598452/1‑149 (MQ=38)
                                               aGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCttt                                                                                                        <  1:2327951/150‑1 (MQ=255)
                                                ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTt                                                                                                       >  1:7361310‑M2/1‑149 (MQ=255)
                                                 gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtc                                                                                                      >  2:2331978‑M2/1‑148 (MQ=255)
                                                 gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtc                                                                                                      >  1:4267414‑M2/1‑148 (MQ=255)
                                                   ttcagaagtaactcgccaacgggccaactgcggccccgcctgtggggcggcctGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgcg                                                                                                   >  1:1791236‑M1/54‑150 (MQ=255)
                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCCTACCGATAAGGCCGCAAAACCGCCCCCCCAAATTAAGCGGCAAAAATGCCAACGTAGGGGCCCTCCCCACCCATTTagc                                                                                                    >  2:1710323/1‑150 (MQ=255)
                                                      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgac                                                                                                 >  1:4698416‑M2/1‑143 (MQ=255)
                                                      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgac                                                                                                 >  2:2192928‑M2/1‑143 (MQ=255)
                                                      actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGTAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                 >  2:3499564‑M1/51‑150 (MQ=255)
                                                      actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                 >  1:4870614‑M1/51‑150 (MQ=255)
                                                      actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                 >  2:6960068‑M1/51‑150 (MQ=255)
                                                       cTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCGGGAGGGCCGCAAAGCCGCCCCACCAAATTAAAAGGCCAAGTAGCAAACGTACGGGCGCTCCGCACCCATTTAGCGACg                                                                                                >  1:8025645/1‑150 (MQ=255)
                                                          gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                                                             >  1:8718239‑M2/1‑139 (MQ=255)
                                                          gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                                                             >  2:2592327‑M2/1‑139 (MQ=255)
                                                          gccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                             >  2:4932212‑M1/47‑150 (MQ=255)
                                                           tgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaag                                                                                            >  2:1893644‑M2/1‑138 (MQ=255)
                                                           tgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaag                                                                                            >  1:1929731‑M2/1‑138 (MQ=255)
                                                           ccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                                                            >  2:2637516‑M1/46‑150 (MQ=255)
                                                             tCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTAGTGCCGCAAAGCCGCCCCAGCAATTTAAGCGGCAAAGCTGCAAACGTAGGGGCGCTCCGCACCCATTtagcgaccccaggc                                                                                          <  2:1213180‑M2/150‑15 (MQ=255)
                                                              acgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGcc                                                                                         >  2:6120307‑M1/43‑150 (MQ=255)
                                                                 tCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCCCAAAGCGGCCCCCGCAAATTCAGCGGCAAAGCTGCAAACGTAGGGGCGCTCCGCACCCGGTtcgcgacgcaaggccgcg                                                                                      <  2:540799‑M2/150‑19 (MQ=255)
                                                                          tGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAg                                                                             <  2:4249154/150‑1 (MQ=255)
                                                                            tAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTGGCGATATTGCCGCAAATCGGCCCCAGCAATTTAATCGGCAAAGCTGCAATAGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGta                                                                           >  2:703827/1‑149 (MQ=37)
                                                                             cgCGGCCTGGGCGCCTGTGGGGCGGCGTGTGGGTGCGGCCTTGCTTCGCTATCGGGCCGAAACGCGGCCCCAGCAATTTCCGCGCAAAAGCTGCAATCGTTGGGGCGATCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTcc                                                                          <  2:6308789/149‑1 (MQ=37)
                                                                                    tGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAGTTCAGCGGCCAAGCTGCAATCGTAGGGGCGCTCCGCACCCGTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCCACcccg                                                                   >  2:6458907/1‑150 (MQ=25)
                                                                                       gcgcCGGGGGGGCGGCCGGGGGGCGCGGCCTTGGGTCGCGTTCGGGCCGCAACGCGGCCCCAGCATTTCCCGCGGCAATGCTGAAATGGTCGGGGCGCCCCGCACCGTTTTCGCGCCGCAAGCCCGCTCCCACAAAGTCCAACCGCGccc                                                                <  2:1487256/150‑1 (MQ=255)
                                                                                             gtgGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACagg                                                          <  2:5691102/150‑1 (MQ=255)
                                                                                                  gcggcCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAgggg                                                     <  2:1374692/150‑1 (MQ=255)
                                                                                                  gcggcCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAgggg                                                     <  2:6538299/150‑1 (MQ=255)
                                                                                                  gcggcCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCCCAAAGCGGCCCCACCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAgggg                                                     >  2:5098862/1‑150 (MQ=255)
                                                                                                  gcggcCTGTGGGAGCGGCCGTGCGGCCCCGTAGGGGCCCCAACCCGCCCCCCCAAATTCAGCGGCCAACCTGCCATCGGGGGGGCGCGACGCACCCCGGTCGCGCCGCCAGGCCGCGCCCACCCACGGCAACCCCGCCCCCCACGAgggg                                                     >  2:4521972/1‑150 (MQ=255)
                                                                                                          tgGGAGCGGGCTTGCGTCGCGATACGGCCCCAAAACGGGCCCCCCAATTTCAGCGGCAAAGCTTCCATGGTAGGGGCGCGCCCCACCCTATTAGCGGCGCCAGGCCGCGCCCACAAAGAGCAACCGCGCCCCACCGGAGGGGCCGCAGtt                                             >  2:7496880/1‑150 (MQ=255)
                                                                                                            ggAGCGGCCTTTCGTCGCGATATTGCCGCAAATCGGCCCCAGCAATTTCAGCGGCAAATCTGCAATCGTAGGGGCGCTCCGCACCCTTTTAGCGAAGCAAGGCCGCTCCCACAAATTCCAACCGCGCCCCACAGGAGGGGCCGCAGGTgg                                           >  2:668838/1‑150 (MQ=37)
                                                                                                              aGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCACCAAATTCAGCGGCAAACCAGCAATACAAGGGGCGCTCCGCACCCATTTCGCGACCCAAGGCCCCTCCCACAAAGTCCAACCCCCCCCCACACGAGGGGCCCCACTTGGcc                                         >  1:3098122/1‑150 (MQ=21)
                                                                                                              AGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGcc                                         <  2:2024545‑M1/150‑1 (MQ=255)
                                                                                                               GCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGccc                                        <  2:5850670‑M1/150‑1 (MQ=255)
                                                                                                                  GCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGAAATCTTAGGGGCGCTCCGCACCATTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGtt                                     >  2:7311343‑M1/1‑150 (MQ=255)
                                                                                                                   ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgac                                    <  1:1679231‑M2/150‑69 (MQ=255)
                                                                                                                   CCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTg                                    <  1:6537947‑M1/150‑1 (MQ=255)
                                                                                                                     ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacac                                  <  1:1059153‑M2/150‑71 (MQ=255)
                                                                                                                     gtGCGTCGCGAGTGGGCCGCAACGCGGCCCCAGCATGTGCCGCGGCCCCGCTGCAATCGGCGGGGCGCGCCGCACCCTTTTCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                  <  2:8372098/149‑1 (MQ=255)
                                                                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                  <  2:1297657‑M1/150‑1 (MQ=255)
                                                                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                  <  2:1160216‑M1/150‑1 (MQ=255)
                                                                                                                         GTCGCGGTAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGt                              <  2:6502890‑M1/150‑1 (MQ=255)
                                                                                                                         GTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGt                              <  2:8693222‑M1/150‑1 (MQ=255)
                                                                                                                           CGCGATATGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTa                            >  2:2815620‑M1/1‑150 (MQ=255)
                                                                                                                            GCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTAc                           <  1:5090796‑M1/150‑1 (MQ=255)
                                                                                                                              GATCGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGAAAAGCTGAAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACtt                         >  2:955320‑M1/1‑150 (MQ=255)
                                                                                                                              GATCGGGCCGCAAAGCGGCCCCAGCAATTACAGCGGCAACGCTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACtt                         >  2:846593‑M1/1‑150 (MQ=255)
                                                                                                                               attgggCCGCAACGCGGCCCCAGCAATTTCCCCGGCAAAGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc                        <  2:8430078/147‑1 (MQ=39)
                                                                                                                               ATAGGGCCGCAACGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTAGGGGCGCCCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc                        >  2:3085016‑M1/1‑150 (MQ=255)
                                                                                                                                  gggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGa                     >  1:7265022/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                    >  2:3672124/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                    >  2:8142554/1‑150 (MQ=255)
                                                                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                    >  1:7041771/1‑150 (MQ=255)
                                                                                                                                   ggCCGAAAAGCGGCCCCAGCACTTTCCGCGGCAAAGCTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGCCGCATGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                    <  2:8228849/150‑1 (MQ=255)
                                                                                                                                    gCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                   >  2:3170159/1‑150 (MQ=255)
                                                                                                                                     ccGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg                  >  2:2492333/1‑150 (MQ=255)
                                                                                                                                       gCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGc                >  2:7108874/1‑150 (MQ=255)
                                                                                                                                        caccgCGGCCCCCTCATTTGCAGCGGAAACGCTGCATGCGGCGGGGCGCCCCGCCCCCTTTCCGCCCCCCACGCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCa               <  1:6245988/146‑1 (MQ=255)
                                                                                                                                        cAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCa               <  2:8590915/150‑1 (MQ=255)
                                                                                                                                          aaGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGc             >  2:1090625/1‑150 (MQ=255)
                                                                                                                                               gCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACg        >  1:8199907/1‑150 (MQ=255)
                                                                                                                                                ccctAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGc       <  2:8655019/150‑1 (MQ=255)
                                                                                                                                                ccccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGc       >  1:5640539/1‑150 (MQ=255)
                                                                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt     >  2:3607297/1‑150 (MQ=255)
                                                                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc    >  1:8622119/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGG  >  NC_002947/4084303‑4084601

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: