Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,407,045 C→T A131V (GCT→GTT)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,0450CT96.5% 1411.9 / 26.3 428A131V (GCT→GTT) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base C (6/9);  new base T (213/198);  total (219/209)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

Only 100 of 430 total aligned reads displayed.

GCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCC  >  NC_000913/4406895‑4407192
                                                                                                                                                      |                                                                                                                                                   
    cGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGTCTACGGCTTGCTGCGGGTTGAAGGGCGTAAAGTTGATTTGTATCTCTCCAGCGAGCAGATGAAAACATGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGAc                                                                                                                                                 <  4:1325966/151‑1 (MQ=255)
    cGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGa                                                                                                                                                  <  2:478781/150‑1 (MQ=255)
     gACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGAcc                                                                                                                                                >  5:939466/1‑151 (MQ=255)
      aCCTGTTGAAAGGTAGCGTTATTGGCCACGGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGGAAAGATGATTTGTATCTCTCCAGCGAGCAGATGTAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCg                                                                                                                                               <  3:970668/151‑1 (MQ=255)
        cTGGTGAAACGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGt                                                                                                                                              >  6:691134/1‑150 (MQ=255)
           gTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAg                                                                                                                                          <  6:1303517/151‑1 (MQ=255)
            tGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAgg                                                                                                                                         >  1:1460324/1‑151 (MQ=255)
             gAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGgt                                                                                                                                        <  5:1541298/151‑1 (MQ=255)
              aaaGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGgtc                                                                                                                                       <  8:1606228/151‑1 (MQ=255)
              aaaGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGCGTGTTGACCGTAAAGgtc                                                                                                                                       >  8:1213274/1‑151 (MQ=255)
                aGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGgtcgt                                                                                                                                     <  5:30562/151‑1 (MQ=255)
                 ggTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGgtcgtc                                                                                                                                    >  6:1374650/1‑151 (MQ=255)
                 ggTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGgtcgtc                                                                                                                                    <  5:857952/151‑1 (MQ=255)
                  gTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGtcgtcg                                                                                                                                   <  5:1176166/151‑1 (MQ=255)
                   tACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTcgtcgt                                                                                                                                  >  3:166592/1‑151 (MQ=255)
                    aCCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTg                                                                                                                                 <  3:1168884/151‑1 (MQ=255)
                      cGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGa                                                                                                                                >  5:167397/1‑150 (MQ=255)
                          aTTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGGGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAgcgc                                                                                                                           >  2:1621603/1‑151 (MQ=255)
                            tGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGt                                                                                                                         >  1:1356249/1‑151 (MQ=255)
                               ccACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTAtt                                                                                                                      >  4:485534/1‑151 (MQ=255)
                                cACCGTGATGGCTACGGGTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATATCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGGAttt                                                                                                                     >  2:716771/1‑150 (MQ=255)
                                cACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTg                                                                                                                     <  3:1508506/151‑1 (MQ=255)
                                      gttGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAATGATGTTTGGTATCTGTCCAGTGAGCAGCTGTAAACCTGCATTCTTGGCGAGCAGGTGCTGGCGCAGCCGCTGGGGGTGGTCCGTCAAGGTCGTCGTGATGCGCGTTTTGTCcgcg                                                                                                               <  4:239320/149‑1 (MQ=255)
                                         ggCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTAc                                                                                                            <  3:109891/151‑1 (MQ=255)
                                              cGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTg                                                                                                       <  3:799450/151‑1 (MQ=255)
                                               ggCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGc                                                                                                      >  5:1416204/1‑151 (MQ=255)
                                                 cTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGcc                                                                                                     >  5:617370/1‑150 (MQ=255)
                                                 cTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGcc                                                                                                     >  7:371098/1‑150 (MQ=255)
                                                   ttCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCaa                                                                                                   >  7:835849/1‑150 (MQ=255)
                                                     cTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCaaaa                                                                                                 >  6:1119337/1‑150 (MQ=255)
                                                         gggTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAcc                                                                                              >  2:1315382/1‑149 (MQ=255)
                                                         gggTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAcc                                                                                              >  8:422171/1‑149 (MQ=255)
                                                              gAATGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGa                                                                                       <  1:1042087/151‑1 (MQ=255)
                                                              gAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagccag                                                                                        <  1:1674411/150‑1 (MQ=255)
                                                              gAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagcc                                                                                          >  2:587428/1‑148 (MQ=255)
                                                              gAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGa                                                                                       <  7:164643/151‑1 (MQ=255)
                                                               aaGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagcc                                                                                          >  1:582880/1‑147 (MQ=255)
                                                                aGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGAtt                                                                                     <  3:832598/151‑1 (MQ=255)
                                                                 gggCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGAttg                                                                                    >  2:1213223/1‑151 (MQ=255)
                                                                    cGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGAttgttg                                                                                 >  3:1436180/1‑151 (MQ=255)
                                                                       aaaGATGATTTGTATCTCTCGAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGt                                                                               <  1:699951/150‑1 (MQ=255)
                                                                        aaGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCg                                                                             >  7:1039853/1‑151 (MQ=255)
                                                                        aaGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCg                                                                             >  3:472099/1‑151 (MQ=255)
                                                                          gatgatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCt                                                                           >  2:914141/1‑151 (MQ=255)
                                                                           atgatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTa                                                                          <  4:762436/151‑1 (MQ=255)
                                                                           atgatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTa                                                                          <  3:505231/151‑1 (MQ=255)
                                                                              atTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACtt                                                                       >  2:378607/1‑151 (MQ=255)
                                                                              atTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACt                                                                        <  6:417899/150‑1 (MQ=255)
                                                                                ttGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACttt                                                                      <  4:906428/150‑1 (MQ=255)
                                                                                ttGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACttt                                                                      <  8:1634248/150‑1 (MQ=255)
                                                                                 tGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTAc                                                                    >  4:1448368/1‑151 (MQ=255)
                                                                                   tATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGAGCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTAccc                                                                  >  5:431243/1‑150 (MQ=255)
                                                                                   tATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCg                                                                  >  8:874147/1‑151 (MQ=255)
                                                                                    aTCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGGCAGATTGTTGGTCGCTACTTTACCg                                                                  <  8:249478/150‑1 (MQ=255)
                                                                                    aTCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCg                                                                  >  8:1319671/1‑150 (MQ=255)
                                                                                    aTCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCg                                                                  <  6:1667774/150‑1 (MQ=255)
                                                                                    aTCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCg                                                                  >  1:1211189/1‑150 (MQ=255)
                                                                                       tctcCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCGAGATTGTTGGTCGCTACTTTACCGAAGc                                                              >  4:1646880/1‑151 (MQ=255)
                                                                                        ctcCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGa                                                                 <  4:999003/147‑1 (MQ=255)
                                                                                         tcCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCgg                                                            >  4:1150963/1‑151 (MQ=255)
                                                                                          ccAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCgg                                                            <  3:682242/150‑1 (MQ=255)
                                                                                          ccAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCg                                                             <  7:975355/149‑1 (MQ=255)
                                                                                          ccAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCg                                                             <  1:1526671/149‑1 (MQ=255)
                                                                                          ccAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCAGCGAACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCggg                                                           >  8:1577510/1‑151 (MQ=255)
                                                                                           cAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGc                                                          >  8:811941/1‑151 (MQ=255)
                                                                                            aGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCg                                                         >  5:529646/1‑151 (MQ=255)
                                                                                             gCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGTCCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCTAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGt                                                        <  3:261407/151‑1 (MQ=255)
                                                                                              cGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTc                                                       <  8:688359/151‑1 (MQ=255)
                                                                                               gAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCg                                                      >  6:453849/1‑151 (MQ=255)
                                                                                                   aGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCggggt                                                  >  6:917158/1‑148 (MQ=255)
                                                                                                     aTGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTtg                                                <  2:1060379/151‑1 (MQ=255)
                                                                                                     aTGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTtg                                                >  5:702144/1‑151 (MQ=255)
                                                                                                       gAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTtgtg                                              >  4:677134/1‑151 (MQ=255)
                                                                                                         aaaCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGt                                            <  6:83319/151‑1 (MQ=255)
                                                                                                           cccTGCAGTCATGGCGATCAGGTGCTGGCGCAGCCGCTGGGTGTTGACCGTAAAGGGCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGGTGGTCGCTACTTTACCGAAGCGGGCGGCGGCTTTGTGGTTc                                          <  7:1205949/150‑1 (MQ=255)
                                                                                                           aCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTc                                          >  8:1067157/1‑151 (MQ=255)
                                                                                                           aCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTc                                          >  2:1504443/1‑151 (MQ=255)
                                                                                                           aCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTc                                          >  7:1086359/1‑151 (MQ=255)
                                                                                                           aCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTc                                          <  5:778067/151‑1 (MQ=255)
                                                                                                              tGCATTCATGGCGATCAGGTGCTGGCTAAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTg                                       >  5:686133/1‑150 (MQ=255)
                                                                                                                    cATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACa                                 <  1:844205/151‑1 (MQ=255)
                                                                                                                        gCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTGCCTGACGACAGCCg                             <  1:960156/151‑1 (MQ=255)
                                                                                                                         cGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGt                            <  5:135479/151‑1 (MQ=255)
                                                                                                                         cGATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGt                            <  5:808325/151‑1 (MQ=255)
                                                                                                                          gATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTc                           >  3:1550561/1‑151 (MQ=255)
                                                                                                                          gATCAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCATGACGACAGCCGTc                           <  7:1461965/151‑1 (MQ=255)
                                                                                                                             cAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTg                         >  1:1106173/1‑150 (MQ=255)
                                                                                                                             cAGGTGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTg                         >  7:1341694/1‑150 (MQ=255)
                                                                                                                              aGGTGCTGGCTCAGGCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCTCGTACTGGTGCGAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAg                       <  1:490129/151‑1 (MQ=255)
                                                                                                                                gTGCTGGCTCAGCCGCTGGGTGTTGAACTAAAAGGTCGTCGTGAAGATAGTATTGTCCGCGTACTGGTGCCAAAAACCGGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGCCAGCCGTCTGAGCt                     >  4:1110140/1‑151 (MQ=255)
                                                                                                                                 tGCTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCtt                    >  1:1274088/1‑151 (MQ=255)
                                                                                                                                  gCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCtt                    <  3:1321386/150‑1 (MQ=255)
                                                                                                                                   cTGGCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCg                  >  5:348628/1‑151 (MQ=255)
                                                                                                                                     ggCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGat                >  5:1371562/1‑151 (MQ=255)
                                                                                                                                      gCTCAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGat                >  3:1028664/1‑150 (MQ=255)
                                                                                                                                       atCAGCCGCTGGGTGTAGACCTTAAAGTTAGACGTGAAGCGCGTATTGTCCGCGTACTGTGGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTGTGCTTTGTGGTTCCTCACTACATCCGTCTGAGCTTCGatat              <  6:148785/150‑1 (MQ=255)
                                                                                                                                         cAGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACCGCCGTCTGAGCTTCGATATCt            >  8:363352/1‑151 (MQ=255)
                                                                                                                                          aGCCGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCtt           >  8:619334/1‑151 (MQ=255)
                                                                                                                                             cGCTGGGTGTTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAAt        >  7:512734/1‑151 (MQ=255)
                                                                                                                                             cGCTGGGTGTTGACCGTAAAGGTCGTCGGGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAAt        <  6:1205603/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCC  >  NC_000913/4406895‑4407192

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