Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,255 G→C 6.0% S196T (AGC→ACC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2550GC 6.0% 341.8 / 8.5 121S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (56/57);  new base C (3/5);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 121 total aligned reads displayed.

CGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGG  >  NC_000913/4262114‑4262395
                                                                                                                                             |                                                                                                                                            
cggtgactgtgaATACGCGTATTGGCATCGACGACCAGCACAGCTATGAATTTCTCTGCTATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGc                                                                                                                                             <  2:227902/143‑1 (MQ=255)
cggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTGTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGc                                                                                                                                             <  2:141028/143‑1 (MQ=255)
    gacggtgaCACCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                                                                         <  1:1007161/143‑1 (MQ=255)
    gacggtgaAAAGGCGTATCGGCATCGACGACCAGGACAGCGATGAAGTTCTCTGCGATTTCATCAATACCGTTGCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                                                                         <  2:671028/143‑1 (MQ=255)
    gacggtgaAAACGCGCATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                                                                         <  1:36665/143‑1 (MQ=255)
      cggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                                                                       >  1:819242/1‑143 (MQ=255)
      cggtgaAAAAGCGTACTGGCAGCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa                                                                                                                                       <  1:527497/143‑1 (MQ=255)
      ccgtgaAAACGCGTATTCGACTAAACGACCATTACAGCTATCAAGTGCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa                                                                                                                                       <  1:473603/141‑1 (MQ=255)
       ggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                                                                      >  2:1223299/1‑143 (MQ=255)
        gtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGGAGAAACCTGCGATGTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCGGGCTTCCGGGGTTAAGCACgaagga                                                                                                                                     >  2:700968/1‑143 (MQ=255)
         tgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaa                                                                                                                                    <  2:455002/143‑1 (MQ=255)
          gaAACCGCGGATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                                                                   <  1:1123380/143‑4 (MQ=255)
           aaaaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTTATCATCCACGCACGTAACGCCTGGCCTTCGGGGATAAGCCCGAAAGaaaa                                                                                                                                  >  2:874646/1‑143 (MQ=255)
           aaaaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGCGGTTAAGCCCGAAAGcaaa                                                                                                                                  >  2:1057543/1‑143 (MQ=255)
            agaCGCGTATTGGCAGCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAGTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccagg                                                                                                                                 <  2:450451/141‑6 (MQ=255)
               cgcgTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGCTAAGCCCGAAAGAAAACCGt                                                                                                                              >  2:1135505/1‑143 (MQ=255)
                 cggaTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCAAGCAAGTACAGTTCAATTTTCGTGGTTAAGCCCGAAAGAAAACCGTGa                                                                                                                            <  1:244171/140‑1 (MQ=255)
                  gtccgggCATCGGCCACCAGGACAGCGATGAATTTCTCTGCGATTTCCTCAATACCGCTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                                                                           <  2:638729/138‑1 (MQ=255)
                   tATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaaa                                                                                                                          <  1:1191186/143‑1 (MQ=255)
                    atgggCACCGCCGACCAGGCCAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAAt                                                                                                                         <  2:1090746/140‑1 (MQ=255)
                     ttGGCATCGACGACCAGGACAGCTATGCATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAACGCCTTGCTTCCGGCGTTCAGCCCGCCAGAAACCCGTCACATc                                                                                                                        >  1:276428/1‑143 (MQ=255)
                        gCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCAAGCACGTAAAGCCTGGCTTTCGGGGTTACGCCCGAAAGAAAACCGTGAAATCccg                                                                                                                     >  1:196583/1‑143 (MQ=255)
                         cATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGCAATCccgc                                                                                                                    >  2:709108/1‑143 (MQ=255)
                           tCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgccg                                                                                                                  <  1:137137/143‑1 (MQ=255)
                           tCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGGTAAGCCCGAAAGCACACAGCGTAACCccgccg                                                                                                                  >  2:195240/1‑143 (MQ=255)
                           tCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAATGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAACTCccgccg                                                                                                                  >  1:160708/1‑143 (MQ=255)
                              acgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAACGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGACCGCAAGAAAACCGTGAAATCCcgccgcgc                                                                                                               >  2:253261/1‑141 (MQ=255)
                                gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCGTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCAGCTCGa                                                                                                             >  1:1202988/1‑143 (MQ=255)
                                    aGGACGGCTGCGGATTTCTCTGCTGTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTTAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                                                         <  2:321165/143‑1 (MQ=255)
                                          gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgt                                                                                                   <  1:267436/143‑1 (MQ=255)
                                           cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGCTGTTCATCATCCACGCACGTAACGCCTGGCTTGCGGGGTGAAGGCCGCACGAAAGCCGTGAAATCCCGCGGCTCGATAATCcgcgcg                                                                                                  >  2:380796/1‑141 (MQ=255)
                                            tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                                                                                 <  2:521067/143‑1 (MQ=255)
                                            tATGAATTTCCCTGCGATTTCAGCAATTCCGTTCCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                                                                                 <  2:673830/143‑1 (MQ=255)
                                               gAATTTCTCTGCGATTTCATCAATCCCGTTTCCGGCAAAGCCGAGTGGGAGATGGTCATCATCCCCGCACGTAGCGCCTGGCTTTCGGGGTTAAGCGCGAAAGAAAACCGTGACATCACGCCGGCGGCTCATCCGCATGTGTa                                                                                              >  1:430909/1‑143 (MQ=255)
                                                aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTAt                                                                                             >  1:509478/1‑143 (MQ=255)
                                                aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCAGCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTAt                                                                                             <  1:1246411/143‑1 (MQ=255)
                                                 aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                                                                            >  1:1135601/1‑143 (MQ=255)
                                                 aTTTCTCGGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                                                                            <  1:103104/143‑1 (MQ=255)
                                                   ttCTCGGCGATTGCATCAATACCGTGCCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCaa                                                                                          <  1:560440/143‑1 (MQ=255)
                                                      tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGCGTGTATCAACTg                                                                                       >  1:407270/1‑143 (MQ=255)
                                                       ttGCGGTTTCATCAATACCGTTTTCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGa                                                                                      <  1:1014000/142‑1 (MQ=255)
                                                           gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTCAGCCCGAAACAAAACCGTGAAATCCCGCCGCTCGATTATACGCGTGTGTATCAACTGAAGCg                                                                                  >  1:95798/1‑143 (MQ=255)
                                                           gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTACGCCCGCAAGAAAACCGTGACATCCCGCCGCCCGATTCTCCGCGTGTGTctcaactcaaccc                                                                                  >  2:1210536/1‑140 (MQ=255)
                                                              ttCGTCAATACCGTTTCCGGCAAAGGCGAGGGGGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                               <  2:502436/143‑1 (MQ=255)
                                                              ttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                               <  1:765870/143‑1 (MQ=255)
                                                              ttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                               >  1:1047427/1‑143 (MQ=255)
                                                                 atcaggaCCGGTTCCGGCAAAGGCGAGTGTGAGATGTTCGTCATCCACGCACGTAAAGCCTGGCGTTCGGGGTTAAGCCGGAAAGAAAACCGTGAAAGCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt                                                                            <  2:540178/137‑1 (MQ=255)
                                                                 atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGGGTGTATCAACTGAAGCGTGACtt                                                                            >  1:656123/1‑143 (MQ=255)
                                                                  tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGCAAACCGTGAAATCCCGCCGCTCGATTCTCCGCGCGTGCATCAACTGCAGCGTGCCtct                                                                           >  2:340163/1‑141 (MQ=255)
                                                                   caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                                                          <  2:1124512/143‑1 (MQ=255)
                                                                     aTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg                                                                        >  1:1161225/1‑143 (MQ=255)
                                                                     aTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTACCAACTAAAGCGTGCCTTTCCg                                                                        >  2:47779/1‑143 (MQ=255)
                                                                      tACCGTTTCTGGCAAAGGCGAGTGTGAGGTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCGGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGc                                                                       <  2:89988/143‑1 (MQ=255)
                                                                       aCCGTTTCCTGCAAAGGCGAGTGTCAGATTTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                                                                      <  1:148869/143‑1 (MQ=255)
                                                                       aCCGTTTCCGGCAAAGGCGAGTATCAGATCATTCTCATTTACTGACGTAAAGCCTGGCTTTCGGGGGTAAGCACTGGATATGACGGTTAATGTCCACCGCGCGAACATCCGCGTTTATAGCAACAGACGCGTGACTCTCCGCa                                                                      >  2:972307/1‑143 (MQ=255)
                                                                        ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTATCGTGGTTAAGCCCGCAAGAAACCCGCGCAATCCCCCCGCCCGACTATCCGCCTGCGCACCCACTGACGCGTGACCTTCCGCac                                                                     >  2:843449/1‑142 (MQ=255)
                                                                          gttgcgtgCAAAGGCGAGTGTGAGATGGTCATCATCCACGCACGTAAAGCATGGGTTTAGCGACGAGCCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                                                                   <  1:1271570/136‑1 (MQ=255)
                                                                          gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTGCGGATCt                                                                   >  2:1166298/1‑143 (MQ=255)
                                                                            ttgcGGCAAACGCGAGGGTGACATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGa                                                                 <  1:829628/140‑1 (MQ=255)
                                                                            ttCCGGCAAAGGCGAGTGTGAGATGTTAATCATCCACGCGCGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTgc                                                                 >  2:848423/1‑142 (MQ=255)
                                                                              ccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                                                               >  2:591383/1‑143 (MQ=255)
                                                                                  cgacggCGAGTGTGAGATGTTCACCCGCCACGCACGTAAACGCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt                                                           <  2:619794/139‑1 (MQ=255)
                                                                                  cAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGGGACTATCCGCACCTGACacagt                                                           >  2:1126013/1‑139 (MQ=255)
                                                                                       gCGAGTATGAGGTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAtt                                                      <  1:1101753/143‑1 (MQ=255)
                                                                                               gggaTGTTCATCATCCGCGCACGTGACGCCTGCACTCCTACGGGAAGGCGGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg                                              <  2:563183/141‑1 (MQ=255)
                                                                                               gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg                                              <  1:857862/143‑1 (MQ=255)
                                                                                                agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtggt                                             >  1:146469/1‑143 (MQ=255)
                                                                                                      tcatcaTCCACGCACGGAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAg                                       <  1:974204/143‑1 (MQ=255)
                                                                                                        ctcatcCACGTACGTAACGCAGGGCTTTCGGGGTTAAGCCGAAAATAAAACCGGGAAATCCAGCCGCTCGATCATCCGCGTGTGTATCAACTGAAGCGTGACGTCCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTc                                     <  1:253261/142‑1 (MQ=255)
                                                                                                           atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                  <  2:1084878/143‑1 (MQ=255)
                                                                                                           atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                  <  1:1281592/143‑1 (MQ=255)
                                                                                                               gcgcacgTAAACCCTGGCTTTCGGGGTTAAGCCCTAAAGAAAACCGTGAACTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTCGTATCAAGTCGCTGCaa                              <  2:1050509/142‑1 (MQ=255)
                                                                                                               acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACATTAAGATATTAAGAAAGTCGATTAACGGTGGTATCAAATCGCTGGaa                              >  2:1239786/1‑143 (MQ=255)
                                                                                                               acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCCGGTGTATCAACTGAAGCGTCGCTTTACGCGTCTGACACTGTCGATTAGCGGTGGTCGCATGTCGCtccga                              >  2:1226809/1‑139 (MQ=255)
                                                                                                                 gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAACGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAGTAACGGTGGTATCAAGACGCTGGaaca                            >  1:1031752/1‑141 (MQ=255)
                                                                                                                  catgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                           <  1:983746/140‑1 (MQ=255)
                                                                                                                  cacggcaaGCCTGGCTTTCGGGGTTACGCCCGAAAGCAACCCGTGAAATCCCGCCGCTCGATTCTCCGCGTGTGCATCAACTGAAGCGTGACTTTCCGCTTCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAagag                           <  1:700968/137‑1 (MQ=255)
                                                                                                                     gtagggCCGGGCTTTCGGGGTTAAGCCCGAAAGAGAACCGTGAAGTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGCGGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGcc                        <  1:257240/138‑1 (MQ=255)
                                                                                                                      tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCACGTCGCTGGCAGAGGCCa                       >  1:30765/1‑143 (MQ=255)
                                                                                                                         tgCCTCGCTTTCGGGCGTGAGCCCGCAAGAAAACCGTGAAAGCCCGCCGCTCGATTAGCCGCGTGTGGATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAGCAAGTCGCTGGAAGAGGCCAAAg                    <  2:1008086/142‑1 (MQ=255)
                                                                                                                          gCCTGGCTTTCGGGGTTAAGGCCGAAAGCAAACCGGGGAAGCCCGCCGGTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGc                   <  2:559748/143‑1 (MQ=255)
                                                                                                                           ccTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                  <  1:703485/143‑1 (MQ=255)
                                                                                                                             tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGCAGCGGCCAAAGcaca                >  1:316050/1‑143 (MQ=255)
                                                                                                                             tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCCAAGcaca                >  1:738146/1‑143 (MQ=255)
                                                                                                                             tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCACAGcaca                >  1:1192524/1‑143 (MQ=255)
                                                                                                                             tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCGCGATTATCCGCGTGTGTATCAACTGAAGCGTCACTTTCCGCATCTGACGATGTCGATTCACGGCGGTCTCAAGTCGCCGGAGGAGGCCAAAGcaca                >  2:307401/1‑143 (MQ=255)
                                                                                                                              ttcTTTCGGGCTTAACCTCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCAGCGTGTCTACCAACTGAAGCGTGACCTTCCGCATCTGACAATGTCGCTTAACGGCGATATCAAGTCGCTGGAAGAGGCCAAACCacac               >  2:842858/3‑143 (MQ=255)
                                                                                                                              ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCacac               >  1:72194/1‑143 (MQ=255)
                                                                                                                                ctcgcGGGGTTCAGCCGGAAAGCAAACGGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTGCCGCATCTGACAAGGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCt             <  2:185064/139‑1 (MQ=255)
                                                                                                                                cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTGCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCt             <  2:360251/143‑1 (MQ=255)
                                                                                                                                cTTTCGGGGTTAAGCCAGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGCAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCt             <  1:942263/143‑1 (MQ=255)
                                                                                                                                 tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTg            >  1:591752/1‑143 (MQ=255)
                                                                                                                                 tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCGGGAAGAGGCCAAAGCACACCTg            <  1:294724/143‑1 (MQ=255)
                                                                                                                                  ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCATGACTTTCCGCATATGACAATGTCGATTAACGGGGGTATCGAGTCGCTGGGAGAGGCCAAAGGACACCTGc           >  2:539344/1‑143 (MQ=255)
                                                                                                                                  ttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc           >  2:781107/1‑143 (MQ=255)
                                                                                                                                    cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa         <  1:1135505/143‑1 (MQ=255)
                                                                                                                                     ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGTAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCACAGCACACCTGCAAc        >  1:423389/1‑143 (MQ=255)
                                                                                                                                      gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACCTTCCGCATCTGACAATGTCGATTAACGGTGGGATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCGACa       >  1:518117/1‑143 (MQ=255)
                                                                                                                                        gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACata     <  1:754315/143‑1 (MQ=255)
                                                                                                                                           aaGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACCCCTGCAACACatgg  >  1:62480/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGG  >  NC_000913/4262114‑4262395

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