Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 521,941 T→G 5.3% L509R (CTC→CGC)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913521,9410TG 5.3% 274.7 / 4.5 75L509R (CTC→CGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base G (2/2);  total (37/38)
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.98e-01

TTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGA  >  NC_000913/521794‑522085
                                                                                                                                                   |                                                                                                                                                
tttCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCt                                                                                                                                                  <  4:1019310/148‑1 (MQ=255)
         ccTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcg                                                                                                                                                >  3:326071/1‑141 (MQ=255)
         ccTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcg                                                                                                                                                <  4:326071/141‑1 (MQ=255)
           tttATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATc                                                                                                                                       >  2:160355/1‑148 (MQ=255)
                     tgGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTccc                                                                                                                             <  2:549254/148‑1 (MQ=255)
                       gCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCgg                                                                                                                           >  4:227225/1‑148 (MQ=255)
                        cACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACAt                                                                                                                                        <  1:1063210/134‑1 (MQ=255)
                        cACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACAt                                                                                                                                        <  3:320221/134‑1 (MQ=255)
                        cACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACAt                                                                                                                                        >  4:320221/1‑134 (MQ=255)
                        cACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACAt                                                                                                                                        >  2:1063210/1‑134 (MQ=255)
                             ctgctgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGt                                                                                                                     <  3:227225/148‑1 (MQ=255)
                               gctgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTccc                                                                                                                             <  2:483434/138‑1 (MQ=255)
                               gctgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTccc                                                                                                                             >  1:483434/1‑138 (MQ=255)
                                 tgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGtt                                                                                                                    >  1:1089658/1‑145 (MQ=255)
                                 tgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGtt                                                                                                                    <  2:1089658/145‑1 (MQ=255)
                                            ggCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcg                                                                                                                                                >  1:782303/1‑106 (MQ=255)
                                            ggCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcg                                                                                                                                                <  2:782303/106‑1 (MQ=255)
                                                    gCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGaa                                                                                                                                           >  2:1119435/1‑103 (MQ=255)
                                                    gCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGaa                                                                                                                                           <  1:1119435/103‑1 (MQ=255)
                                                       cGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCggt                                                                                                           <  4:622794/132‑1 (MQ=255)
                                                       cGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCgga                                                                                                           >  3:622794/1‑131 (MQ=255)
                                                            gCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTg                                                                                           >  3:47116/1‑143 (MQ=255)
                                                            gCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTg                                                                                           <  4:47116/143‑1 (MQ=255)
                                                                    cagcagcTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   >  3:184243/1‑143 (MQ=255)
                                                                    cagcagcTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   <  4:184243/143‑1 (MQ=255)
                                                                    cagcagcTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACg                                                                              <  1:160355/148‑1 (MQ=255)
                                                                                  aGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 <  3:381340/147‑1 (MQ=255)
                                                                                  aGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 >  4:381340/1‑147 (MQ=255)
                                                                                   gAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATaa                                                                          <  3:421114/137‑1 (MQ=255)
                                                                                   gAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATaa                                                                          >  4:421114/1‑137 (MQ=255)
                                                                                    aaaGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGa                                                                                                                                    <  4:427043/78‑1 (MQ=255)
                                                                                    aaaGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGa                                                                                                                                    >  3:427043/1‑78 (MQ=255)
                                                                                          cGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGgcg                                                                                                 <  3:3023/107‑1 (MQ=255)
                                                                                          cGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGgcg                                                                                                 >  4:3023/1‑107 (MQ=255)
                                                                                                   ttttAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGTGCGCGGCTGACGGCg                                                                                     >  4:837146/1‑110 (MQ=255)
                                                                                                   ttttAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGTGCGCGGCTGACGGCg                                                                                     <  3:837146/110‑1 (MQ=255)
                                                                                                          taaCATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        >  2:999103/1‑116 (MQ=255)
                                                                                                          taaCATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        <  1:999103/116‑1 (MQ=255)
                                                                                                            aCATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGa                                                     <  4:394777/133‑1 (MQ=255)
                                                                                                            aCATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGa                                                     >  3:394777/1‑133 (MQ=255)
                                                                                                               tCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGa                                                                                          <  1:1124663/93‑1 (MQ=255)
                                                                                                               tCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGa                                                                                          >  2:1124663/1‑93 (MQ=255)
                                                                                                               tCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCg                                                                                     <  3:1140071/98‑1 (MQ=255)
                                                                                                               tCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCg                                                                                     >  4:1140071/1‑98 (MQ=255)
                                                                                                                 gCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTaa                                                                                >  1:560689/1‑101 (MQ=255)
                                                                                                                 gCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTaa                                                                                <  2:560689/101‑1 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   <  2:199036/94‑1 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   >  1:199036/1‑94 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   >  3:683257/1‑142 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   >  1:21122/1‑142 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   <  4:683257/142‑1 (MQ=255)
                                                                                                                     caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   <  2:21122/142‑1 (MQ=255)
                                                                                                                       gAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACg                                                                              <  4:494027/97‑1 (MQ=255)
                                                                                                                       gAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACg                                                                              >  3:494027/1‑97 (MQ=255)
                                                                                                                        aaCAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGaa                                                                >  2:309838/1‑110 (MQ=255)
                                                                                                                        aaCAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGaa                                                                <  1:309838/110‑1 (MQ=255)
                                                                                                                          cAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGa                                                                     <  4:839079/103‑1 (MQ=255)
                                                                                                                          cAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGa                                                                     >  3:839079/1‑103 (MQ=255)
                                                                                                                           aGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCgg                                                                                                            >  4:626955/1‑63 (MQ=255)
                                                                                                                           aGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCgg                                                                                                            <  3:626955/63‑1 (MQ=255)
                                                                                                                                 cACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCa                 >  4:889131/1‑148 (MQ=255)
                                                                                                                                   ccGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTa                                                            <  3:21981/103‑1 (MQ=255)
                                                                                                                                   ccGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTa                                                            >  4:21981/1‑103 (MQ=255)
                                                                                                                                    cGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGgcg                                              <  2:238900/116‑1 (MQ=255)
                                                                                                                                    cGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGgcg                                              >  1:238900/1‑116 (MQ=255)
                                                                                                                                          aaGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             >  1:669687/1‑127 (MQ=255)
                                                                                                                                          aaGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             <  2:669687/127‑1 (MQ=255)
                                                                                                                                           gatgTTCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCtt                         <  4:675857/127‑1 (MQ=255)
                                                                                                                                           gatgTTCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCtt                         >  3:675857/4‑130 (MQ=255)
                                                                                                                                               ttCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATg                                                   >  4:642792/1‑100 (MQ=255)
                                                                                                                                               ttCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATg                                                   <  3:642792/100‑1 (MQ=255)
                                                                                                                                               ttCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        >  1:840742/1‑79 (MQ=255)
                                                                                                                                               ttCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        <  2:840742/79‑1 (MQ=255)
                                                                                                                                                tCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGa  <  3:889131/148‑1 (MQ=255)
                                                                                                                                                 ccTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             >  2:61387/1‑120 (MQ=255)
                                                                                                                                                 ccTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             <  1:61387/120‑1 (MQ=255)
                                                                                                                                                   tCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 >  1:383139/1‑82 (MQ=255)
                                                                                                                                                   tCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 <  2:383139/82‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGA  >  NC_000913/521794‑522085

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