Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 521,944 C→A 5.5% A510E (GCG→GAG)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913521,9440CA 5.5% 260.7 / 4.3 73A510E (GCG→GAG) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base A (2/2);  total (36/37)
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.75e-01

TTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGA  >  NC_000913/521805‑522085
                                                                                                                                           |                                                                                                                                             
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)
                                         gCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGaa                                                                                                                                           <  1:1119435/103‑1 (MQ=255)
                                         gCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGaa                                                                                                                                           >  2:1119435/1‑103 (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                                                                                   <  4:184243/143‑1 (MQ=255)
                                                         cagcagcTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   >  3:184243/1‑143 (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                                                                                                                                    >  3:427043/1‑78 (MQ=255)
                                                                         aaaGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGa                                                                                                                                    <  4:427043/78‑1 (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                                                                                <  2:560689/101‑1 (MQ=255)
                                                                                                      gCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTaa                                                                                >  1:560689/1‑101 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   >  1:199036/1‑94 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAt                                                                                   <  2:199036/94‑1 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   >  3:683257/1‑142 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   >  1:21122/1‑142 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   <  2:21122/142‑1 (MQ=255)
                                                                                                          caGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGa                                   <  4:683257/142‑1 (MQ=255)
                                                                                                            gAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACg                                                                              <  4:494027/97‑1 (MQ=255)
                                                                                                            gAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACg                                                                              >  3:494027/1‑97 (MQ=255)
                                                                                                             aaCAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGaa                                                                <  1:309838/110‑1 (MQ=255)
                                                                                                             aaCAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGaa                                                                >  2:309838/1‑110 (MQ=255)
                                                                                                               cAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGa                                                                     >  3:839079/1‑103 (MQ=255)
                                                                                                               cAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGa                                                                     <  4:839079/103‑1 (MQ=255)
                                                                                                                aGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCgg                                                                                                            <  3:626955/63‑1 (MQ=255)
                                                                                                                aGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCgg                                                                                                            >  4:626955/1‑63 (MQ=255)
                                                                                                                      cACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCa                 >  4:889131/1‑148 (MQ=255)
                                                                                                                        ccGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTa                                                            >  4:21981/1‑103 (MQ=255)
                                                                                                                        ccGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTa                                                            <  3:21981/103‑1 (MQ=255)
                                                                                                                         cGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGgcg                                              <  2:238900/116‑1 (MQ=255)
                                                                                                                         cGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGgcg                                              >  1:238900/1‑116 (MQ=255)
                                                                                                                               aaGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             <  2:669687/127‑1 (MQ=255)
                                                                                                                               aaGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             >  1:669687/1‑127 (MQ=255)
                                                                                                                                gatgTTCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCtt                         <  4:675857/127‑1 (MQ=255)
                                                                                                                                gatgTTCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCtt                         >  3:675857/4‑130 (MQ=255)
                                                                                                                                    ttCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATg                                                   <  3:642792/100‑1 (MQ=255)
                                                                                                                                    ttCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATg                                                   >  4:642792/1‑100 (MQ=255)
                                                                                                                                    ttCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        >  1:840742/1‑79 (MQ=255)
                                                                                                                                    ttCCGCGAGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGc                                                                        <  2:840742/79‑1 (MQ=255)
                                                                                                                                     tCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGa  <  3:889131/148‑1 (MQ=255)
                                                                                                                                      ccTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             <  1:61387/120‑1 (MQ=255)
                                                                                                                                      ccTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa                             >  2:61387/1‑120 (MQ=255)
                                                                                                                                        tCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 >  1:383139/1‑82 (MQ=255)
                                                                                                                                        tCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGa                                                                 <  2:383139/82‑1 (MQ=255)
                                                                                                                                         cgcgGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTaa                                         <  3:918469/105‑1 (MQ=255)
                                                                                                                                         cgcgGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTaa                                         >  4:918469/1‑105 (MQ=255)
                                                                                                                                           |                                                                                                                                             
TTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGA  >  NC_000913/521805‑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: