Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 520,603 T→C 5.8% V63A (GTG→GCG)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913520,6030TC 5.8% 206.6 / 3.0 69V63A (GTG→GCG) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base C (2/2);  total (33/36)
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.40e-01

CTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCT  >  NC_000913/520458‑520750
                                                                                                                                                 |                                                                                                                                                   
cTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGcc                                                                                                                                                  <  6:235679/149‑1 (MQ=255)
 tATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCg                                                                                                                                                 >  8:159498/1‑149 (MQ=255)
   ttAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGc                                                                                                                                                >  7:120756/1‑148 (MQ=255)
      aTTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGc                                                                                                                                            >  7:230021/1‑149 (MQ=255)
       ttGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCg                                                                                                                                           >  1:331918/1‑149 (MQ=255)
         gTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTg                                                                                                                                         >  7:163276/1‑149 (MQ=255)
          tCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTgg                                                                                                                                        <  5:273823/149‑1 (MQ=255)
                tggcGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCACTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGCGCCGCAAGCGTGGCTggag                                                                                                                                  >  6:234014/1‑149 (MQ=255)
                  gcgcTAAGCCTGGGGGTGGCCTGCGTGCTGGCGCTGGGCAATATGAGCGATCGCATGGAGAAGGGGTTAAGCCAGCAAAGACGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGTAGTGCCGCAAGCGTGGCTggag                                                                                                                                  <  8:277656/147‑1 (MQ=255)
                    gcTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAg                                                                                                                              >  3:226995/1‑149 (MQ=255)
                    gcTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCGTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAg                                                                                                                              <  7:248030/149‑1 (MQ=255)
                      tAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAgcg                                                                                                                            >  2:199700/1‑149 (MQ=255)
                      tAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAgcg                                                                                                                            >  8:191025/1‑149 (MQ=255)
                        aGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAgcg                                                                                                                            <  6:75791/147‑1 (MQ=255)
                           cTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCaaaa                                                                                                                       <  3:193160/149‑1 (MQ=255)
                                ggTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGcgcg                                                                                                                  <  3:1849/149‑1 (MQ=255)
                                      cTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTg                                                                                                             >  6:131175/1‑148 (MQ=255)
                                      cTGCGAGCTGGCGATGGGCAATATACGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGACGTTAGTTTATGGCGGGAGATCGGTCGTTGCGCACTTCACGCGAATTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGTCCTGa                                                                                                            <  3:204255/149‑1 (MQ=255)
                                         cGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGaaa                                                                                                          <  7:159498/148‑1 (MQ=255)
                                         cGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGaa                                                                                                           <  2:27395/147‑1 (MQ=255)
                                              tGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAACGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCgg                                                                                                     >  3:234355/1‑148 (MQ=255)
                                                gcgcTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGc                                                                                                    >  7:113700/1‑147 (MQ=255)
                                                gcgcTGGGCAATATCAGCGAGCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGc                                                                                                    <  3:180466/147‑1 (MQ=255)
                                                  gcTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTgga                                                                                                                                   >  4:3717/1‑114 (MQ=255)
                                                        cAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGa                                                                                          >  7:278918/1‑149 (MQ=255)
                                                        cAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGCGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGa                                                                                          >  5:105067/1‑149 (MQ=255)
                                                        cAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCACAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGa                                                                                          <  2:139722/149‑1 (MQ=255)
                                                         aaTATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGAc                                                                                         <  4:236630/149‑1 (MQ=255)
                                                           tatCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACtt                                                                                       >  5:125093/1‑149 (MQ=255)
                                                           tatCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACtt                                                                                       >  8:28284/1‑149 (MQ=255)
                                                             tCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGCTCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTc                                                                                     >  4:206655/1‑149 (MQ=255)
                                                                gCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCt                                                                                                                                      <  8:4251/97‑1 (MQ=255)
                                                                   aTCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCCCGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCcac                                                                                >  5:255516/1‑148 (MQ=255)
                                                                   aTCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCcaca                                                                               <  8:163276/149‑1 (MQ=255)
                                                                      gCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATg                                                                            >  3:43321/1‑149 (MQ=255)
                                                                       cATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCca                                                                                 <  5:223902/143‑1 (MQ=255)
                                                                       cATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCca                                                                                 >  6:223902/1‑143 (MQ=255)
                                                                        aTGGAGAAGGGCTTAAGCCAGCAAAGCCGGGAGTTAATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCGGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGAc                                                                          <  4:212616/149‑1 (MQ=255)
                                                                           gagaAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCt                                                                        <  8:113700/148‑1 (MQ=255)
                                                                             gaAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTg                                                                     >  2:46801/1‑149 (MQ=255)
                                                                                gggCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCat                                                                  >  6:331312/1‑148 (MQ=255)
                                                                                 ggCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAAt                                                                             >  4:69513/1‑137 (MQ=255)
                                                                                 ggCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAAt                                                                             <  3:69513/137‑1 (MQ=255)
                                                                                    ttAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCa                                                                   >  4:82622/1‑144 (MQ=255)
                                                                                    ttAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCa                                                                   <  3:82622/144‑1 (MQ=255)
                                                                                      aaGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGCGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGa                                                              <  6:105067/147‑1 (MQ=255)
                                                                                        gCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGacaca                                                          <  4:237278/149‑1 (MQ=255)
                                                                                               aaGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAGTGACCTTTGCAGGCGACACACCGCAGc                                                   >  4:130543/1‑149 (MQ=255)
                                                                                                         tttATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGt                                         <  1:244431/149‑1 (MQ=255)
                                                                                                             tGGCGGTCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGATGACTTTCGCCACAATTACCTTTGCAGGCGCCACACCGCAGCTg                                                 >  8:314306/1‑137 (MQ=255)
                                                                                                             tGGCGGGCGATCGGGCGTTGCGCATTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTTAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTg                                                 <  7:314306/137‑1 (MQ=255)
                                                                                                                  ggcgATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGt                                <  3:277688/149‑1 (MQ=255)
                                                                                                                  ggcgATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGt                                <  5:211070/149‑1 (MQ=255)
                                                                                                                  ggcgATCGGGCGTTGCGCAGTTCACGCGAAGCGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCggg                                <  5:234014/149‑2 (MQ=255)
                                                                                                                     gATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGga                             <  4:226995/149‑1 (MQ=255)
                                                                                                                      aTCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTATCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgat                            <  6:255516/149‑1 (MQ=255)
                                                                                                                         gggCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgatga                          >  1:148829/1‑148 (MQ=255)
                                                                                                                         gggCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGgatga                          <  6:125093/148‑1 (MQ=255)
                                                                                                                           gCGTTGCGCAGTTAACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGTCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGatat                       >  4:50171/1‑149 (MQ=255)
                                                                                                                            cGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATc                      <  3:40883/149‑1 (MQ=255)
                                                                                                                                 cgcAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTAcc                  <  7:28284/148‑1 (MQ=255)
                                                                                                                                   cAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGa               >  2:9434/1‑149 (MQ=255)
                                                                                                                                    aGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGAt              <  3:206655/149‑1 (MQ=255)
                                                                                                                                    aGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGAt              <  7:191025/149‑1 (MQ=255)
                                                                                                                                    aGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGAt              >  3:257685/1‑149 (MQ=255)
                                                                                                                                         aCGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATg         >  2:135384/1‑149 (MQ=255)
                                                                                                                                          cgcgAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATgg        <  2:148829/149‑1 (MQ=255)
                                                                                                                                           gcgAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGc       <  1:46801/149‑1 (MQ=255)
                                                                                                                                                gTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCt  >  8:286236/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTATTAATTGTCTGGCTGGCGCTAAGCCTGGCGGTGGCCTGCGTGCTGGCGCTGGGCAATATCAGCGATCGCATGGAGAAGGGCTTAAGCCAGCAAAGCCGTGAGTTTATGGCGGGCGATCGGGCGTTGCGCAGTTCACGCGAAGTGCCGCAAGCGTGGCTGGAGGAAGCGCAAAAGCGCGGCCTGAAAGTCGGCAAGCAGCTGACTTTCGCCACAATGACCTTTGCAGGCGACACACCGCAGCTGGCGAACGTCAAAGCGGTGGATGATATCTACCCGATGTATGGCGATCT  >  NC_000913/520458‑520750

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