Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,353 C→G 8.4% intergenic (‑91/+28) 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,3530CG 8.4% 305.6 / 15.6 95intergenic (‑91/+28)sbmC/dacDDNA gyrase inhibitor/D‑alanyl‑D‑alanine carboxypeptidase DacD
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base G (2/6);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATT  >  NC_000913/2081204‑2081498
                                                                                                                                                     |                                                                                                                                                 
ccAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCAc                                                                                                                                                   <  2:1504501/150‑1 (MQ=255)
ccAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCAc                                                                                                                                                   >  1:267929/1‑150 (MQ=255)
    cGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCt                                                                                                                                               >  2:714879/1‑150 (MQ=255)
       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:402953/150‑3 (MQ=255)
                 tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                  <  1:440094/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                                                                                                              <  1:538178/150‑1 (MQ=255)
                                     ccTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTa                                                                                                              <  2:1388260/150‑1 (MQ=255)
                                            aaTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                       <  1:1425987/150‑1 (MQ=255)
                                             aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                      >  1:1419553/1‑150 (MQ=255)
                                             aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                      >  2:1592728/1‑150 (MQ=255)
                                             aTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                      >  1:456290/1‑150 (MQ=255)
                                             aTCTCGTAGTTCATGACGCCTCCATTGATAATGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaa                                                                                                      >  2:1055897/1‑150 (MQ=255)
                                              tctcGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaaa                                                                                                     <  2:245682/150‑1 (MQ=255)
                                                 cGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataa                                                                                                  <  1:391750/150‑1 (MQ=255)
                                                 cGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataa                                                                                                  >  1:1152772/1‑150 (MQ=255)
                                                  gTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataat                                                                                                 <  1:273283/150‑1 (MQ=255)
                                                   tAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATc                                                                                                <  2:1057646/150‑1 (MQ=255)
                                                     gTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCAc                                                                                              >  2:1300804/1‑150 (MQ=255)
                                                        cATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCa                                                                                           <  2:1271625/150‑1 (MQ=255)
                                                        cATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCa                                                                                           <  1:1225923/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                                                                                                                                    >  2:226350/1‑101 (MQ=38)
                                                                cTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgct                                                                                                                                    <  1:226350/101‑1 (MQ=38)
                                                                  ccATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGaaaa                                                                                 >  1:24977/1‑150 (MQ=255)
                                                                  ccATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGaaaa                                                                                 >  2:432867/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:760653/150‑1 (MQ=255)
                                                                   cATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                <  2:422389/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                                                  <  1:729801/150‑1 (MQ=255)
                                                                                                 tGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTa                                                  >  2:849251/1‑150 (MQ=255)
                                                                                                   cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAAc                                                >  2:1097091/1‑150 (MQ=255)
                                                                                                   cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAAc                                                >  1:1256330/1‑150 (MQ=255)
                                                                                                      aaTTAACACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGGCTGCTGCAAAAGAAAGGGCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             <  1:500342/150‑1 (MQ=255)
                                                                                                      aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             >  2:1049389/1‑150 (MQ=255)
                                                                                                      aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             <  1:494820/150‑1 (MQ=255)
                                                                                                      aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             >  1:253873/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:391380/1‑150 (MQ=255)
                                                                                                                   aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                >  1:1273855/1‑150 (MQ=255)
                                                                                                                   aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                >  2:13030/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)
                                                                                                                                                     |                                                                                                                                                 
CCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATT  >  NC_000913/2081204‑2081498

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: