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

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

GATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTAC  >  NC_000913/2081211‑2081506
                                                                                                                                                     |                                                                                                                                                  
gaTGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGcc                                                                                                                                                    <  1:772956/150‑1 (MQ=255)
  ttgAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAg                                                                                                                                                  >  2:897570/3‑150 (MQ=255)
   ggAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGt                                                                                                                                                 <  2:1268316/150‑1 (MQ=255)
   ggAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGt                                                                                                                                                 <  2:320992/150‑1 (MQ=255)
    gAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTc                                                                                                                                                <  2:381575/150‑1 (MQ=255)
       aCCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCTTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGTATGAGTGTTGATTATGCGAGCAGACTCGCACTCGTGCCAGTctgc                                                                                                                                             <  2:1092659/150‑1 (MQ=255)
        ccTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTctgct                                                                                                                                            >  1:168147/1‑150 (MQ=255)
          tGCACCGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          <  2:718916/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          <  2:1277645/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          <  2:1817949/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          <  2:749850/150‑3 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCtgctgc                                                                                                                                          >  2:424961/1‑150 (MQ=255)
          tGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCtgctgc                                                                                                                                          <  2:1875242/150‑1 (MQ=255)
           gCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCa                                                                                                                                         <  1:1158861/150‑1 (MQ=255)
                ggTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAaaaga                                                                                                                                    <  1:1792325/150‑1 (MQ=255)
                             tCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCtt                                                                                                                       >  1:590631/1‑150 (MQ=255)
                               cTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAaaag                                                                                                                                     >  2:73979/1‑134 (MQ=255)
                               cTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAaaag                                                                                                                                     <  1:73979/134‑1 (MQ=255)
                                   ttAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTAtggtg                                                                                                                 <  2:550847/150‑1 (MQ=255)
                                    tAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTAtggtgg                                                                                                                >  1:1326921/1‑150 (MQ=255)
                                     aaTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                               <  1:1616420/150‑1 (MQ=255)
                                     aaTCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                               <  1:138380/150‑1 (MQ=255)
                                       tctcGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGaaa                                                                                                             <  1:891201/150‑1 (MQ=255)
                                         tcGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCAGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaata                                                                                                           <  1:1826240/150‑1 (MQ=255)
                                         tcGTAGTTCATGACGCCTCCATTGATACTGTTTTAATATACAGTATTGTTGCAAATTTAAAACCACAAGGAATGAGTGTTGTTTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGTCCTTATGGGGGAaata                                                                                                           <  2:1769328/150‑1 (MQ=255)
                                           gTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataat                                                                                                         <  1:1092107/150‑1 (MQ=255)
                                           gTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAaataat                                                                                                         <  1:333036/150‑1 (MQ=255)
                                            tAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATc                                                                                                        <  1:962030/150‑1 (MQ=255)
                                              tttCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCAc                                                                                                      <  1:859548/149‑1 (MQ=255)
                                              gTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCAc                                                                                                      >  2:570804/1‑150 (MQ=255)
                                                tCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                               <  1:389960/139‑1 (MQ=255)
                                                tCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGa                                                                                                               >  2:389960/1‑139 (MQ=255)
                                                tCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCAGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTc                                                                                                    >  1:1097383/1‑150 (MQ=255)
                                                 cATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCa                                                                                                   >  1:757746/1‑150 (MQ=255)
                                                  aTGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAg                                                                                                  >  2:525487/1‑150 (MQ=255)
                                                          tCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgag                                                                                             <  2:863801/147‑1 (MQ=255)
                                                          tCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCgag                                                                                             >  1:863801/1‑147 (MQ=255)
                                                           ccATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAATCGACAAGGAATGAGTGTTGATTATGAGAGCAGAGTCGCACTCCTGGCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATGAGTCAGGAGAGaaaa                                                                                         <  2:1480850/150‑1 (MQ=255)
                                                                  tACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACAtgctg                                                                                  >  1:1311266/1‑150 (MQ=255)
                                                                  tACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACAtgctg                                                                                  <  2:1581212/150‑1 (MQ=255)
                                                                   aCTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACAtgctgc                                                                                 <  2:1075269/150‑1 (MQ=255)
                                                                    cTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGcc                                                                                >  2:1424315/1‑150 (MQ=255)
                                                                    cTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGAGATAAAACATGCTGcc                                                                                <  2:1555630/150‑1 (MQ=255)
                                                                      gTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCtt                                                                              >  1:1979747/1‑150 (MQ=255)
                                                                         tttATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTccc                                                                           >  2:1023455/1‑150 (MQ=255)
                                                                          ttATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTcccc                                                                          <  1:1424315/150‑1 (MQ=255)
                                                                              ataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGGCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACa                                                                      <  2:140577/150‑1 (MQ=255)
                                                                              ataCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACa                                                                      <  1:1575250/150‑1 (MQ=255)
                                                                               taCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAg                                                                     >  1:1645361/1‑150 (MQ=255)
                                                                                aCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCg                                                                         <  2:402314/145‑1 (MQ=255)
                                                                                aCAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCg                                                                         >  1:402314/1‑145 (MQ=255)
                                                                                 cAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          <  2:1940265/79‑3 (MQ=25)
                                                                                 cAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTGGCAGGAGTGCGAGTctgctcg                                                                                                                                          >  1:1940265/1‑77 (MQ=25)
                                                                                 cAGTATAGTTGCAAAATCAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGAt                                                                   >  1:149258/1‑150 (MQ=255)
                                                                                      taGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAg                                                              <  2:1668279/150‑1 (MQ=255)
                                                                                      taGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAg                                                              >  1:537602/1‑150 (MQ=255)
                                                                                      taGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAg                                                              <  1:221432/150‑1 (MQ=255)
                                                                                           gCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTaa                                                         >  2:92071/1‑150 (MQ=255)
                                                                                           gCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTaa                                                         >  1:1288831/1‑150 (MQ=255)
                                                                                            cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAAc                                                        <  1:92071/150‑1 (MQ=255)
                                                                                            cAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTACTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAAACACTCACGCGAGAAAACATGCTGCCTTCCCAGACAGAATCCAGGGTAAc                                                        >  1:1833285/1‑150 (MQ=255)
                                                                                                aTTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGc                                                    <  2:168147/150‑1 (MQ=255)
                                                                                                  tAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCgg                                                  <  2:786269/150‑1 (MQ=255)
                                                                                                  tAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCgg                                                  <  2:1047922/150‑1 (MQ=255)
                                                                                                        cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                            >  2:1922038/1‑150 (MQ=255)
                                                                                                        cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                            >  2:1982539/1‑150 (MQ=255)
                                                                                                        cacaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTg                                            >  2:15202/1‑150 (MQ=255)
                                                                                                         acaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTgc                                           >  1:1771375/1‑150 (MQ=255)
                                                                                                         acaAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTgc                                           <  2:100792/150‑1 (MQ=255)
                                                                                                          caAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTgcg                                          >  2:1778652/1‑150 (MQ=255)
                                                                                                          caAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGGCTGCTGCAAAAGAAAGGGCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTgcg                                          <  2:1174197/150‑1 (MQ=255)
                                                                                                            aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                        >  1:1496148/1‑150 (MQ=255)
                                                                                                            aGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGcc                                        >  1:1066771/1‑150 (MQ=255)
                                                                                                               taTGAGTGTTGTTTATGAGAGGAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCAGACAGATTCCAGGGTAACCAGCTGCCAGTTCGCCAcc                                     <  2:537602/149‑1 (MQ=255)
                                                                                                               aaTGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCAcc                                     <  2:1311266/150‑1 (MQ=255)
                                                                                                                   aGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGtt                                 <  1:1778652/150‑1 (MQ=255)
                                                                                                                   aGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGtt                                 <  1:1600960/150‑1 (MQ=255)
                                                                                                                     tgtTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTa                               <  1:897570/150‑1 (MQ=255)
                                                                                                                     tgtTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTa                               <  1:563091/150‑1 (MQ=255)
                                                                                                                        tGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCa                            >  1:1525812/1‑150 (MQ=255)
                                                                                                                         gATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCAc                           >  1:962270/1‑150 (MQ=255)
                                                                                                                              tGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCg                      <  1:1260723/150‑1 (MQ=255)
                                                                                                                               gCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGt                     <  1:1880667/150‑1 (MQ=255)
                                                                                                                                gcaGCAGACTGGCAGGAGTGCGAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTa                    >  2:143175/3‑150 (MQ=255)
                                                                                                                                cGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTa                    >  1:65427/1‑150 (MQ=255)
                                                                                                                                cGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTa                    >  1:1299676/1‑150 (MQ=255)
                                                                                                                                 gAGCAGACTCGCACTGCTGCCAGTCTGCTGCAAAAGAAAGGGCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTaa                   <  2:808936/150‑1 (MQ=255)
                                                                                                                                   gCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAg                 <  2:149258/150‑1 (MQ=255)
                                                                                                                                   gCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAg                 <  1:1347886/150‑1 (MQ=255)
                                                                                                                                   gCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCACGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAg                 >  2:1359024/1‑150 (MQ=255)
                                                                                                                                       aCTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCa             <  1:570804/150‑1 (MQ=255)
                                                                                                                                          cGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAAtt          >  2:1573243/1‑150 (MQ=255)
                                                                                                                                          cGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCATTATGGTGGAAATAATCACTCAGGCGAGAAAACATGATGCCTTCCACGACAGATTCCAGGGTAACCAGCGGCAAGTGCGCCACATGTTTATCACGGTAGTAAAGTTCAAtt          >  2:1565923/1‑150 (MQ=255)
                                                                                                                                             aCTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCGGAGAGATTCCAGGGTAACCAGCGGCCAGTTCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTcc       <  2:1097383/150‑1 (MQ=255)
                                                                                                                                               tCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTcccc     <  2:1525812/150‑1 (MQ=255)
                                                                                                                                                  tGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTAc  >  1:523163/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTTTTATATACAGTATAGTTGCAAAATTAAAACCACAAGGAATGAGTGTTGATTATGCGAGCAGACTCGCACTCCTGCCAGTCTGCTGCAAAAGAAAGGTCAGGCCTTATGGTGGAAATAATCACTCAGGCGAGAAAACATGCTGCCTTCCCCGACAGATTCCAGGGTAACCAGCGGCCAGTGCGCCACCTGTTTATCACGGTCGTAAAGTTCAATTTCCCCTAC  >  NC_000913/2081211‑2081506

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: