Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,558,291 C→T G8G (GGC→GGT rtcR → DNA‑binding transcriptional activator RtcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,558,2910CT98.5% 230.9 / ‑6.4 65G8G (GGC→GGTrtcRDNA‑binding transcriptional activator RtcR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/28);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACG  >  NC_000913/3558050‑3558535
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aTTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGt                                                                                                                                                                                                                                             >  1:1188502/1‑251 (MQ=255)
aTTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGt                                                                                                                                                                                                                                             <  2:1187664/251‑1 (MQ=255)
        tGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTgg                                                                                                                                                                                                                                                 <  2:810056/239‑2 (MQ=255)
        tGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGACCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGAAAACGTAGATTTCACGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGAGGCTTTTGGTTTTGt                                                                                                                                                                                                                                                 >  1:810056/1‑239 (MQ=255)
         ggTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTACGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTg                                                                                                                                                                                                                                    <  1:1133789/251‑1 (MQ=255)
                                     ttCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATg                                                                                                                                                                                                                             <  2:208130/230‑1 (MQ=255)
                                     ttCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATg                                                                                                                                                                                                                             >  1:208130/1‑230 (MQ=255)
                                         ctTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTc                                                                                                                                                                                                    >  2:934116/1‑251 (MQ=255)
                                           tttCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCa                                                                                                                                                                                                   <  2:625579/250‑1 (MQ=255)
                                              cGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGCTTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGACATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGACGGACGCGGCAGTCAGCGCTGGTCaacat                                                                                                                                                                                               >  2:577050/1‑248 (MQ=255)
                                                  gttTGCTGTCATGATAATTGCCACCGTCGTGCCAGACACTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTACGCGCTGGTCAAAATGGCg                                                                                                                                                                                           >  2:1179052/1‑251 (MQ=255)
                                                    ttGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTc                                                                                                                                                                                         >  1:646531/1‑251 (MQ=255)
                                                     tGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTCCGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTAGGCCGGGCGCGGCAGTCAGCTCTGGTCAAAATAGCGTc                                                                                                                                                                                         <  1:828158/250‑1 (MQ=255)
                                                        tGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGACGGGCGCGGCAGTCAGGGATGGTCAACATggactcagac                                                                                                                                                                                     >  2:878754/1‑243 (MQ=255)
                                                              gATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACActct                                                                                                                                                                                <  2:467830/250‑1 (MQ=255)
                                                                     gCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGc                                                                                                                                                                        >  2:971127/1‑251 (MQ=255)
                                                                       aaCCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCa                                                                                                                                                                      <  2:646531/251‑1 (MQ=255)
                                                                            tcgtGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCa                                                                                                                                                                   >  1:556557/1‑249 (MQ=255)
                                                                               tGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAAt                                                                                                                                                              >  1:924217/1‑251 (MQ=255)
                                                                                gCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATc                                                                                                                                                             <  2:556557/251‑1 (MQ=255)
                                                                                        aTTAAATAAGCAGCTTAATTTTTTAATTAATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTc                                                                                                                                                      >  1:1102548/1‑250 (MQ=255)
                                                                                              tAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTAAGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCatcgatc                                                                                                                                               >  2:819411/1‑251 (MQ=255)
                                                                                               aaGCAGCTTAGTTTTTTAATTCATTGCTTTAAAAAAGATTATCTCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCCATTTACCGATAATGAAATATAGTCTTTTATAAGGATATATAAGATGCGTAAAAAAGTGGCTTTTGCTTTTGTCGATACCGTACTGGATTATGCCGGGCGCGGCAGTAAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGACTCGTTGGTCatcgntcg                                                                                                                                              >  2:1155718/1‑251 (MQ=255)
                                                                                                    gcTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGTAATCTAAGATGCGTAAAACAGGGGCTTTTGGTTTTGTCGGTCCCGTACTGGATTATGCGGGGCGCGGCAGTCGGCGCTGGTCAAAATGGCGTCCGAAACTCTGTTTATGCCACAAAGAATGGTTGTTCatcgtcggcnggg                                                                                                                                         >  2:1148206/1‑242 (MQ=255)
                                                                                                            ttttAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAAttgtt                                                                                                                                  <  1:197577/250‑1 (MQ=255)
                                                                                                                     aTTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGcc                                                                                                                        <  1:971127/251‑1 (MQ=255)
                                                                                                                     aTTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGcc                                                                                                                        >  2:156833/1‑251 (MQ=255)
                                                                                                                      ttGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGccc                                                                                                                       >  2:397417/1‑251 (MQ=255)
                                                                                                                            tAAAAAAGATTATATCGTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGGTGGTCATCGATCGACTGGAATTGTTGCACGACGCCcgctcg                                                                                                                  <  2:1188502/250‑1 (MQ=255)
                                                                                                                                     ttATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTAt                                                                                                        >  2:921092/1‑251 (MQ=255)
                                                                                                                                         atCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTg                                                                                                     >  1:888080/1‑250 (MQ=255)
                                                                                                                                            tttACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAAc                                                                                                 >  2:251539/1‑251 (MQ=255)
                                                                                                                                            tttACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGAGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAAc                                                                                                 >  1:677680/1‑251 (MQ=255)
                                                                                                                                             ttACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGa                                                                                                    <  2:888080/247‑1 (MQ=255)
                                                                                                                                             ttACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACg                                                                                                >  2:1000530/1‑251 (MQ=255)
                                                                                                                                               aCGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCt                                                                                              >  1:42878/1‑251 (MQ=255)
                                                                                                                                                     gTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAAcg                                                                                        >  2:672488/1‑251 (MQ=255)
                                                                                                                                                        aCCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATCATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGGCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAAGTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACTCTTAAACGCGa                                                                                     <  1:1155718/251‑1 (MQ=255)
                                                                                                                                                            gagaTTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATc                                                                                 <  1:1000530/251‑1 (MQ=255)
                                                                                                                                                                        aaaGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCg                                                                     >  1:1142046/1‑251 (MQ=255)
                                                                                                                                                                                   aCCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAAc                                                             >  2:283804/1‑248 (MQ=255)
                                                                                                                                                                                   aCCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGa                                                          <  1:251539/251‑1 (MQ=255)
                                                                                                                                                                                      gATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTc                                                                                                                                                      <  2:196768/156‑1 (MQ=255)
                                                                                                                                                                                      gATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTc                                                                                                                                                      >  1:196768/1‑156 (MQ=255)
                                                                                                                                                                                          aTGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTGGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACTCTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGGTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAgtggtg                                                   <  1:1179052/251‑1 (MQ=255)
                                                                                                                                                                                            gAAATATCGTCTTTTCTAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAgtggtg                                                   >  2:170877/1‑249 (MQ=255)
                                                                                                                                                                                             aaaTATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGc                                                >  1:47268/1‑251 (MQ=255)
                                                                                                                                                                                                 atCGTCTTTTATAAGGATATCTAAGATGCGTAAAGCAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGtt                                             <  2:42878/250‑1 (MQ=255)
                                                                                                                                                                                                 atCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGcgctcgct                                                                                                          <  1:255651/189‑1 (MQ=255)
                                                                                                                                                                                                 atCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGcgctcgct                                                                                                          >  2:255651/1‑189 (MQ=255)
                                                                                                                                                                                                  tCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTg                                            >  1:45816/1‑250 (MQ=255)
                                                                                                                                                                                                   cGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTga                                           <  1:672488/250‑1 (MQ=255)
                                                                                                                                                                                                    gTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGttt                                                                       <  2:5720/221‑1 (MQ=255)
                                                                                                                                                                                                    gTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGttt                                                                       >  1:5720/1‑221 (MQ=255)
                                                                                                                                                                                                    gTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTgaga                                         <  1:934116/251‑1 (MQ=255)
                                                                                                                                                                                                          tataAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAACATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAAc                                   >  2:879689/1‑251 (MQ=255)
                                                                                                                                                                                                             aaGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCTTTGAGATTGAACTg                                 <  1:322565/250‑1 (MQ=255)
                                                                                                                                                                                                              aGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGc                                >  2:874890/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCAGAACCCGtt                      <  2:330437/251‑2 (MQ=255)
                                                                                                                                                                                                                        aGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTgg                     <  2:45816/251‑1 (MQ=255)
                                                                                                                                                                                                                         gATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTggg                    <  1:921092/251‑1 (MQ=255)
                                                                                                                                                                                                                         gATGCGTAAAACAGTGGCTTTTGGTTTTGACGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGGTCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCCGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTggg                    >  2:889829/1‑251 (MQ=255)
                                                                                                                                                                                                                              gTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTc               <  1:1148206/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAag           >  2:558082/1‑249 (MQ=255)
                                                                                                                                                                                                                                     aGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAgg        >  2:648501/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tttggttttGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTAc   <  1:874890/249‑1 (MQ=255)
                                                                                                                                                                                                                                            tttggttttGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACg  <  2:924217/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
ATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACG  >  NC_000913/3558050‑3558535

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