Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,903,534 G→A 5.3% V151V (GTG→GTA manY → mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,903,5340GA 5.3% 309.8 / 5.0 94V151V (GTG→GTAmanYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (42/47);  new base A (2/3);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTAT  >  NC_000913/1903391‑1903674
                                                                                                                                               |                                                                                                                                            
gCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTa                                                                                                                                          <  3:348620/148‑1 (MQ=255)
            gTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCa                                                                                                                              >  3:26511/1‑148 (MQ=255)
                 gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGtt                                                                                                                         >  4:362146/1‑148 (MQ=255)
                  aCCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGt                                                                                                                          <  8:235992/146‑1 (MQ=255)
                      tcatcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGc                                                                                                                                   <  8:92214/133‑1 (MQ=255)
                      tcatcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGc                                                                                                                                   >  7:92214/1‑133 (MQ=255)
                       catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACGTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGct                                                                                                                    <  8:431651/147‑1 (MQ=255)
                       catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctg                                                                                                                   >  2:343330/1‑148 (MQ=255)
                       catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctg                                                                                                                   >  5:435986/1‑148 (MQ=255)
                        atcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCa                                                                                                                              <  4:228350/136‑1 (MQ=255)
                        atcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCa                                                                                                                              >  3:228350/1‑136 (MQ=255)
                        atcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                                                      <  6:128345/144‑1 (MQ=255)
                        atcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                                                      >  5:128345/1‑144 (MQ=255)
                          cgttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtct                                                                                                                >  7:131922/1‑148 (MQ=255)
                            ttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctgt                                                                                                              <  2:330259/148‑1 (MQ=255)
                             tcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTAGCTATTCCGGCCGTCATCGTTGCGCTGTCTGtt                                                                                                             >  8:196768/1‑148 (MQ=255)
                              cgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGtt                                                                                                             <  4:146115/147‑1 (MQ=255)
                               gtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCActgt                                                                                                                  >  4:249677/1‑141 (MQ=255)
                               gtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCActgt                                                                                                                  <  3:249677/141‑1 (MQ=255)
                                tACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGt                                                                                                          <  4:33660/148‑1 (MQ=255)
                                 aCTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCAt                                                                                                                             <  4:301831/128‑1 (MQ=255)
                                 aCTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCAt                                                                                                                             >  3:301831/1‑128 (MQ=255)
                                    aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAc                                                                                                        >  7:367148/1‑146 (MQ=255)
                                    aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAc                                                                                                        >  4:73771/1‑146 (MQ=255)
                                    aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAc                                                                                                        >  6:107825/1‑146 (MQ=255)
                                     ttACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAg                                                                                                     >  3:146871/1‑148 (MQ=255)
                                        ccGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                                                   >  6:44767/1‑147 (MQ=255)
                                        ccGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGa                                                                                                  <  7:132018/148‑1 (MQ=255)
                                        ccGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGa                                                                                                  <  3:362146/148‑1 (MQ=255)
                                         cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                                                   >  7:442175/1‑146 (MQ=255)
                                         cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                                                   >  4:382146/1‑146 (MQ=255)
                                          gTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGa                                                                                                  <  8:442175/146‑1 (MQ=255)
                                          gTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAg                                                                                                >  6:218807/1‑148 (MQ=255)
                                          gTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGCCACCCTGACAGCGATTTCCTGGATCCACGATTCTTCTCTGTTCCTGCAAGCAAAGCGTGTAGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGTAg                                                                                                >  1:192530/1‑148 (MQ=255)
                                            tGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTg                                                                                                                        >  2:394144/1‑122 (MQ=255)
                                            tGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTg                                                                                                                        <  1:394144/122‑1 (MQ=255)
                                              cTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                             <  1:343330/147‑1 (MQ=255)
                                                ttCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAg                                                                                           >  8:50451/1‑147 (MQ=255)
                                                   cAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                                                        <  8:367148/147‑1 (MQ=255)
                                                    aGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctg                                                                                                               <  1:286936/123‑1 (MQ=255)
                                                    aGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctg                                                                                                               >  2:286936/1‑123 (MQ=255)
                                                    aGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACAt                                                                                      <  5:218807/148‑1 (MQ=255)
                                                     gCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATg                                                                                     <  7:50451/148‑1 (MQ=255)
                                                      cACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGc                                                                                    <  4:146871/148‑1 (MQ=255)
                                                      cACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCATGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGc                                                                                    >  7:328633/1‑148 (MQ=255)
                                                         gCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctgt                                                                                                              >  2:94275/1‑119 (MQ=255)
                                                         gCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctgt                                                                                                              <  1:94275/119‑1 (MQ=255)
                                                         gCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                                                                 <  8:328243/148‑1 (MQ=255)
                                                          cTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGaa                                                                                >  4:358752/1‑148 (MQ=255)
                                                            gCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctg                                                                                                               >  2:187100/1‑115 (MQ=255)
                                                            gCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctg                                                                                                               <  1:187100/115‑1 (MQ=255)
                                                              ggATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCg                                                                            <  4:26511/148‑1 (MQ=255)
                                                                     gctgctGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCg                                                                      >  4:137934/1‑147 (MQ=255)
                                                                           gATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGaa                                                                                         >  3:46507/1‑122 (MQ=255)
                                                                           gATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGaa                                                                                         <  4:46507/122‑1 (MQ=255)
                                                                            aTAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                                                                 >  6:400832/1‑129 (MQ=255)
                                                                            aTAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                                                                 <  5:400832/129‑1 (MQ=255)
                                                                            aTAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtg                                                                <  8:131922/146‑1 (MQ=255)
                                                                            aTAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTAGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                              <  7:196768/148‑1 (MQ=255)
                                                                               aCGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                            <  5:44767/147‑1 (MQ=255)
                                                                                   cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                                                       <  3:137934/148‑1 (MQ=255)
                                                                                      ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                             <  6:237477/107‑1 (MQ=255)
                                                                                      ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                             >  5:237477/1‑107 (MQ=255)
                                                                                      ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGt                                                    <  6:435986/148‑1 (MQ=255)
                                                                                        tGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCt                                                  >  4:253923/1‑148 (MQ=255)
                                                                                            aGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCa                                                         <  3:331159/137‑1 (MQ=255)
                                                                                            aGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGaa                                               >  1:243349/1‑147 (MQ=255)
                                                                                             gCGATTTCCTGGATGCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTAGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGATGTACAGACCATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAata                                             <  2:192530/148‑1 (MQ=255)
                                                                                             gCGATTTCCAGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGAGAGGCTATTCCGGCCGTCATCGTAGCCCTCTCTGTTGGTACCAGCGAAGTAAAGAACAACCTGAATGCGATTCAGGAAGTGGTGACCaa                                                        >  4:331159/1‑137 (MQ=255)
                                                                                               gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                           >  6:146151/1‑148 (MQ=255)
                                                                                                  ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc                                        >  8:123606/1‑148 (MQ=255)
                                                                                                        gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCATGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc                                  <  8:328633/148‑1 (MQ=255)
                                                                                                                tttCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                             <  1:273886/81‑1 (MQ=255)
                                                                                                                tttCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgt                          >  7:379131/1‑148 (MQ=255)
                                                                                                                tttCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgt                          <  3:73771/148‑1 (MQ=255)
                                                                                                                   cttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtgg                        <  7:215695/147‑1 (MQ=255)
                                                                                                                     tctcTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTCCGCTCTCTGTTGGTACCAGCGAAGTACATAACATGCTGAATGCGATTCCGGAAGTGGTGACAAATGGTCTGAAAATCGCCTGTGGCATGATCGtggtgg                     >  8:344430/1‑148 (MQ=255)
                                                                                                                      ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                       <  8:379131/145‑1 (MQ=255)
                                                                                                                       tctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTggtggt                    <  5:146151/147‑1 (MQ=255)
                                                                                                                           ttCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGt                                                    >  2:333376/1‑111 (MQ=255)
                                                                                                                           ttCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGt                                                    <  1:333376/111‑1 (MQ=255)
                                                                                                                           ttCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt               <  5:107825/148‑1 (MQ=255)
                                                                                                                             ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa              <  3:253923/147‑1 (MQ=255)
                                                                                                                             ccTGCAAGCAATGCGTGTAGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa              <  3:324229/147‑1 (MQ=255)
                                                                                                                              cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATg            <  3:358752/148‑1 (MQ=255)
                                                                                                                               tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGc           >  2:416952/1‑148 (MQ=255)
                                                                                                                                 caagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGa         <  1:416952/148‑1 (MQ=255)
                                                                                                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCt                                                  <  4:67213/100‑1 (MQ=255)
                                                                                                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCt                                                  >  3:67213/1‑100 (MQ=255)
                                                                                                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt               <  6:283403/135‑1 (MQ=255)
                                                                                                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt               >  5:283403/1‑135 (MQ=255)
                                                                                                                                         gCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCGGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCTATGGTTAt  >  6:12497/1‑147 (MQ=255)
                                                                                                                                          cGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggtt                   >  4:402823/1‑129 (MQ=255)
                                                                                                                                          cGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggtt                   <  3:402823/129‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
GCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTAT  >  NC_000913/1903391‑1903674

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