Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,360 C→G 8.3% intergenic (‑98/+21) sbmC ← / ← dacD DNA gyrase inhibitor/D‑alanyl‑D‑alanine carboxypeptidase DacD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,081,3600CG 8.3% 304.0 / 15.6 96intergenic (‑98/+21)sbmC/dacDDNA gyrase inhibitor/D‑alanyl‑D‑alanine carboxypeptidase DacD
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base G (2/6);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTACCC  >  NC_000913/2081211‑2081508
                                                                                                                                                     |                                                                                                                                                    
gaTGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGcc                                                                                                                                                      <  2:1457180/150‑1 (MQ=255)
gaTGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGcc                                                                                                                                                      <  1:1279702/150‑1 (MQ=255)
  tGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAg                                                                                                                                                    <  2:1548173/150‑1 (MQ=255)
    gAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTc                                                                                                                                                  <  2:515891/150‑1 (MQ=255)
      aaCCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTctg                                                                                                                                                >  1:567330/1‑150 (MQ=255)
        ccTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgct                                                                                                                                              <  1:864981/150‑1 (MQ=255)
        ccTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgct                                                                                                                                              <  1:1568261/150‑1 (MQ=255)
        ccTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGAAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTctgct                                                                                                                                              <  2:808330/150‑1 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                            <  1:440094/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                            <  1:402953/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCGTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTTTGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                            <  1:1127138/150‑3 (MQ=255)
            cAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCaa                                                                                                                                          >  1:704553/1‑150 (MQ=255)
             aaCGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCaaa                                                                                                                                         >  1:439039/1‑150 (MQ=255)
                   aCGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAaaagaaag                                                                                                                                   >  2:1603205/1‑150 (MQ=255)
                       ttCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCa                                                                                                                               <  1:1228670/150‑1 (MQ=255)
                        tctcTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAg                                                                                                                              <  1:315109/150‑1 (MQ=255)
                         ctctTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAgg                                                                                                                             >  1:1388260/1‑150 (MQ=255)
                           ctTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGcc                                                                                                                           >  2:497341/1‑150 (MQ=255)
                           ctTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGcc                                                                                                                           >  1:423116/1‑150 (MQ=255)
                              ccTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTa                                                                                                                        <  2:1388260/150‑1 (MQ=255)
                              ccTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTa                                                                                                                        <  1:538178/150‑1 (MQ=255)
                                     aaTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                                 <  1:1425987/150‑1 (MQ=255)
                                      aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                                >  1:456290/1‑150 (MQ=255)
                                      aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                                >  1:1419553/1‑150 (MQ=255)
                                      aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                                >  2:1592728/1‑150 (MQ=255)
                                      aTCTCGTAGTTCATGACGCCTCCATTGATAATGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                                >  2:1055897/1‑150 (MQ=255)
                                       tctcGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaaa                                                                                                               <  2:245682/150‑1 (MQ=255)
                                          cGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataa                                                                                                            >  1:1152772/1‑150 (MQ=255)
                                          cGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataa                                                                                                            <  1:391750/150‑1 (MQ=255)
                                           gTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataat                                                                                                           <  1:273283/150‑1 (MQ=255)
                                            tAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATc                                                                                                          <  2:1057646/150‑1 (MQ=255)
                                              gTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCAc                                                                                                        >  2:1300804/1‑150 (MQ=255)
                                                 cATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCa                                                                                                     <  1:1225923/150‑1 (MQ=255)
                                                 cATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCa                                                                                                     <  2:1271625/150‑1 (MQ=255)
                                                  aTGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAg                                                                                                    <  2:423116/150‑1 (MQ=255)
                                                     aCGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCg                                                                                                 <  2:446414/150‑1 (MQ=255)
                                                      cGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCga                                                                                                >  2:1242518/1‑150 (MQ=255)
                                                       gCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgag                                                                                               <  1:1130798/150‑1 (MQ=255)
                                                        ccTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgaga                                                                                              >  1:736161/1‑150 (MQ=255)
                                                        ccTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgaga                                                                                              >  2:1583476/1‑150 (MQ=255)
                                                         cTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgct                                                                                                                                              <  1:226350/101‑1 (MQ=38)
                                                         cTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgct                                                                                                                                              >  2:226350/1‑101 (MQ=38)
                                                           ccATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGaaaa                                                                                           >  2:432867/1‑150 (MQ=255)
                                                           ccATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGaaaa                                                                                           >  1:24977/1‑150 (MQ=255)
                                                            caactaTACTGTATATAAAAACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                          >  2:876739/6‑150 (MQ=255)
                                                            caactaTACTGTATATAAAAACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                          >  1:1342914/6‑150 (MQ=255)
                                                            caactaTACTGTATATAAAAACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                          >  1:574827/6‑150 (MQ=255)
                                                            cATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                          <  2:422389/150‑1 (MQ=255)
                                                            cATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                          <  2:760653/150‑1 (MQ=255)
                                                                   aCTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACAtgctgc                                                                                   <  1:861038/150‑1 (MQ=255)
                                                                         tttATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTccc                                                                             <  2:456290/150‑1 (MQ=255)
                                                                          ttATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTcccc                                                                            <  1:497341/150‑1 (MQ=255)
                                                                           tatataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCg                                                                           >  2:761005/1‑150 (MQ=255)
                                                                           tatataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACGCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCg                                                                           >  1:1432516/1‑150 (MQ=255)
                                                                            atataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCTGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGa                                                                          <  2:870112/150‑1 (MQ=255)
                                                                            atataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGa                                                                          <  2:1197264/150‑1 (MQ=255)
                                                                             tataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGAc                                                                         <  2:439039/150‑1 (MQ=255)
                                                                              ataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACa                                                                        >  2:19006/1‑150 (MQ=255)
                                                                                 cAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGAt                                                                     >  1:364652/1‑150 (MQ=255)
                                                                                  aGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGAtt                                                                    >  2:351079/1‑150 (MQ=255)
                                                                                      taGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAg                                                                <  2:617778/150‑1 (MQ=255)
                                                                                        gTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAggg                                                              <  2:752032/150‑1 (MQ=255)
                                                                                          tGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTa                                                            >  2:849251/1‑150 (MQ=255)
                                                                                          tGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTa                                                            <  1:729801/150‑1 (MQ=255)
                                                                                            cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAAc                                                          >  1:1256330/1‑150 (MQ=255)
                                                                                            cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAAc                                                          >  2:1097091/1‑150 (MQ=255)
                                                                                               aaTTAACACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGGCTGCTGCAAAAGAAAGGGCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                       <  1:500342/150‑1 (MQ=255)
                                                                                               aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                       <  1:494820/150‑1 (MQ=255)
                                                                                               aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                       >  1:253873/1‑150 (MQ=255)
                                                                                               aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                       >  2:1049389/1‑150 (MQ=255)
                                                                                               aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                       >  2:1190124/1‑150 (MQ=255)
                                                                                                  tAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCgg                                                    <  2:253873/150‑1 (MQ=255)
                                                                                                    aaaCCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGcc                                                  <  1:876739/150‑1 (MQ=255)
                                                                                                     aaCCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCa                                                 >  2:124087/1‑150 (MQ=255)
                                                                                                      aCCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTTCAGGGTAACCAGCGGCCAg                                                <  1:714879/150‑1 (MQ=255)
                                                                                                        cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                              >  2:477151/1‑150 (MQ=255)
                                                                                                        cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                              >  1:1542775/1‑150 (MQ=255)
                                                                                                         acaAGGAATGAGTTTTGATTATGCGAGAAGAATCCCACTCCTGCCAGTCTGCAGCAAAACAAAGGTCAGCCAATATGGAGGAAATAATAAATCAGGCTAGAAAACATGCTGCCTTCCCCGACAGATTCAAGGTTAACAAGCGCCCAGTgc                                             >  2:1603246/1‑150 (MQ=255)
                                                                                                            aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                          >  2:13030/1‑150 (MQ=255)
                                                                                                            aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                          >  2:391380/1‑150 (MQ=255)
                                                                                                            aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                          >  1:1273855/1‑150 (MQ=255)
                                                                                                                aTGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGCCAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCt                                      >  1:1434647/1‑150 (MQ=255)
                                                                                                                   aGTGTTGATTCTGCGAGCAGACTCGCACTCCTGCCAGTCTGCTTCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGtt                                   <  2:1632566/150‑1 (MQ=255)
                                                                                                                     tgtTGATTATGCGAGCAGACTCGCACTCCTGCCAGGCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTa                                 <  2:1464936/150‑1 (MQ=255)
                                                                                                                        tGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCa                              >  2:1101169/1‑150 (MQ=255)
                                                                                                                          aTTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACg                            <  2:378484/150‑1 (MQ=255)
                                                                                                                           ttATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACgg                           <  1:849251/150‑1 (MQ=255)
                                                                                                                            tATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGt                          <  2:515607/150‑1 (MQ=255)
                                                                                                                             aTGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTc                         <  2:704553/150‑1 (MQ=255)
                                                                                                                                gcaGCAGACTGGCAGGAGTGCGAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTa                      >  1:1550950/3‑150 (MQ=255)
                                                                                                                                    cAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGt                  >  1:75707/1‑150 (MQ=255)
                                                                                                                                      gACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTc                <  2:369021/150‑1 (MQ=255)
                                                                                                                                         tCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAAt             <  2:1273855/150‑1 (MQ=255)
                                                                                                                                          cGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCCGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGGAAAGTTCAAtt            >  2:1302725/1‑150 (MQ=255)
                                                                                                                                              cTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGTCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCAc                                        <  1:1604846/118‑1 (MQ=255)
                                                                                                                                              cTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGTCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCAc                                        >  2:1604846/1‑118 (MQ=255)
                                                                                                                                               tCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTcccc       <  1:1320160/150‑1 (MQ=255)
                                                                                                                                                   gCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTAcc   <  1:1055897/150‑1 (MQ=255)
                                                                                                                                                    ccAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTAccc  <  2:75707/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTACCC  >  NC_000913/2081211‑2081508

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