Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,865,563 C→G 5.7% L237L (CTG→CTC)  yidE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,865,5630CG 5.7% 275.4 / 6.1 82L237L (CTG→CTC) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base G (2/3);  total (38/44)
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.08e-01

TGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCT  >  NC_000913/3865415‑3865710
                                                                                                                                                    |                                                                                                                                                   
tGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTc                                                                                                                                                     <  1:244353/149‑1 (MQ=255)
tGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTc                                                                                                                                                     <  6:50866/149‑1 (MQ=255)
tGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTc                                                                                                                                                     <  2:102516/149‑1 (MQ=255)
    tttCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAc                                                                                                                                                 >  4:296673/1‑149 (MQ=255)
    tttCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAc                                                                                                                                                 >  1:54125/1‑149 (MQ=255)
      tCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAc                                                                                                                                                 <  8:37493/147‑1 (MQ=255)
       cGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgc                                                                                                                                               <  7:179088/148‑1 (MQ=255)
        gTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgcg                                                                                                                                              <  4:167376/148‑1 (MQ=255)
         tGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTCTGAAACGCGaa                                                                                                                                            <  5:225478/149‑2 (MQ=255)
          ggACAGCGAAGTATCGACCTCCTGACCACTCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGc                                                                                                                                           <  5:145134/149‑1 (MQ=255)
           gACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCTCATTATGTAAATCCACTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATTGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCa                                                                                                                                          <  1:105722/149‑1 (MQ=255)
                    gTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAAtttt                                                                                                                                 <  1:253834/149‑1 (MQ=255)
                     tATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTg                                                                                                                                >  3:315155/1‑149 (MQ=255)
                       tCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCa                                                                                                                                          <  6:167155/137‑1 (MQ=255)
                       tCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCa                                                                                                                                          >  5:167155/1‑137 (MQ=255)
                          aCCTCCTGTCCAATCACCAGAGTCGCATTATGTAAATCCGCTGGCTGACCCACCAGCTGCAGCATATCGCCCAGTTGGATAATTGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGcc                                                                                                                           <  8:174319/149‑1 (MQ=255)
                           cctcctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACAAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGTAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCg                                                                                                                          >  6:186724/1‑149 (MQ=255)
                              cctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCCGGCGAAGGAAC‑TTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTga                                                                                                                      <  6:281229/149‑1 (MQ=255)
                               ctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCTTTg                                                                                                                       <  2:179114/148‑1 (MQ=255)
                               ctGACCAATCAACAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTCGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTga                                                                                                                      <  8:83495/149‑1 (MQ=255)
                                tGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGtt                                                                                                                        >  3:113324/1‑146 (MQ=255)
                                      aTCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCgg                                                                                                               <  5:83144/149‑1 (MQ=255)
                                         acgaGTTGCGCATTATGTAAATCCGCTGGTTGAACCACCAGATGCATCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTa                                                                                                             <  4:98048/145‑1 (MQ=255)
                                                gcgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                      <  2:54125/148‑1 (MQ=255)
                                                gcgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTGGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                      >  8:122945/1‑148 (MQ=255)
                                                 cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTaa                                                                                                    <  6:248805/149‑1 (MQ=255)
                                                 cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTa                                                                                                     <  8:63377/148‑1 (MQ=255)
                                                   cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATg                                                                                                  <  7:323069/149‑1 (MQ=255)
                                                   cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGAGAATTTTGTCGCCGGTGAGAATCGGTACATCTTTAATg                                                                                                  >  1:282068/1‑149 (MQ=255)
                                                       aTGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCGGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCt                                                                                               >  1:285320/1‑148 (MQ=255)
                                                       aTGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCt                                                                                               >  8:284154/1‑148 (MQ=255)
                                                          tAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAAt                                                                                           <  8:229563/149‑1 (MQ=255)
                                                          tAAATCCGCTGGCGGAAGCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGATCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTatat                                                                                           <  5:264908/149‑1 (MQ=255)
                                                            aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCa                                                                                         <  4:61395/149‑1 (MQ=255)
                                                              tCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATg                                                                                       >  4:135966/1‑149 (MQ=255)
                                                                   tggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGAAGGt                                                                                  >  2:109748/1‑149 (MQ=255)
                                                                     gctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGt                                                                                  <  7:254568/147‑1 (MQ=255)
                                                                      ctgACCCACCAGATGCAGCAAATCGCCAAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACTCGATCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAg                                                                               <  1:307318/149‑1 (MQ=255)
                                                                            caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                         >  8:190941/1‑149 (MQ=255)
                                                                             accaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCTGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGCACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                        <  3:242119/149‑1 (MQ=255)
                                                                             accaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTAGACTTTTAGAATCGTTACATCTTTAATTGCTAAATAATGCAGGTTAGGGttat                                                                         >  8:334206/1‑146 (MQ=255)
                                                                              ccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                        >  7:225131/1‑148 (MQ=255)
                                                                                 gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGTTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCtt                                                                                                       >  5:236245/1‑114 (MQ=255)
                                                                                 gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCtt                                                                                                       <  6:236245/114‑1 (MQ=255)
                                                                                 gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcg                                                                    >  7:330908/1‑149 (MQ=255)
                                                                                  aTGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgc                                                                   <  8:225131/149‑1 (MQ=255)
                                                                                  aTGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAAATTTTAGGGTTTCTTCGCTGTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgc                                                                   >  6:157247/1‑149 (MQ=255)
                                                                                     cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                      >  1:1167/1‑111 (MQ=255)
                                                                                     cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                      <  2:1167/111‑1 (MQ=255)
                                                                                     cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAgg                                                                              >  2:181647/1‑135 (MQ=255)
                                                                                     cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAgg                                                                              <  1:181647/135‑1 (MQ=255)
                                                                                      agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGaa                                                                >  8:243205/1‑148 (MQ=255)
                                                                                      agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGaa                                                                >  7:282274/1‑148 (MQ=255)
                                                                                        caAATCGCCCGGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAatat                                                             <  2:285320/149‑1 (MQ=255)
                                                                                           aTCGCCCAGTTGGATACTGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTGAGGGTGCTCAACGCGAATATTGa                                                          <  7:334206/149‑1 (MQ=255)
                                                                                                   gTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTCTGAAACGCGaa                                                                                                                                            <  7:31568/59‑2 (MQ=255)
                                                                                                   gTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTCTGAAACGCGaa                                                                                                                                            >  8:31568/1‑58 (MQ=255)
                                                                                                    ttGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGa                                                 >  4:234856/1‑149 (MQ=255)
                                                                                                    ttGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGa                                                 >  1:275901/1‑149 (MQ=255)
                                                                                                      ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGAt                                                >  6:134521/1‑148 (MQ=255)
                                                                                                      ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGAt                                                <  3:296673/148‑1 (MQ=255)
                                                                                                        aTAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGATTATTGATAGTCTTGATCAg                                             >  3:284749/1‑149 (MQ=255)
                                                                                                           aTGGTATCTGGCGAAGGAACTGTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcgc                                          >  6:303955/1‑149 (MQ=255)
                                                                                                                       gAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGt                              <  7:190941/149‑1 (MQ=255)
                                                                                                                         aGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGc                            <  8:282274/149‑1 (MQ=255)
                                                                                                                           gAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGt                          <  7:274538/149‑1 (MQ=255)
                                                                                                                             aCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTg                         >  2:301047/1‑148 (MQ=255)
                                                                                                                              cTTTTAGGGTTTCTTCGCTTATCAGACGCGAGCAGATAATTATGTCGCCGCTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGaa                       <  7:122945/149‑1 (MQ=255)
                                                                                                                               ttttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTg                         >  3:202281/1‑146 (MQ=255)
                                                                                                                                tttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGCTCAGCGCGCCGCCATTGGTGCTGGAAGa                     >  2:108781/1‑149 (MQ=255)
                                                                                                                                tttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                     >  3:155365/1‑149 (MQ=255)
                                                                                                                                 ttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCTGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGAc                    >  6:55634/1‑149 (MQ=255)
                                                                                                                                 ttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTACTGGCTAAATCCTGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                     <  1:237261/148‑1 (MQ=255)
                                                                                                                                 ttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                     <  7:284154/148‑1 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                     <  4:118124/111‑1 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                     >  3:118124/1‑111 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGa                                                 <  3:182062/115‑1 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGa                                                 >  4:182062/1‑115 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtg               <  3:234856/149‑1 (MQ=255)
                                                                                                                                      gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtg               >  5:338512/1‑149 (MQ=255)
                                                                                                                                       tttCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa                                                          <  4:13100/105‑1 (MQ=255)
                                                                                                                                       tttCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa                                                          >  3:13100/1‑105 (MQ=255)
                                                                                                                                           ctcGCGTCTGAAACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTgctgct          >  6:143539/2‑149 (MQ=255)
                                                                                                                                              gcgGTTGAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGGTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAg       <  4:315155/149‑1 (MQ=255)
                                                                                                                                                 tttCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt    <  1:109748/149‑1 (MQ=255)
                                                                                                                                                    cAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCt  <  4:202281/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCT  >  NC_000913/3865415‑3865710

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