Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,408,758 T→C 5.7% L702P (CTG→CCG) ‡ rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,408,7580TC 5.7% 405.4 / 8.7 122L702P (CTG→CCG) ‡rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base T (66/49);  new base C (2/5);  total (68/54)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 127 total aligned reads displayed.

AGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAGGC  >  NC_000913/4408610‑4408906
                                                                                                                                                    |                                                                                                                                                    
aGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCt                                                                                                                                                      <  2:798798/149‑1 (MQ=255)
  cGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTgg                                                                                                                                                    <  1:377484/149‑1 (MQ=255)
    tCACTGGCTTTGGCTTCTTCGTCCGTCTGGAGGACTTGTTCATGGATGGTCTGGTCCGTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGTTAGGGCATCGCCTGAGGCGGGAATCGAGCTGCCAGACTTATCGCCTGGGc                                                                                                                                                  >  4:82046/1‑149 (MQ=255)
         ggctttggcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaa                                                                                                                                             <  1:122467/149‑3 (MQ=255)
            tttggcttCTTCGTCCGTCTGGACGACTTGTTCGTTGATGGTCTGGTCCATGTCTCTTCGCTGGACCATGACTACTATCGCTTTGACCAGGTCGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGt                                                                                                                                          >  1:843944/1‑149 (MQ=255)
            tttggcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaagtc                                                                                                                                          <  3:248149/149‑6 (MQ=255)
             ttggcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaagtct                                                                                                                                         <  3:598851/149‑7 (MQ=255)
                gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                                                                                                                      >  4:795662/1‑149 (MQ=255)
                  ttcttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGt                                                                                                                                    >  1:345398/1‑149 (MQ=255)
                  ttcttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCAGCCAGACTTATCGCCCAGGCGAtaa                                                                                                                                             <  2:820685/140‑3 (MQ=255)
                  ttcttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCCGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaa                                                                                                                                             >  1:820685/1‑138 (MQ=255)
                      tcGTCCGTCTGGACGACTTGTTCATTGATGGTCTAGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                                                                                                                                <  3:528286/149‑1 (MQ=255)
                         tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCAGGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTc                                                                                                                             >  1:734768/1‑149 (MQ=255)
                          ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaa                                                                                                                                             <  3:636330/132‑3 (MQ=255)
                          ccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCCAGGCGAtaa                                                                                                                                             >  4:636330/1‑130 (MQ=255)
                             tcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                                                                                                                         >  4:536924/1‑149 (MQ=255)
                              cTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTTAGGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                                                                                                                        <  1:748476/149‑1 (MQ=255)
                              cTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                                                                                                                        <  1:519057/149‑1 (MQ=255)
                                    gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                                                                                                                  >  3:345100/1‑149 (MQ=255)
                                    gacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                                                                                                                  >  3:647603/1‑149 (MQ=255)
                                       ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                                                                                                                               <  3:157526/117‑1 (MQ=255)
                                       ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat                                                                                                               >  4:53661/1‑149 (MQ=255)
                                         gTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATgg                                                                                                             <  2:764018/149‑1 (MQ=255)
                                             aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa                                                                                                         <  1:405419/149‑1 (MQ=255)
                                             aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa                                                                                                         >  1:497915/1‑149 (MQ=255)
                                             aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa                                                                                                         >  2:880382/1‑149 (MQ=255)
                                             aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa                                                                                                         >  2:405419/1‑149 (MQ=255)
                                              ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg                                                                                                        >  3:865841/1‑149 (MQ=255)
                                              ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAg                                                                                                        >  1:270946/1‑149 (MQ=255)
                                               tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc                                                                                                       <  4:255737/149‑1 (MQ=255)
                                               tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc                                                                                                       <  2:734768/149‑1 (MQ=255)
                                                   ggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCaa                                                                                                   >  2:469438/1‑149 (MQ=255)
                                                     tctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCaaaa                                                                                                 <  1:858323/149‑1 (MQ=255)
                                                      ctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAAt                                                                                                <  2:863699/149‑1 (MQ=255)
                                                       tggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATc                                                                                               <  1:898025/149‑1 (MQ=255)
                                                         gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGa                                                                                             >  2:164127/1‑149 (MQ=255)
                                                            cATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACtt                                                                                          <  2:270946/149‑1 (MQ=255)
                                                             aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTg                                                                                                                                         <  3:458155/101‑1 (MQ=255)
                                                             aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTg                                                                                                                                         >  4:458155/1‑101 (MQ=255)
                                                             aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACttt                                                                                         >  2:295402/1‑149 (MQ=255)
                                                              tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTa                                                                                        >  4:392907/1‑149 (MQ=255)
                                                               gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAg                                                                                       >  4:459421/1‑149 (MQ=255)
                                                                tctcTTCGCTGGACGATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGGTGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGc                                                                                      >  3:879889/1‑149 (MQ=255)
                                                                tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgcg                                                                                                      >  1:85153/1‑133 (MQ=255)
                                                                tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGc                                                                                      <  1:612891/149‑1 (MQ=255)
                                                                tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGc                                                                                      >  4:66765/1‑149 (MQ=255)
                                                                  tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCt                                                                                    >  1:69551/1‑149 (MQ=255)
                                                                    ttCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGa                                                                                  >  4:576287/1‑149 (MQ=255)
                                                                       gCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGAtct                                                                               >  1:354297/1‑149 (MQ=255)
                                                                       gCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGAtct                                                                               >  4:50370/1‑149 (MQ=255)
                                                                                tGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAAc                                                                      <  1:535020/149‑1 (MQ=255)
                                                                                 gACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAAcg                                                                     >  3:378756/1‑149 (MQ=255)
                                                                                  actactATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAAcgc                                                                    >  3:422477/1‑149 (MQ=255)
                                                                                  actactATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAAcgc                                                                    >  1:821898/1‑149 (MQ=255)
                                                                                   ctactaTCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGAc                                                                                            >  3:385933/1‑124 (MQ=255)
                                                                                    tactaTCGCTTTGACCAGGTGGGGCAACGCCTGATGGGGGAATCCAGCGACCAGACTTATCGCCTGGGCGATGGCGTGGAAGTTCGCGGCGAAGCGGATAATATGGACGAGCGCAAAGTCGACTTTAGCCTGATCTCCAGCGAACgcgc                                                                  >  2:330355/1‑149 (MQ=255)
                                                                                      ctaTCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCAc                                                                <  4:341143/149‑1 (MQ=255)
                                                                                      ctaTCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCAc                                                                >  4:114630/1‑149 (MQ=255)
                                                                                       taTCGCTTTGACTAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCTCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATTTCCAGCGAACCCGCAcc                                                               >  4:745339/1‑149 (MQ=255)
                                                                                         tCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACcgc                                                             >  2:843001/1‑149 (MQ=255)
                                                                                         tCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACcgc                                                             <  3:665377/149‑1 (MQ=255)
                                                                                          cGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACcgcg                                                            <  1:81780/149‑1 (MQ=255)
                                                                                           gCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCgcgc                                                           >  3:578164/1‑149 (MQ=255)
                                                                                            cTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCa                                                          >  2:725717/1‑149 (MQ=255)
                                                                                             tttGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCaa                                                         >  2:864531/1‑149 (MQ=255)
                                                                                             tttGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCaa                                                         >  1:89885/1‑149 (MQ=255)
                                                                                             tttGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGCCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCaa                                                         >  4:600283/1‑149 (MQ=255)
                                                                                                 aCCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTc                                                     >  1:89520/1‑149 (MQ=255)
                                                                                                   cAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCgg                                                   >  4:208330/1‑149 (MQ=255)
                                                                                                   cAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCgg                                                   <  1:553823/149‑1 (MQ=255)
                                                                                                    aGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGt                                                  >  3:856622/1‑149 (MQ=255)
                                                                                                    aGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGt                                                  <  2:150424/149‑1 (MQ=255)
                                                                                                           gCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACgg                                           >  1:855156/1‑149 (MQ=255)
                                                                                                             aaCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGgcg                                         >  1:224049/1‑149 (MQ=255)
                                                                                                                 ccTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGcgcgcg                                     >  2:559086/1‑149 (MQ=255)
                                                                                                                  cTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCga                                    >  4:359663/1‑149 (MQ=255)
                                                                                                                   tGATGGGGGAATCCAGCGGACAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCCCGTCGAAGCGGCTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCgag                                   >  4:610925/1‑149 (MQ=255)
                                                                                                                      tGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGaaa                                >  4:804562/1‑149 (MQ=255)
                                                                                                                      tGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGaaa                                >  3:566625/1‑149 (MQ=255)
                                                                                                                      tGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCCTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGaaa                                <  2:280552/149‑1 (MQ=255)
                                                                                                                         gggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCg                             >  2:12911/1‑149 (MQ=255)
                                                                                                                         gggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCg                             <  4:353011/149‑1 (MQ=255)
                                                                                                                         gggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCg                             >  3:364746/1‑149 (MQ=255)
                                                                                                                         gggAATCCAGCGCCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCg                             <  3:600283/149‑1 (MQ=255)
                                                                                                                          ggAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGa                            >  2:106556/1‑149 (MQ=255)
                                                                                                                               ccAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGaaaaaa                       <  4:647603/149‑1 (MQ=255)
                                                                                                                                cAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAg                      <  2:497915/149‑1 (MQ=255)
                                                                                                                                  gCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGc                    <  1:164127/149‑1 (MQ=255)
                                                                                                                                    ggCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGa                  <  3:144949/149‑1 (MQ=255)
                                                                                                                                        aGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCa              >  4:466695/1‑149 (MQ=255)
                                                                                                                                         taCTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCGTGATCTCCAGCGAAGGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAg             >  4:764780/2‑149 (MQ=255)
                                                                                                                                         gACTTATCGCCTGGGCGATCGCGTGGCAGTTCGCGTCGAAGCGGTTAAAAAGGACGAGCGCAAAAGCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAg             <  2:843944/149‑1 (MQ=255)
                                                                                                                                            ttATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTa          >  2:886521/1‑149 (MQ=255)
                                                                                                                                                 gCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTaaaaaa     <  4:378756/149‑1 (MQ=255)
                                                                                                                                                 gCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCCACGCAGGTaaaaaa     <  3:66765/149‑1 (MQ=255)
                                                                                                                                                  ccTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAg    <  3:466695/149‑1 (MQ=255)
                                                                                                                                                   cTGGGCGATCGCGTGGAAGGTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACTTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAgg   <  4:211229/149‑1 (MQ=255)
                                                                                                                                                    tGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAggc  >  3:469709/1‑149 (MQ=255)
                                                                                                                                                    tGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAggc  <  3:53661/149‑1 (MQ=255)
                                                                                                                                                    tGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAggc  <  3:536924/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAAAATCGACTTTAGCCTGATCTCCAGCGAACGCGCACCGCGCAACGTCGGTAAAACGGCGCGCGAGAAAGCGAAAAAAGGCGATGCAGGTAAAAAAGGC  >  NC_000913/4408610‑4408906

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