Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,355 C→A 8.2% intergenic (‑93/+26) 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,3550CA 8.2% 325.5 / 13.3 98intergenic (‑93/+26)sbmC/dacDDNA gyrase inhibitor/D‑alanyl‑D‑alanine carboxypeptidase DacD
Reads supporting (aligned to +/- strand):  ref base C (45/45);  new base A (2/6);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCC  >  NC_000913/2081208‑2081503
                                                                                                                                                   |                                                                                                                                                    
cGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCt                                                                                                                                                    >  2:714879/1‑150 (MQ=255)
   gaTGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGcc                                                                                                                                                 <  1:1279702/150‑1 (MQ=255)
   gaTGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGcc                                                                                                                                                 <  2:1457180/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: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                                                                                                <  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                                                                                         >  2:1583476/1‑150 (MQ=255)
                                                           ccTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgaga                                                                                         >  1:736161/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: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                                                       <  1:729801/150‑1 (MQ=255)
                                                                                             tGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTa                                                       >  2:849251/1‑150 (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                                                  >  2:1049389/1‑150 (MQ=255)
                                                                                                  aaTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAg                                                  <  1:494820/150‑1 (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                                         >  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)
                                                                                                                                                 cTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGTCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCAc                                   <  1:1604846/118‑1 (MQ=255)
                                                                                                                                                 cTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGTCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCAc                                   >  2:1604846/1‑118 (MQ=255)
                                                                                                                                                  tCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTcccc  <  1:1320160/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCC  >  NC_000913/2081208‑2081503

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: