Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,349 C→G 8.2% intergenic (‑87/+32) 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,3490CG 8.2% 316.9 / 15.7 97intergenic (‑87/+32)sbmC/dacDDNA gyrase inhibitor/D‑alanyl‑D‑alanine carboxypeptidase DacD
Reads supporting (aligned to +/- strand):  ref base C (44/45);  new base G (2/6);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

GGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATT  >  NC_000913/2081201‑2081498
                                                                                                                                                    |                                                                                                                                                     
gggCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCg                                                                                                                                                      <  1:574812/150‑1 (MQ=255)
   ccAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCAc                                                                                                                                                   >  1:267929/1‑150 (MQ=255)
   ccAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCAc                                                                                                                                                   <  2:1504501/150‑1 (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: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: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                                                                                           <  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                                                                                    >  2:1583476/1‑150 (MQ=255)
                                                                  ccTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgaga                                                                                    >  1:736161/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:574827/6‑150 (MQ=255)
                                                                      caactaTACTGTATATAAAAACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAAc                                                                                >  1:1342914/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:253873/1‑150 (MQ=255)
                                                                                                         aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             <  1:494820/150‑1 (MQ=255)
                                                                                                         aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             >  2:1190124/1‑150 (MQ=255)
                                                                                                         aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                             >  2:1049389/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                                    >  1:1542775/1‑150 (MQ=255)
                                                                                                                  cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                    >  2:477151/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)
                                                                                                                                                    |                                                                                                                                                     
GGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATT  >  NC_000913/2081201‑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: ‑