Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,919,474 T→C 5.3% Y655Y (TAT→TAC)  yebT → MCE domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,919,4740TC 5.3% 236.4 / 2.6 75Y655Y (TAT→TAC) yebTMCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base C (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTC  >  NC_000913/1919330‑1919620
                                                                                                                                                |                                                                                                                                                  
ccTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCg                                                                                                                                                >  5:286574/1‑149 (MQ=255)
 cTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCgg                                                                                                                                               >  6:361461/1‑149 (MQ=255)
  tCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCgg                                                                                                                                               >  5:165418/1‑148 (MQ=255)
    aGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTa                                                                                                                                             <  6:233257/148‑1 (MQ=255)
      cGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTg                                                                                                                                          <  5:39854/149‑1 (MQ=255)
         tGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGata                                                                                                                                       >  8:202548/1‑149 (MQ=255)
            cAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCg                                                                                                                                    >  2:315929/1‑149 (MQ=255)
              gcgcCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCggg                                                                                                                                  <  4:290948/149‑1 (MQ=255)
                     tCAGCGTAAAGGCGATAAGAGTATTCTGTATGCTTCCGAAACTGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATc                                                                                                                            <  4:5818/148‑1 (MQ=255)
                          gTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACg                                                                                                                      >  1:105756/1‑149 (MQ=255)
                           tAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGc                                                                                                                     >  4:182780/1‑149 (MQ=255)
                             aaGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTg                                                                                                                   >  5:185882/1‑149 (MQ=255)
                             aaGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTg                                                                                                                   >  5:179359/1‑149 (MQ=255)
                               ggCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGa                                                                                                                 >  5:449761/1‑149 (MQ=255)
                                     taGCGTATTCTGTATGCTTAGGTAACATCGGCCCGTGCGGTTTTCTGGCTGATTACGCTTCACGGTTTCGATGCCGGAAAAGTGGCGGTCGTGATGCCATTTCGCTATCTCGGTATTGATATCGTGCAAATCCAGACGCTGGATGTGa                                                                                                            <  6:286574/147‑1 (MQ=255)
                                      aGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGAtt                                                                                                          <  7:390180/149‑1 (MQ=255)
                                           aTTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcg                                                                                                      >  3:88295/1‑148 (MQ=255)
                                            ttCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcgca                                                                                                    >  2:112979/1‑149 (MQ=255)
                                            ttCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcgca                                                                                                    >  1:83725/1‑149 (MQ=255)
                                            ttCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcgc                                                                                                     <  4:88295/148‑1 (MQ=255)
                                                     cTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGaa                                                                                           <  5:361461/149‑1 (MQ=255)
                                                      ttCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGGTGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTTTCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAg                                                                                          >  6:376746/1‑149 (MQ=255)
                                                      ttCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTGGATATCTCGGTTTTGATATCGGGCAAATCCAGCGGAGGGATATGATTACCGCACGCCATGAAg                                                                                          >  2:127265/1‑149 (MQ=255)
                                                      ttCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATg                                                                                             >  3:22904/1‑146 (MQ=255)
                                                      ttCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAg                                                                                          >  4:23740/1‑149 (MQ=255)
                                                        ccGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCTCGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAg                                                                                          <  4:234625/147‑1 (MQ=255)
                                                        ccGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGaa                                                                                           <  4:22904/146‑1 (MQ=255)
                                                        ccGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGt                                                                                         <  7:202548/148‑1 (MQ=255)
                                                        ccGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTGCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGAGGCCAATTCGCTATCTCGGTATTGATATCGGGCAAAGCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTa                                                                                        >  6:127947/1‑149 (MQ=255)
                                                         cGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTa                                                                                        >  8:303039/1‑148 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTa                                                                                                         <  8:278469/130‑1 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTa                                                                                                         >  7:278469/1‑130 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcgcacgc                                                                                                 <  7:400547/138‑1 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACcgcacgc                                                                                                 >  8:400547/1‑138 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACa                                                                                      <  3:175531/149‑1 (MQ=255)
                                                          gAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGATATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACa                                                                                      <  8:424875/149‑1 (MQ=255)
                                                           aaaCAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAg                                                                                     <  8:288636/149‑1 (MQ=255)
                                                           aaaCAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAg                                                                                     <  7:303039/149‑1 (MQ=255)
                                                             aCAGCGCCCCGTGCGGATGGCGGGCAGATAACGCTTAACGCTTTCGATGCCGGAAAACTGGCGCTAGGGATGCCAAGTCGCTATCTCGGTATTAATAACCCGCAAATCAACCCGCTGGATCTGAAAACCGCCCGCAAGTAAGTAAAgg                                                                                    >  6:56594/1‑148 (MQ=255)
                                                                     ccGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCgg                                                                           <  6:449761/149‑1 (MQ=255)
                                                                       gTGCGGTTTGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTa                                                                         <  6:179359/149‑1 (MQ=255)
                                                                       gTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTa                                                                         <  6:185882/149‑1 (MQ=255)
                                                                       gTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTa                                                                         <  3:113426/149‑1 (MQ=255)
                                                                            gTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTa                                                                    <  3:23740/149‑1 (MQ=255)
                                                                                gCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTAccc                                                                 <  2:23620/148‑1 (MQ=255)
                                                                                gCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCg                                                                <  8:312066/149‑1 (MQ=255)
                                                                                gCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCg                                                                <  3:182780/149‑1 (MQ=255)
                                                                                  gggCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGa                                                              <  2:105756/149‑1 (MQ=255)
                                                                                    gCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGaa                                                             >  1:365137/1‑148 (MQ=255)
                                                                                         ttACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTc                                                       >  1:322578/1‑149 (MQ=255)
                                                                                         ttACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTc                                                       <  7:267083/149‑1 (MQ=255)
                                                                                             gCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTc                                                       <  6:345499/145‑1 (MQ=255)
                                                                                             gCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTc                                                       >  5:345499/1‑145 (MQ=255)
                                                                                                   cGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTg                                             <  1:112979/149‑1 (MQ=255)
                                                                                                          gATGCCGGAATACTGGCGGTCGGGAGGCCAACTCGCAATCGCGGTATTGATATCGGGCAAAGCCAGACGCGGGATGTGATTACCGCACGCAAGTAAGTACAGGCAACGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCgg                                      <  5:56594/149‑1 (MQ=255)
                                                                                                           aTGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggt                                     <  2:83725/149‑1 (MQ=255)
                                                                                                            tGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtg                                    <  2:365137/149‑1 (MQ=255)
                                                                                                               cGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCgg                                                                           <  3:193341/107‑1 (MQ=255)
                                                                                                               cGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCgg                                                                           >  4:193341/1‑107 (MQ=255)
                                                                                                                 gAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTAcg                               >  5:377265/1‑149 (MQ=255)
                                                                                                                    aaCTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTCCCGCACGCAATGAAGTACAGGCACAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGGGGTACgcgc                            <  1:315929/149‑1 (MQ=255)
                                                                                                                         gCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTtct                        <  5:376746/148‑1 (MQ=255)
                                                                                                                             tCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTg                   <  6:165418/149‑1 (MQ=255)
                                                                                                                               gggATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGt                 >  7:386167/1‑149 (MQ=255)
                                                                                                                                 gATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTca               >  4:324012/1‑149 (MQ=255)
                                                                                                                                      cAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCg          >  7:35280/1‑149 (MQ=255)
                                                                                                                                       aaTTCGCTACCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAat                                                            >  6:213064/1‑98 (MQ=255)
                                                                                                                                       aaTTCGCTACCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAat                                                            <  5:206295/98‑1 (MQ=255)
                                                                                                                                       aaTTCGCTACCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAat                                                            <  5:213064/98‑1 (MQ=255)
                                                                                                                                       aaTTCGCTACCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAat                                                            >  6:206295/1‑98 (MQ=255)
                                                                                                                                             cTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt   <  2:322578/149‑1 (MQ=255)
                                                                                                                                              tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtgg                                   <  6:66828/116‑1 (MQ=255)
                                                                                                                                              tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtgg                                   >  5:66828/1‑116 (MQ=255)
                                                                                                                                              tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTc  >  7:382403/1‑149 (MQ=255)
                                                                                                                                              tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTc  >  1:470602/1‑149 (MQ=255)
                                                                                                                                                tctcGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTc  <  5:127947/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCTCAGCGGTGCCAGCGCCAGTCAGCGTAAAGGCGATAAGCGTATTCTGTATGCTTCCGAAACAGCGGCCCGTGCGGTTGGCGGGCAGATTACGCTTCACGCTTTCGATGCCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTC  >  NC_000913/1919330‑1919620

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