Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,825,328 G→A S83S (TCC→TCT)  mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,825,3280GA88.7% 122.7 / 4.5 56S83S (TCC→TCT) mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (20/27);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGAA  >  NC_000913/2825082‑2825573
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
ccAGCGGGGTTCAACGCCGATAATCCCGACGATAATGTCCGGCGGTACTCCATCCACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTGGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGCCCTGATGTCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTgagaga                                                                                                                                                                                                                                                   <  1:283278/251‑1 (MQ=255)
         tccaacccctaTAACCCCGACGATAATTACCGGCGGTACTCCATACCCCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTGCGGTCCTTATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGGGATAGAATTTcc                                                                                                                                                                                                                                           <  2:218032/240‑1 (MQ=255)
          gcAACGCCGATAATCCCGAAGATAATTTCCGGTGGTACGCCAGACACCTGACACGCGCGATTCAACGCATTTTCATACTGATTACAGAAAACCACACCGTTCTGCATGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGTCGTTAGGACCTGATGGCGGTTTCACCGATGTTGTTGGTGCCTGGTTATCCATCAGCCGCAGTCCCGAATCC‑AGACGGTTCGCCTGAGAGAGAACTTCCt                                                                                                                                                                                                                                          <  2:156433/249‑1 (MQ=255)
                ccGATACTCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACTCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCTCGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAActg                                                                                                                                                                                                                                   <  2:143272/251‑1 (MQ=255)
                ccGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGCGAGAGAATTTCCTGCAActg                                                                                                                                                                                                                                   >  1:87136/1‑251 (MQ=255)
                  gATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAATCCACACCGTTCTGCACGTTTTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAActgct                                                                                                                                                                                                                                 <  1:231227/251‑1 (MQ=255)
                       caccGAAGATAATTCCCCGGGGTACTCCATCCACCTGCCCCGCGCGATTCAACTCATCTTCATGCTGATTCCAGAAATCCGCGCCGTTCTGCACGGGGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGTCCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGTTTATCCATCAGCCGCAGTACCGAATGC‑AGGCTCTTCTCCTGAGAGAGAATGTCCTGCAACTGCTGACGa                                                                                                                                                                                                                            <  2:427848/249‑1 (MQ=255)
                               aTAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAAcc                                                                                                                                                                                                                    <  1:128901/251‑1 (MQ=255)
                                tAATTTCCGGCGGTACTCCATACCCCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGGGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑CGACGCTTCGCCTGAGGGGGACTTTCCTGCAACTGCGGACGATCGAAAc                                                                                                                                                                                                                     >  2:342237/1‑249 (MQ=255)
                                tAATTTCCGGCGGTACTCCATAAACCTGCCACGCGCGATTCCACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAACTTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGGTTTCCATCCGCAGCCGGACCGCATCCCAG‑CGCTTCGCCTGAGAGAGAATTTCCTCCCACTACTGACgcgagcacccc                                                                                                                                                                                                                   >  2:190406/1‑241 (MQ=255)
                                 aaTTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCgg                                                                                                                                                                                                                  >  1:113982/1‑250 (MQ=255)
                                      ccGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGtgttt                                                                                                                                                                                                              <  1:52639/250‑1 (MQ=255)
                                         gcggTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTc                                                                                                                                                                                                          <  1:110808/251‑1 (MQ=255)
                                          cggTACTCCGTACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCAATCAGCCGCAGTACCGAATCC‑AAACGCTTCGCCTGCGTGAGAATTTCCTTCAACTGCTGACGATCGAAAACGTGTTTGTTc                                                                                                                                                                                                          >  1:342121/1‑250 (MQ=255)
                                                 ccATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGAAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGAAACTGCTGAAGATCGAAACCGAGTTTGTTCACCAtttt                                                                                                                                                                                                  >  1:223846/1‑251 (MQ=255)
                                                  cATACCCCTGTCACGCGCGATTCAACGCATCTTCATACGGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTGTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTg                                                                                                                                                                                                 <  1:22653/251‑1 (MQ=255)
                                                  cATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGCTGGGGATGGAGCCTGGTTATCCATCAACCGCAGTACCGAATCC‑AGACGCTTCGCCTGCGAGAGAATTTCCTACACCTTCTGACGAACCAAACCGGGTTTGTTCACCCttttt                                                                                                                                                                                                 >  2:328583/1‑250 (MQ=255)
                                                            ccACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAAct                                                                                                                                                                                       >  1:399023/1‑251 (MQ=255)
                                                               cgcgcgATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGGTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAActgct                                                                                                                                                                                    >  1:234339/1‑251 (MQ=255)
                                                                   cgATTAAACGCAGCTTCATACTGATTCCAGAAAACCACACCGTTCTGCCCGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGATCTGCTGGGc                                                                                                                                                                                <  2:79635/251‑1 (MQ=255)
                                                                            gattctccATACTGATTACAAAAACCCACACCGTGCTGCTCGCTGTCCGGCGTAATAAATTTTTTGCCATAACGGAGCCATGCGACGCTCGGAACTGATGGCGGTTTCCCCGATGCGGATGGTGCCTGGTTATTCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCAGCAACTGCTGACGATTGAATACGTGTTTGTATACCATTTTGTCGATGAACTGCTGGGCATTCGGGtt                                                                                                                                                                       <  2:282081/244‑1 (MQ=255)
                                                                            gCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATACCGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGACTGAGCGAGAATTTCCTGCAACTGCTGCCGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGtt                                                                                                                                                                       >  1:424203/1‑251 (MQ=255)
                                                                                  tCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGTCTCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGttatt                                                                                                                                                                    <  1:310109/248‑1 (MQ=255)
                                                                                                    aCCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCAATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTg                                                                                                                                               <  1:47733/251‑1 (MQ=255)
                                                                                                       acacCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGGGCCTGGTTATCCATCAGCCGCAGTACCGGATCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGAACGAAACCGGGTTTGTTCACCATTTTGGCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCAt                                                                                                                                            >  2:187245/1‑251 (MQ=255)
                                                                                                       acacCGTTCTGCACGCTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAAGCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCAt                                                                                                                                            <  2:98542/251‑1 (MQ=255)
                                                                                                          ccGTTCTGCATGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATcac                                                                                                                                         >  1:39338/1‑251 (MQ=255)
                                                                                                             ttCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGGGTTTGTTCCCCATTTTGTCGATGAACTGCTGGGCATTCGGGTTTTTAGCGAAATCACACCCCCTCTGCATccaatt                                                                                                                                      >  2:249953/1‑246 (MQ=255)
                                                                                                                    cGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGACTGGTTATACATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGAACGAAACCGTGTTTGTTCACCCTTTTGTCGATGGACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCAAACTGCATCACATTGtgctgc                                                                                                                               >  1:93146/1‑251 (MQ=255)
                                                                                                                      ttGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCAACGATGTGGGTGGGGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCAACATTTTGGCGATGAACTGCTGGGCATTCGGGGTATTAGCGAACTCGCCGCCCATCTGCACCACATTGTGCTGCgg                                                                                                                             >  2:54917/1‑251 (MQ=255)
                                                                                                                        gTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGATTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCt                                                                                                                           >  1:317527/1‑251 (MQ=255)
                                                                                                                               gTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGAGGTTGGTGCCTGGTTATCCATAAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGGGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGACCTGCTGGGCATTAGGGTTATTATCGAAATCGCCGCCCATCTGCATCCCCTTGTGCTGCGGCCCAAGc                                                                                                                      >  2:156981/1‑249 (MQ=255)
                                                                                                                                    aaaTTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAg                                                                                                                <  1:156912/250‑1 (MQ=255)
                                                                                                                                        tttttGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAgccgc                                                                                                            <  1:46675/250‑1 (MQ=255)
                                                                                                                                          tttGCGATAACGGAGCCATGCGCCGTTCGAACCTGATGGCGGTTTCCCCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑CGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCACCGCCCATCTGCATCCCATTGTGCTGCGGCTCAAGCAGGAAgccgcc                                                                                                           >  1:375640/1‑249 (MQ=255)
                                                                                                                                              cGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCCCCATTTTGTCGATGAACTGCTGGGCCTTCGGGTTATTAGCGAAATAGCCGCCCATCTGCCTCACATTGTTCTGCGGCTCAAGCAGGAAACCGCCCGACgg                                                                                                     >  2:244530/1‑251 (MQ=255)
                                                                                                                                                                       gACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCa                                                                            <  1:187245/251‑1 (MQ=255)
                                                                                                                                                                       gACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTTACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGtctc                                                                             >  1:135177/1‑250 (MQ=255)
                                                                                                                                                                            gATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCACCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAgg                                                                       <  1:145655/251‑1 (MQ=255)
                                                                                                                                                                                     ttCACCGATGTGGTTGGTGCCTGGTTATCCCTCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGGTCCGGTGGGCGTATCAGTATCAGTAGGTTTTGGCtt                                                              >  2:369139/1‑251 (MQ=255)
                                                                                                                                                                                       cACCGATGTTGTTGGTGCCTGGGTATCCATCATCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGATAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTtgc                                                            <  2:172968/251‑1 (MQ=255)
                                                                                                                                                                                        aCCGATGTGGTTGGTTCCTGGTTAGCCAGCAGCCGCATTACCGAATCC‑CGACGCTGCGCCTGATAGATAATTTCCTTCAACTGCGGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGATCTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTtgct                                                           <  2:183857/251‑1 (MQ=255)
                                                                                                                                                                                           gATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTgctgct                                                        <  1:369975/251‑1 (MQ=255)
                                                                                                                                                                                                     gtgacTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGgca                                              <  2:399023/247‑1 (MQ=255)
                                                                                                                                                                                                              tATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAacac                                      <  1:337512/250‑1 (MQ=255)
                                                                                                                                                                                                                 ccATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTTGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACaaaa                                  <  1:465380/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACGTTGTGCTGCGGCTCAAGCAGGACGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGATACAGGCGGCCAGCAACACaaa                                   >  2:298386/1‑250 (MQ=255)
                                                                                                                                                                                                                            aGTACCGAATCC‑AGACGCTTCGCCTGAGGGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCa                        >  1:182167/1‑250 (MQ=255)
                                                                                                                                                                                                                             gTACCGAATCC‑AGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGGCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCTGCCAGACGGCGGTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAAt                      <  2:39338/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaccc‑AGACGCTTCGCCTGAGAGATAATTTCCTTCAACTGCTGACGATCGAAACCGTGTTTTGTCCCCATTTTGTCGATGAACTGCTGGGCATTCGGGTGATTATCGAAATCGCCGCCCATCTGCAGCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCTTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGTGGCAATCAACACAAAAAGGGGAAGTAATGTTACa                <  2:416749/248‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCC‑CGACGCTTCGCCTTAGTGAGAATTTCCTGCAACTGCTGCTGATCGAAACCGTTGTTGTTCACCATTTTGTCGATGAACTGTTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCAGCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACa                <  2:303585/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGACGCTTCGGCTGAGAGAGAAATTCCTGCAACTGCTGACTAGCGAAACCGGGGTTGTTCACCATTTTGTCGATGATCTGCGGGGCATTTGGGTTATTAGCGAAATCGCCGCCCATCTGCATCTCATGGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCTTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGTAAGCAATGTTACATAac            <  2:87136/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTCTCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAacg           >  2:188830/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGc       >  2:303486/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCt      <  1:370348/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGCCTGAGAGAGTATTTCCTGCAACTGCGGACGATCGAAACCGTTGTTGTTCATCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATGGGGCTGCGGCTCAAGCAGGAAGCCTCTAGACGGCGTTCCGGTGGTCGTAGCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCACCACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGa   <  2:375640/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGaa  <  1:342237/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCC‑AGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGAA  >  NC_000913/2825082‑2825573

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