Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,408,287 A→C Q545P (CAG→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,2870AC98.7% 241.5 / ‑3.8 76Q545P (CAG→CCG) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (34/41);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGT  >  NC_000913/4408054‑4408519
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
tatCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCgg                                                                                                                                                                                                                                        <  1:627617/236‑1 (MQ=255)
  tctcGGCCGCGCGTTTCGTTAAAAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGGCAAGCCGAGCACCGAAGCCATTCCCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCg                                                                                                                                                                                                                                      <  2:199837/236‑1 (MQ=255)
     cggcggCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGCAGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAACGCTGCAAACCATGCTGCTGCGCTCGATGAAAACGGCGAtt                                                                                                                                                                                                                                   >  1:529448/1‑236 (MQ=255)
        cggcgCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGGGACTACCCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACACTGCTGCTGCGCGCGATGAAACAGGCGCTTTa                                                                                                                                                                                                                                 >  1:695860/1‑235 (MQ=255)
         ggcgCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTAc                                                                                                                                                                                                                                >  1:330492/1‑235 (MQ=255)
           cgcgTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGAt                                                                                                                                                                                                                             >  2:767291/1‑236 (MQ=255)
           cgcgTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGAACGTCCTGATGCAGAAATGCTGCAAACCAAGACGCATCGCTCGCTTGAACCGGCGATTTAAGag                                                                                                                                                                                                                             >  2:680988/1‑235 (MQ=255)
                     tGAGAAAGCGAAAGATCCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAAcc                                                                                                                                                                                                                   <  1:519344/236‑1 (MQ=255)
                           aGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTc                                                                                                                                                                                                             >  1:38244/1‑236 (MQ=255)
                                         cACTGTGCCGTTGTCACGACAAGCCGAGCACGGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTgg                                                                                                                                                                                                 <  2:529448/234‑1 (MQ=255)
                                             gTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGc                                                                                                                                                                                                >  1:511461/1‑231 (MQ=255)
                                             gTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGc                                                                                                                                                                                                >  1:483658/1‑231 (MQ=255)
                                             gTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTg                                                                                                                                                                                            <  1:20296/235‑1 (MQ=255)
                                             gTTCCGTATCCACGGCAAGCCGAGCAACGAAGCGATTCCCTAGTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTg                                                                                                                                                                                            <  2:184556/235‑1 (MQ=255)
                                                      tCACGACAAGCCGAGCACCGAAGCGATTACCTCTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATg                                                                                                                                                                                  >  1:369933/1‑235 (MQ=255)
                                                        aCGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATgcg                                                                                                                                                                                <  2:18096/236‑1 (MQ=255)
                                                          gACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCa                                                                                                                                                                              <  1:767291/236‑1 (MQ=255)
                                                                   aGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTgg                                                                                                                                                                                                 <  1:467550/208‑1 (MQ=255)
                                                                          aaGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTAAGATCCAGAAAACCGTGGTCACTTTTGGCCGGCATTGCa                                                                                                                                                                                          >  2:467550/1‑208 (MQ=255)
                                                                                       tttCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAg                                                                                                                                                 <  1:327828/236‑1 (MQ=255)
                                                                                          ccGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGAcc                                                                                                                                              <  1:130217/236‑1 (MQ=255)
                                                                                          ccGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGAcc                                                                                                                                              <  1:580567/236‑1 (MQ=255)
                                                                                             ttCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCTAGGGGACTTTGGCCTGGGATTGCAGTCCCATGCGCACTTTAATTACCCGCTTCTTCGTAATccgagccgga                                                                                                                                           >  2:145300/1‑228 (MQ=255)
                                                                                             ttCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGCGGTCCCTTTGGCCTGGCATTTCAGTCCTATGCGCACCTTACTTCGCCGATTTGCCGTTAttccaaccta                                                                                                                                            >  2:599899/1‑226 (MQ=255)
                                                                                               cAGTGCTGGCGTAGCTGGGGCTGGACCTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACg                                                                                                                                         <  1:185228/236‑1 (MQ=255)
                                                                                                       gCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGTGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGCAACCGGCGATTTa                                                                                                                                                                                                                                 >  2:702312/1‑140 (MQ=255)
                                                                                                       gCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTa                                                                                                                                                                                                                                 <  1:702312/140‑1 (MQ=255)
                                                                                                       gCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATtcgtcg                                                                                                                                                         <  2:16568/212‑1 (MQ=255)
                                                                                                       gCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATtcgtcg                                                                                                                                                         >  1:16568/1‑212 (MQ=255)
                                                                                                        ccgAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACcg                                                                                                                                 <  1:384019/233‑1 (MQ=255)
                                                                                                        cGGAGTGGGGGCTGGAACTGCCGGGCGGTACCAAGCCGGACCCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACcg                                                                                                                                 <  2:244767/235‑1 (MQ=255)
                                                                                                           aGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCgcgc                                                                                                                              <  2:190079/235‑1 (MQ=255)
                                                                                                            gCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTCATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGCTATTCCGACCTGACGCTGCACCGCGCCa                                                                                                                            >  2:73383/1‑236 (MQ=255)
                                                                                                                       aaCTGCCGGGCGGTAACAAGCCGGAACCGCGTGCCTACGCGGAGCTGCTGGCGTCGGTTGGCGATCGTCCTGATGCAGAAATGCAGCAAACCATGCTTCTCCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCCCTTTGGCCTGGCATTGCATTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGCGGCTGCACCGCGCCATTTAata                                                                                                                     >  2:179693/1‑232 (MQ=255)
                                                                                                                                      aCAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTCATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCa                                                                                                   <  1:73383/235‑1 (MQ=255)
                                                                                                                                        aaGCCGGAACCGCGTGACTACGGGGAGCTGCGGGAGTCGGTTCCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAgg                                                                                                 <  2:115757/235‑1 (MQ=255)
                                                                                                                                        aaGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAgg                                                                                                 <  1:196894/235‑1 (MQ=255)
                                                                                                                                              gAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAACATTCTGCTCGCGAAAGAGCAAGGGGAtaa                                                                                          >  2:499212/1‑234 (MQ=255)
                                                                                                                                              gAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAgg                                                                                                 <  2:483658/229‑1 (MQ=255)
                                                                                                                                              gAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATc                                                                                           >  1:241936/1‑235 (MQ=255)
                                                                                                                                                 ccGGGTGACTACGCGGAGCGGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCAt                                                                                            <  2:697687/231‑1 (MQ=255)
                                                                                                                                                    cgTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCa                                                                                     >  2:18868/1‑235 (MQ=255)
                                                                                                                                                    cgTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACACTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAAGGGCATcctggga                                                                                     >  2:258661/1‑229 (MQ=255)
                                                                                                                                                    cgTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCCATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAACAGCAGGTGCATCAGGGCa                                                                                     >  2:150768/1‑235 (MQ=255)
                                                                                                                                                      tGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAaca                                                                                  <  1:12705/236‑1 (MQ=255)
                                                                                                                                                            cgcgGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACt                                                                             >  1:433009/1‑235 (MQ=255)
                                                                                                                                                            cgcgGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACt                                                                             >  2:11158/1‑235 (MQ=255)
                                                                                                                                                                      ctgGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGc                                                                                                    >  1:53403/1‑202 (MQ=255)
                                                                                                                                                                      ctgGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACGTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGc                                                                                                    <  2:53403/202‑1 (MQ=255)
                                                                                                                                                                       tgGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCggcggc                                                                 <  1:18868/236‑1 (MQ=255)
                                                                                                                                                                        ggAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCggcggc                                                                 <  2:241936/235‑1 (MQ=255)
                                                                                                                                                                           gTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCggcggc                                                                 >  2:515393/1‑232 (MQ=255)
                                                                                                                                                                              ggTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCat                                                           >  1:436854/1‑235 (MQ=255)
                                                                                                                                                                               gTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGCCGGCTACCatta                                                         >  1:178694/1‑236 (MQ=255)
                                                                                                                                                                                    cGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCg                                                     <  1:11158/235‑1 (MQ=255)
                                                                                                                                                                                          tCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCa                                                            >  1:47316/1‑222 (MQ=255)
                                                                                                                                                                                          tCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCa                                                            <  2:47316/222‑1 (MQ=255)
                                                                                                                                                                                                     aaaTGCTGCAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGCGGCATCAGGGCAACACCACTGAAACCGGCGGCTCCCATTATTCGCTGGAagcaaggct                                       >  2:559224/1‑225 (MQ=255)
                                                                                                                                                                                                      aaTGCTGCAAGCCATGCTGCTTCGCTTGAGGAAACCGGCGATTTCCGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  2:264616/236‑1 (MQ=255)
                                                                                                                                                                                                         gctgcAATCCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  1:145300/233‑1 (MQ=255)
                                                                                                                                                                                                         gctgcAAGCCATGCTGCTGCGCTTGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  1:258661/233‑1 (MQ=255)
                                                                                                                                                                                                         gctgcAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  2:38244/233‑1 (MQ=255)
                                                                                                                                                                                                         gctgcAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  1:515393/233‑1 (MQ=255)
                                                                                                                                                                                                         gctgcAAACCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAAc                                  <  1:599899/233‑1 (MQ=255)
                                                                                                                                                                                                           tgcAACCCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGt                             <  2:436854/236‑1 (MQ=255)
                                                                                                                                                                                                                aCCATGCTGCTGCGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCa                        <  1:680988/236‑1 (MQ=255)
                                                                                                                                                                                                                    tgctgctgcGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATtcgtc                                                                                                                                                          <  1:429384/102‑1 (MQ=255)
                                                                                                                                                                                                                    tgctgctgcGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATtcgtc                                                                                                                                                          >  2:429384/1‑102 (MQ=255)
                                                                                                                                                                                                                    tgctgctgcGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGCGCAGGGGCATCAGGGCAACACCAGGTAAACCGGCGGCTACCATTATTCGATGGACAAGATGCTGGAACTggtctagaacggt                    >  2:431076/1‑225 (MQ=255)
                                                                                                                                                                                                                      ctgctgcGCTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCCCTGTTc                  >  2:194829/1‑236 (MQ=255)
                                                                                                                                                                                                                        gctgcGCTCGATGAAACCGGCGATTTACGATCCAGAAACCCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGa                <  2:351457/236‑1 (MQ=255)
                                                                                                                                                                                                                            cgcTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCAt                                                                                            >  1:289187/1‑156 (MQ=255)
                                                                                                                                                                                                                            cgcTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCAt                                                                                            <  2:289187/156‑1 (MQ=255)
                                                                                                                                                                                                                            cgcTCGATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGc            >  1:543415/1‑236 (MQ=255)
                                                                                                                                                                                                                                 gATGAAACCGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAAc       >  2:515524/1‑236 (MQ=255)
                                                                                                                                                                                                                                      aaCCGGCGATTTACGATCCAGAAACCCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAAcgtcgt  <  2:330492/236‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
TATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGT  >  NC_000913/4408054‑4408519

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