Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,053 G→C H5D (CAT→GAT)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,0530GC98.0% 167.2 / ‑6.0 51H5D (CAT→GAT) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (22/28);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCA  >  NC_000913/4616807‑4617286
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
gccgttCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAATACATAGCGCGGTTTGTACGGCGCATGACCGTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCTCGGGGCTGGTGACAATTTGCTGGCAACGTGGTTGCAGTGCATGTTCATCAGaa                                                                                                                                                                                                                                       <  2:95458/247‑1 (MQ=255)
      cGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGATTGTACGGCGCATGACCTTCAAACCTATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCCGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGCGCCAGAACATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGc                                                                                                                                                                                                                                 >  1:96175/1‑251 (MQ=255)
       ggAACCGTTCGCCAGAAAAAGGGCGTAATCGGGTAAGACCTAGCGCGGTTTGTCCGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGc                                                                                                                                                                                                                                 >  2:613300/1‑250 (MQ=255)
         aaCCGTTCGCCAGAAAACGGTCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCAGGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATGTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtc                                                                                                                                                                                                                              <  2:338963/251‑1 (MQ=255)
           ccGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtctt                                                                                                                                                                                                                            <  1:95428/251‑1 (MQ=255)
            cGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtcttc                                                                                                                                                                                                                           >  1:216781/1‑251 (MQ=255)
             gTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtcttct                                                                                                                                                                                                                          >  1:376334/1‑251 (MQ=255)
                      aaaaCGGCCGTAATCGGGTAAGACCTAGCACGGTTTGTACAGCGCATGACCTTCAAACAAATCGCAGATTACACCTTGATCCAGCGCGAGGCGGGCTACGGCAGGAAGCTGTGGGTAAGGACGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGACACGGGCTAAGCACTGGGCTGGTGAAAATTTGCTGGCCACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCtttt                                                                                                                                                                                                                 >  2:36154/1‑251 (MQ=255)
                               gTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTTCACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAAGGGCGCTTCTGCTCCGGGCTAAGCACTTGGCTGGTGACAATTTGCTTGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAg                                                                                                                                                                                                        <  1:382534/251‑1 (MQ=255)
                               gTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAg                                                                                                                                                                                                        <  2:413350/251‑1 (MQ=255)
                                 aaTCGGGTAAGACATAGCGCTGTTTGTACGGCGCAGTACCTTCTAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGTGTTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGGTCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAgg                                                                                                                                                                                                       <  2:512024/250‑1 (MQ=255)
                                                  cgcgGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGATGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGATGAGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTa                                                                                                                                                                                        >  1:316151/1‑248 (MQ=255)
                                                   gcgGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCAGCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCGACGTTGTTTCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTaagaa                                                                                                                                                                                    <  2:402673/251‑1 (MQ=255)
                                                     ggTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTGGCTGGCAACGTGGGTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGt                                                                                                                                                                                  <  1:293802/251‑1 (MQ=255)
                                                     ggTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAAGAAGCTGTGGGTAAGGCAGAGTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTTGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCGTCTTTTCCTGTTATGCCGAAGGTGATGCGCCATTGTAAGAAGt                                                                                                                                                                                  <  2:153584/251‑1 (MQ=255)
                                                             cGGCGCATGAACTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTACTGTGCCAGAAAATGGCGCTTCTGCTCCGGGATAAGCACTGGGCTGGTGACAATTTGCTGGAAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGTAGTTTCGTgtt                                                                                                                                                                          >  1:7900/1‑249 (MQ=255)
                                                                 gcATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTc                                                                                                                                                                      >  1:132616/1‑251 (MQ=255)
                                                                     gACCTTCAAACATATCGCAGATTACTCCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCAACAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACt                                                                                                                                                                   <  1:115832/250‑1 (MQ=255)
                                                                      aCCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAGGAAGTTTCGTGATGTTCACttt                                                                                                                                                                 <  2:452517/251‑1 (MQ=255)
                                                                               cATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGGAAGAAGTTTCGTGATGTTCACTTTGATCCTGAt                                                                                                                                                        <  1:187017/251‑1 (MQ=255)
                                                                               aatatCGAAGATTACACCTTCATCCAGCGGGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGTCTAAGCACTGGGCTGGTGACAAGTTGCTGGCAACGTTGTTGCAGTGCAGTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGTTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGAt                                                                                                                                                        <  2:69662/250‑1 (MQ=255)
                                                                                            acacCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGccaccc                                                                                                                                           >  1:518666/1‑250 (MQ=255)
                                                                                                       cAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattc                                                                                                                                >  1:181362/1‑251 (MQ=255)
                                                                                                       cAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattc                                                                                                                                <  2:79525/251‑1 (MQ=255)
                                                                                                          cgcgcgGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattcatt                                                                                                                             <  1:39762/251‑1 (MQ=255)
                                                                                                          cgcgcgGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattcatt                                                                                                                             >  1:327131/1‑251 (MQ=255)
                                                                                                                gcggGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGAGGAGTTTGCCACCACTGACGCAGTCATTTGAAAg                                                                                                                       <  2:132616/251‑1 (MQ=255)
                                                                                                                    gCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGaa                                                                                                                   <  1:615887/251‑1 (MQ=255)
                                                                                                                        cGGCAGGAAGCTGTGGGTAAGGCTGATTGTTTTCTTCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGGTGCAGGTCAGTTTCAGCAGAAGTGGGCATCGTCGTTTCCGTTTATGCCGAAGGTGATGCGCCATTTTaaga                                                                                                                                                                                     <  2:609536/181‑1 (MQ=255)
                                                                                                                        cGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCCGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTaaga                                                                                                                                                                                     >  1:609536/1‑181 (MQ=255)
                                                                                                                          gCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCGTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTAttt                                                                                                             <  1:404984/251‑1 (MQ=255)
                                                                                                                                       ggTAAGGCAGATTGTTTTCGGCTTCCAGTGCCAGAAAATTGCGCTGCTGCTCCGGGCGAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTGTCATCAGAAGTGGGCATCTTCTTTTCCTGTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTGTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCa                                                                                                <  1:613556/251‑1 (MQ=255)
                                                                                                                                                 aTTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGAGCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGAt                                                                                      <  1:613300/251‑1 (MQ=255)
                                                                                                                                                           gCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTg                                                                             <  1:125483/250‑1 (MQ=255)
                                                                                                                                                           gCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGt                                                                            >  1:565572/1‑251 (MQ=255)
                                                                                                                                                                    gCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGAtt                                                                   >  1:181774/1‑251 (MQ=255)
                                                                                                                                                                    gCCAGAAAAGGGCGTTCCCGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTCGCGACGTTGTTGCAGAGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTATAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGCCGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGCTGCAAACTTGTAAGTAGAtt                                                                   >  2:231492/1‑251 (MQ=255)
                                                                                                                                                                       aGAAAATGGCGCTTTTGCGCCGGGCTAAGCACTGTGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTTATCAGAAGTTGGCATCTTCTTTTCCGTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTGCTCTTTGATCCTGATGAGGGTGTCACTACTTACGCATTCATTTGAAAGTGAATGATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGTTTTcc                                                                <  2:628536/251‑1 (MQ=255)
                                                                                                                                                                                     gtGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGTCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCCCTGACGCAGTCTTTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATtgtgctgc                                                   >  2:211869/2‑246 (MQ=255)
                                                                                                                                                                                          ccGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTatat                                               <  2:216781/249‑1 (MQ=255)
                                                                                                                                                                                                tAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAgtgt                                       >  1:390905/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTACAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAgtgtg                                      <  1:159289/251‑1 (MQ=255)
                                                                                                                                                                                                    cACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTg                                   >  1:17534/1‑251 (MQ=255)
                                                                                                                                                                                                        gggCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCg                                 <  1:570671/249‑1 (MQ=255)
                                                                                                                                                                                                              gTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGt                             >  2:155285/1‑247 (MQ=255)
                                                                                                                                                                                                               ggACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTTTTGCGGAGTAGATg                        <  1:161969/250‑1 (MQ=255)
                                                                                                                                                                                                                        gCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATAc               >  1:21068/1‑251 (MQ=255)
                                                                                                                                                                                                                          tGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCAGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACt              >  1:78825/1‑250 (MQ=255)
                                                                                                                                                                                                                               ccGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGCTGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTAGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTTTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTCACaaa        <  2:145926/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gttgttGCAGTGCATTTTCATCAGAAGTGGGAATCTTCTTTTCCTTTTATGCCGAAGGTACTGCGACATTGTAAGAAGTTTCGTGATGTTCACCTTGATACTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAATTAGATTTCCTTAATTGAGATGTATATCGAAGTGTGTTGCGGAGTAGCTGTTAGAATACTAACAAACt      >  2:91159/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCa  >  1:257466/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
GCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCA  >  NC_000913/4616807‑4617286

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