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→TCTmltBmembrane‑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

CCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGAA  >  NC_000913/2825082‑2825573
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
ccAGCGGGGTTCAACGCCGATAATCCCGACGATAATGTCCGGCGGTACTCCATCCACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTGGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGCCCTGATGTCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTgagaga                                                                                                                                                                                                                                                   <  1:283278/251‑1 (MQ=255)
         tccaacccctaTAACCCCGACGATAATTACCGGCGGTACTCCATACCCCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTGCGGTCCTTATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGGGATAGAATTTcc                                                                                                                                                                                                                                           <  2:218032/240‑1 (MQ=255)
          gcAACGCCGATAATCCCGAAGATAATTTCCGGTGGTACGCCAGACACCTGACACGCGCGATTCAACGCATTTTCATACTGATTACAGAAAACCACACCGTTCTGCATGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGTCGTTAGGACCTGATGGCGGTTTCACCGATGTTGTTGGTGCCTGGTTATCCATCAGCCGCAGTCCCGAATCCAGACGGTTCGCCTGAGAGAGAACTTCCt                                                                                                                                                                                                                                          <  2:156433/249‑1 (MQ=255)
                ccGATACTCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACTCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCTCGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAActg                                                                                                                                                                                                                                   <  2:143272/251‑1 (MQ=255)
                ccGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGCGAGAGAATTTCCTGCAActg                                                                                                                                                                                                                                   >  1:87136/1‑251 (MQ=255)
                  gATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAATCCACACCGTTCTGCACGTTTTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAActgct                                                                                                                                                                                                                                 <  1:231227/251‑1 (MQ=255)
                       caccGAAGATAATTCCCCGGGGTACTCCATCCACCTGCCCCGCGCGATTCAACTCATCTTCATGCTGATTCCAGAAATCCGCGCCGTTCTGCACGGGGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGTCCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGTTTATCCATCAGCCGCAGTACCGAATGCAGGCTCTTCTCCTGAGAGAGAATGTCCTGCAACTGCTGACGa                                                                                                                                                                                                                            <  2:427848/249‑1 (MQ=255)
                               aTAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAAcc                                                                                                                                                                                                                    <  1:128901/251‑1 (MQ=255)
                                tAATTTCCGGCGGTACTCCATACCCCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGGGCCTGGTTATCCATCAGCCGCAGTACCGAATCCCGACGCTTCGCCTGAGGGGGACTTTCCTGCAACTGCGGACGATCGAAAc                                                                                                                                                                                                                     >  2:342237/1‑249 (MQ=255)
                                tAATTTCCGGCGGTACTCCATAAACCTGCCACGCGCGATTCCACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAACTTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGGTTTCCATCCGCAGCCGGACCGCATCCCAGCGCTTCGCCTGAGAGAGAATTTCCTCCCACTACTGACgcgagcacccc                                                                                                                                                                                                                   >  2:190406/1‑241 (MQ=255)
                                 aaTTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCgg                                                                                                                                                                                                                  >  1:113982/1‑250 (MQ=255)
                                      ccGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGtgttt                                                                                                                                                                                                              <  1:52639/250‑1 (MQ=255)
                                         gcggTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTc                                                                                                                                                                                                          <  1:110808/251‑1 (MQ=255)
                                          cggTACTCCGTACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCAATCAGCCGCAGTACCGAATCCAAACGCTTCGCCTGCGTGAGAATTTCCTTCAACTGCTGACGATCGAAAACGTGTTTGTTc                                                                                                                                                                                                          >  1:342121/1‑250 (MQ=255)
                                                 ccATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGAAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGAAACTGCTGAAGATCGAAACCGAGTTTGTTCACCAtttt                                                                                                                                                                                                  >  1:223846/1‑251 (MQ=255)
                                                  cATACCCCTGTCACGCGCGATTCAACGCATCTTCATACGGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTGTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTg                                                                                                                                                                                                 <  1:22653/251‑1 (MQ=255)
                                                  cATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGCTGGGGATGGAGCCTGGTTATCCATCAACCGCAGTACCGAATCCAGACGCTTCGCCTGCGAGAGAATTTCCTACACCTTCTGACGAACCAAACCGGGTTTGTTCACCCttttt                                                                                                                                                                                                 >  2:328583/1‑250 (MQ=255)
                                                            ccACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAAct                                                                                                                                                                                       >  1:399023/1‑251 (MQ=255)
                                                               cgcgcgATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGGTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAActgct                                                                                                                                                                                    >  1:234339/1‑251 (MQ=255)
                                                                   cgATTAAACGCAGCTTCATACTGATTCCAGAAAACCACACCGTTCTGCCCGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGATCTGCTGGGc                                                                                                                                                                                <  2:79635/251‑1 (MQ=255)
                                                                            gattctccATACTGATTACAAAAACCCACACCGTGCTGCTCGCTGTCCGGCGTAATAAATTTTTTGCCATAACGGAGCCATGCGACGCTCGGAACTGATGGCGGTTTCCCCGATGCGGATGGTGCCTGGTTATTCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCAGCAACTGCTGACGATTGAATACGTGTTTGTATACCATTTTGTCGATGAACTGCTGGGCATTCGGGtt                                                                                                                                                                       <  2:282081/244‑1 (MQ=255)
                                                                            gCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATACCGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGACTGAGCGAGAATTTCCTGCAACTGCTGCCGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGtt                                                                                                                                                                       >  1:424203/1‑251 (MQ=255)
                                                                                  tCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGTCTCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGttatt                                                                                                                                                                    <  1:310109/248‑1 (MQ=255)
                                                                                                    aCCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCAATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTg                                                                                                                                               <  1:47733/251‑1 (MQ=255)
                                                                                                       acacCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGGGCCTGGTTATCCATCAGCCGCAGTACCGGATCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGAACGAAACCGGGTTTGTTCACCATTTTGGCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCAt                                                                                                                                            >  2:187245/1‑251 (MQ=255)
                                                                                                       acacCGTTCTGCACGCTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAAGCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCAt                                                                                                                                            <  2:98542/251‑1 (MQ=255)
                                                                                                          ccGTTCTGCATGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATcac                                                                                                                                         >  1:39338/1‑251 (MQ=255)
                                                                                                             ttCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGGGTTTGTTCCCCATTTTGTCGATGAACTGCTGGGCATTCGGGTTTTTAGCGAAATCACACCCCCTCTGCATccaatt                                                                                                                                      >  2:249953/1‑246 (MQ=255)
                                                                                                                    cGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGACTGGTTATACATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGAACGAAACCGTGTTTGTTCACCCTTTTGTCGATGGACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCAAACTGCATCACATTGtgctgc                                                                                                                               >  1:93146/1‑251 (MQ=255)
                                                                                                                      ttGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCAACGATGTGGGTGGGGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCAACATTTTGGCGATGAACTGCTGGGCATTCGGGGTATTAGCGAACTCGCCGCCCATCTGCACCACATTGTGCTGCgg                                                                                                                             >  2:54917/1‑251 (MQ=255)
                                                                                                                        gTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGATTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCt                                                                                                                           >  1:317527/1‑251 (MQ=255)
                                                                                                                               gTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGAGGTTGGTGCCTGGTTATCCATAAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGGGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGACCTGCTGGGCATTAGGGTTATTATCGAAATCGCCGCCCATCTGCATCCCCTTGTGCTGCGGCCCAAGc                                                                                                                      >  2:156981/1‑249 (MQ=255)
                                                                                                                                    aaaTTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAg                                                                                                                <  1:156912/250‑1 (MQ=255)
                                                                                                                                        tttttGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAgccgc                                                                                                            <  1:46675/250‑1 (MQ=255)
                                                                                                                                          tttGCGATAACGGAGCCATGCGCCGTTCGAACCTGATGGCGGTTTCCCCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCCGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCACCGCCCATCTGCATCCCATTGTGCTGCGGCTCAAGCAGGAAgccgcc                                                                                                           >  1:375640/1‑249 (MQ=255)
                                                                                                                                              cGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCCCCATTTTGTCGATGAACTGCTGGGCCTTCGGGTTATTAGCGAAATAGCCGCCCATCTGCCTCACATTGTTCTGCGGCTCAAGCAGGAAACCGCCCGACgg                                                                                                     >  2:244530/1‑251 (MQ=255)
                                                                                                                                                                       gACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCa                                                                            <  1:187245/251‑1 (MQ=255)
                                                                                                                                                                       gACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTTACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGtctc                                                                             >  1:135177/1‑250 (MQ=255)
                                                                                                                                                                            gATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCACCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAgg                                                                       <  1:145655/251‑1 (MQ=255)
                                                                                                                                                                                     ttCACCGATGTGGTTGGTGCCTGGTTATCCCTCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGGTCCGGTGGGCGTATCAGTATCAGTAGGTTTTGGCtt                                                              >  2:369139/1‑251 (MQ=255)
                                                                                                                                                                                       cACCGATGTTGTTGGTGCCTGGGTATCCATCATCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGATAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTtgc                                                            <  2:172968/251‑1 (MQ=255)
                                                                                                                                                                                        aCCGATGTGGTTGGTTCCTGGTTAGCCAGCAGCCGCATTACCGAATCCCGACGCTGCGCCTGATAGATAATTTCCTTCAACTGCGGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGATCTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTtgct                                                           <  2:183857/251‑1 (MQ=255)
                                                                                                                                                                                           gATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTgctgct                                                        <  1:369975/251‑1 (MQ=255)
                                                                                                                                                                                                     gtgacTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGgca                                              <  2:399023/247‑1 (MQ=255)
                                                                                                                                                                                                              tATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAacac                                      <  1:337512/250‑1 (MQ=255)
                                                                                                                                                                                                                 ccATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTTGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACaaaa                                  <  1:465380/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACGTTGTGCTGCGGCTCAAGCAGGACGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGATACAGGCGGCCAGCAACACaaa                                   >  2:298386/1‑250 (MQ=255)
                                                                                                                                                                                                                            aGTACCGAATCCAGACGCTTCGCCTGAGGGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCa                        >  1:182167/1‑250 (MQ=255)
                                                                                                                                                                                                                             gTACCGAATCCAGACGCTTCGCCTGAGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGGCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCTGCCAGACGGCGGTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAAt                      <  2:39338/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aacccAGACGCTTCGCCTGAGAGATAATTTCCTTCAACTGCTGACGATCGAAACCGTGTTTTGTCCCCATTTTGTCGATGAACTGCTGGGCATTCGGGTGATTATCGAAATCGCCGCCCATCTGCAGCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCTTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGTGGCAATCAACACAAAAAGGGGAAGTAATGTTACa                <  2:416749/248‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTCCCGACGCTTCGCCTTAGTGAGAATTTCCTGCAACTGCTGCTGATCGAAACCGTTGTTGTTCACCATTTTGTCGATGAACTGTTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCAGCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACa                <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAACCACACCGTTCTGCACGTTGTCCGGCGTAATAAATTTTTTGCGATAACGGAGCCATGCGCCGTTCGGACCTGATGGCGGTTTCACCGATGTGGTTGGTGCCTGGTTATCCATCAGCCGCAGTACCGAATCCAGACGCTTCGCCTGGGAGAGAATTTCCTGCAACTGCTGACGATCGAAACCGTGTTTGTTCACCATTTTGTCGATGAACTGCTGGGCATTCGGGTTATTAGCGAAATCGCCGCCCATCTGCATCACATTGTGCTGCGGCTCAAGCAGGAAGCCGCCAGACGGCGTTCCGGTGGTCGTATCAGTCTCAGTAGGTTTTGGCTTGCTGCTACAGGCGGCAAGCAACACAAAAAGGGGAAGCAATGTTACATAACGACGCTTGAA  >  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: