Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,238,817 A→T 5.1% I158N (ATC→AAC)  mglA ← methyl‑galactoside ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,238,8170AT 5.1% 280.7 / 4.7 78I158N (ATC→AAC) mglAmethyl‑galactoside ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (43/31);  new base T (2/2);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

TGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCCACAT  >  NC_000913/2238670‑2238965
                                                                                                                                                   |                                                                                                                                                    
tGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGat                                                                                                                                                     <  4:504760/149‑1 (MQ=255)
tGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGat                                                                                                                                                     <  4:285405/149‑1 (MQ=255)
tGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGat                                                                                                                                                     <  4:429559/149‑1 (MQ=255)
 gCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGatc                                                                                                                                                    <  2:82717/149‑1 (MQ=255)
   gagaTATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGatcat                                                                                                                                                  >  4:258993/1‑149 (MQ=255)
   gagaTATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGatcat                                                                                                                                                  <  1:206032/149‑1 (MQ=255)
    agaTATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGatcatc                                                                                                                                                 <  4:198782/149‑1 (MQ=255)
    agaTATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGatcatc                                                                                                                                                 <  2:143277/149‑1 (MQ=255)
      atatAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTg                                                                                                                                               >  4:314231/1‑149 (MQ=255)
            aCAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTg                                                                                                                                         >  2:544538/1‑149 (MQ=255)
             cAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGc                                                                                                                                        >  2:674701/1‑149 (MQ=255)
                 aCCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGaaa                                                                                                                                    <  4:138434/149‑1 (MQ=255)
                     cAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGa                                                                                                                                <  3:412298/149‑1 (MQ=255)
                      aGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGAt                                                                                                                               <  3:258993/149‑1 (MQ=255)
                             ctctTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGc                                                                                                                        <  2:281435/149‑1 (MQ=255)
                               ctTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCg                                                                                                                      <  3:417536/149‑1 (MQ=255)
                                  ttAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGAcg                                                                                                                   >  1:634770/1‑149 (MQ=255)
                                           gAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCAcgc                                                                                                          >  2:413859/1‑149 (MQ=255)
                                           gAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCAcgc                                                                                                          >  3:439214/1‑149 (MQ=255)
                                           gAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCAcgc                                                                                                          >  3:99529/1‑149 (MQ=255)
                                              taataGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGTTCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGa                                                                                                       >  2:595025/1‑149 (MQ=255)
                                              taataGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGTTCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGa                                                                                                       >  4:595583/1‑149 (MQ=255)
                                               aataGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGAt                                                                                                      >  1:710298/1‑149 (MQ=255)
                                               aataGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGAt                                                                                                      >  4:557090/1‑149 (MQ=255)
                                                ataGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATc                                                                                                     >  1:14183/1‑149 (MQ=255)
                                                 taGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGc                                                                                                                        >  2:197331/1‑129 (MQ=255)
                                                 taGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGc                                                                                                                        <  1:197331/129‑1 (MQ=255)
                                                  aGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGa                                                                                                   >  3:33900/1‑149 (MQ=255)
                                                     gAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGatat                                                                                                <  2:14183/149‑1 (MQ=255)
                                                     gAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGatat                                                                                                <  1:544538/149‑1 (MQ=255)
                                                         aGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAat                                                                                            >  3:112961/1‑149 (MQ=255)
                                                             gatTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAAt                                                                                                                            <  1:620663/113‑1 (MQ=255)
                                                             gatTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAAt                                                                                                                            >  2:620663/1‑113 (MQ=255)
                                                              atTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCa                                                                                       >  2:287691/1‑149 (MQ=255)
                                                              atTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCa                                                                                       >  4:402572/1‑149 (MQ=255)
                                                                 gACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGtt                                                                                    <  3:314231/149‑1 (MQ=255)
                                                                    ctctTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGt                                                                                 >  3:260686/1‑149 (MQ=255)
                                                                     tctTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTc                                                                                <  1:755717/149‑1 (MQ=255)
                                                                            cGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATc                                                                         <  4:439214/149‑1 (MQ=255)
                                                                            cGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATc                                                                         <  4:99529/149‑1 (MQ=255)
                                                                             ggTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCg                                                                        >  2:44769/1‑149 (MQ=255)
                                                                              gTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGc                                                                       >  1:90611/1‑149 (MQ=255)
                                                                              gTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGc                                                                       <  2:15885/149‑1 (MQ=255)
                                                                                tAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCtt                                                                     >  2:474923/1‑149 (MQ=255)
                                                                                    gaagaaGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCg                                                                        <  1:713620/142‑1 (MQ=255)
                                                                                    gaagaaGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCg                                                                        >  2:713620/1‑142 (MQ=255)
                                                                                    gaagaaGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGt                                                                 >  4:407525/1‑149 (MQ=255)
                                                                                      agaagTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGTTCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGttt                                                               <  1:595025/149‑1 (MQ=255)
                                                                                      agaagTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGTTCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGttt                                                               <  3:595583/149‑1 (MQ=255)
                                                                                       gaagTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTc                                                              >  4:155341/1‑149 (MQ=255)
                                                                                         agTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTcgc                                                            <  2:710298/149‑1 (MQ=255)
                                                                                              gTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTAc                                                       >  4:594366/1‑149 (MQ=255)
                                                                                                  aTCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACAttt                                                   >  3:584288/1‑149 (MQ=255)
                                                                                                  aTCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACAttt                                                   >  4:537582/1‑149 (MQ=255)
                                                                                                   tCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACAtttt                                                  <  3:557090/149‑1 (MQ=255)
                                                                                                   tCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACAtttt                                                  >  1:317101/1‑149 (MQ=255)
                                                                                                     cATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGt                                                >  4:113864/1‑149 (MQ=255)
                                                                                                      aTAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTc                                               >  1:156561/1‑149 (MQ=255)
                                                                                                      aTAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTc                                               >  2:231480/1‑149 (MQ=255)
                                                                                                       tAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTcc                                              >  4:302057/1‑149 (MQ=255)
                                                                                                       tAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACACGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTcc                                              <  2:395503/149‑1 (MQ=255)
                                                                                                            acaATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTcc                                              <  4:736733/144‑1 (MQ=255)
                                                                                                            acaATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTcc                                              >  3:736733/1‑144 (MQ=255)
                                                                                                            acaATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATc                                         <  2:317101/149‑1 (MQ=255)
                                                                                                              aaTTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGa                                       >  2:572079/1‑149 (MQ=255)
                                                                                                              aaTTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGa                                       >  4:164600/1‑149 (MQ=255)
                                                                                                              aaTTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGa                                       <  4:584288/149‑1 (MQ=255)
                                                                                                                  ttCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACaaa                                   <  3:236921/149‑1 (MQ=255)
                                                                                                                       gTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGc                              <  1:674701/149‑1 (MQ=255)
                                                                                                                       gTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGc                              <  2:90611/149‑1 (MQ=255)
                                                                                                                        ttATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGcc                             >  2:389922/1‑149 (MQ=255)
                                                                                                                           taGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCttt                          >  2:79283/1‑149 (MQ=255)
                                                                                                                             ggAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTtgg                        >  3:721899/1‑149 (MQ=255)
                                                                                                                              gAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACa                                 >  4:386599/1‑139 (MQ=255)
                                                                                                                              gAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACa                                 <  3:386599/139‑1 (MQ=255)
                                                                                                                              gAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTtggt                       >  2:357118/1‑149 (MQ=255)
                                                                                                                                 aaCGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTtggtgg                     >  2:138949/1‑148 (MQ=255)
                                                                                                                                 aaCGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTtggtgg                     <  1:138949/148‑1 (MQ=255)
                                                                                                                                   cGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTAGGTGGGat                  <  4:67284/149‑1 (MQ=255)
                                                                                                                                     cTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGatat                <  4:721899/149‑1 (MQ=255)
                                                                                                                                        tGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGc             <  4:260686/149‑1 (MQ=255)
                                                                                                                                        tGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGc             >  3:285873/1‑149 (MQ=255)
                                                                                                                                         ggCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcc            <  1:474923/149‑1 (MQ=255)
                                                                                                                                              tttCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGc       >  4:351908/1‑149 (MQ=255)
                                                                                                                                                 cGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCcac    <  3:113864/149‑1 (MQ=255)
                                                                                                                                                  gATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCcaca   >  3:477147/1‑149 (MQ=255)
                                                                                                                                                   atcatcTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCCACAt  <  3:402572/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCCACAT  >  NC_000913/2238670‑2238965

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