Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,384,808 C→T K305K (AAG→AAA)  yfbK ← IPR002035/DUF3520 domain‑containing protein YfbK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,384,8080CT95.3% 167.7 / ‑5.9 64K305K (AAG→AAA) yfbKIPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (35/26);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 8.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAT  >  NC_000913/2384562‑2385044
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTTCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGGCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTt‑c                                                                                                                                                                                                                                            <  2:240014/250‑1 (MQ=255)
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTt‑cg                                                                                                                                                                                                                                           <  1:855375/251‑1 (MQ=255)
       aCTTTCTGCTCTTCAGAGATGGTAACAATGTAGCTGTAGTTGCCGTTACCAACATCTGCAATTCGCACCATCATCGCCTCGGTGTAATTGCGATTCCCCACGCCAAACGTCGACAGAGTGACACCAGTCTCCCGCTGTTTTTTGACCATTGATTCTATCTATTTTGGAGCGTCAATGCCAGCGTTAAAGTCACCGCCAGTGGCGCAAAAAATGCGATTGATGCCGCCCTGAATAATCCCTTT‑CTTCGCct                                                                                                                                                                                                                                     <  2:537101/250‑1 (MQ=255)
            cTGCGCTTCAGAGGGGGTATCACTGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTAGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACAATTGATTAAATCGATTTTGGATCGTCAATGCCAAAGTTAAAGTCACCGTCAGTGGCTAATATAATGCGATTGATGACGCCCTTAATAAACCCTTT‑CGTCGCctgctgc                                                                                                                                                                                                                                >  2:561104/1‑249 (MQ=255)
            cTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGAt                                                                                                                                                                                                                               >  1:165564/1‑251 (MQ=255)
               cgcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATa                                                                                                                                                                                                                              >  1:108015/1‑249 (MQ=255)
               cgcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATa                                                                                                                                                                                                                              >  1:637858/1‑249 (MQ=255)
                gcTTCAGAGAGGGAATCAATGTAGCTGTAGTTGCCGTTACCACCATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACCGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTT‑CGTCGCCTGCTGATAA‑Gc                                                                                                                                                                                                                          >  2:359351/1‑251 (MQ=255)
                    cAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGt                                                                                                                                                                                                                      >  1:584530/1‑251 (MQ=255)
                        gagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTcc                                                                                                                                                                                                                   >  1:853928/1‑250 (MQ=255)
                                  aTGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCAcc                                                                                                                                                                                                        <  1:435665/251‑1 (MQ=255)
                                                     tACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCg                                                                                                                                                                                     >  1:566776/1‑251 (MQ=255)
                                                       ccAACATCGGCAATTCGCACCATCATTGCCCTGTTGTAATTGCTATTCCCCTCGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTAATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGc                                                                                                                                                                                   <  1:576928/251‑1 (MQ=255)
                                                        cAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCa                                                                                                                                                                                                    >  2:546955/1‑233 (MQ=255)
                                                        cAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCa                                                                                                                                                                                                    <  1:546955/233‑1 (MQ=255)
                                                                 cAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCgc                                                                                                                                                                          >  1:32784/1‑249 (MQ=255)
                                                                           catcatTGCCTCTTTGTAATTGCTATCCCCCCCGCCAACCGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCgg                                                                                                                                                               <  1:550782/251‑1 (MQ=255)
                                                                              catTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑TCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCAt                                                                                                                                                            >  1:800542/1‑251 (MQ=255)
                                                                                       gTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCg                                                                                                                                                   >  1:398427/1‑251 (MQ=255)
                                                                                       gTTGTAATTACTATTACACACGCCAAACGTCGACAGAGTAACACAAGACTCCAGCTGTTTTTTGACCATTGATTCAATCGATTTAGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGCAGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAACCCGGCACCGCCATTGGTACTGACTTCGACATCCAGCGGATCAATTGCGGCATTCATTTCCg                                                                                                                                                   >  2:44134/1‑251 (MQ=255)
                                                                                          gTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCt                                                                                                                                                <  1:652253/251‑1 (MQ=255)
                                                                                                  tATTCCCCCCGCCAAACGTCGACAGAGTAACCCCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCAAACGTTAAAGTCCCCGTCAGTGGCTAATAAAAGGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGGCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGa                                                                                                                                           <  1:362430/248‑1 (MQ=255)
                                                                                                  tATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTAATTTCCGCCTTATGa                                                                                                                                           >  2:362430/1‑248 (MQ=255)
                                                                                                   atcccccACGCCAAACGCCGAAAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGAGGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTc                                                                                                                                        <  1:354794/247‑1 (MQ=255)
                                                                                                     tCCCCACGCCAACCGAAGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATCCAATCGATTTTGGAGCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAAGAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCGGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTccc                                                                                                                                      <  1:300386/250‑1 (MQ=255)
                                                                                                        ccACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCCGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCgg                                                                                                                                    >  2:864748/1‑249 (MQ=255)
                                                                                                           cGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                                >  1:235725/1‑250 (MQ=255)
                                                                                                            gCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                                >  2:318967/1‑249 (MQ=255)
                                                                                                                 aCGTCGACAGAGTAACAGCAGCCTCCCGCTGGTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTAACCGTCAGTGTCTAATAAAATGCGCTTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCAATTGGTCCTGCCTTCGGCATCCAGCGCATCAATTGCGGCATTAATTTCCGCCTTATTTCTCCCgtctatatattgt                                                                                                                         >  2:166062/1‑239 (MQ=255)
                                                                                                                  cGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgc                                                                                                                         >  1:188506/1‑250 (MQ=255)
                                                                                                                    ccgaagtatgacacccaaccccccGCTGTTTTATGCCCATTGTTTCAACCGATTGTGGTTCGTAATTCCCAACTTTAATGCCACCGTCGGTGGTTAATAAAAGGCAATTGTTGCCCCCCTTAATAACCCCCTT‑CCACGCCTGTTGATAA‑GCCAGTTCTGGCCCGGCCCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaat                                                                                                                      <  2:258678/230‑1 (MQ=255)
                                                                                                                     caacagagaaaccccagcaccccGGTGTTTTTTGCCAATTGATTCAATCGATTTTGGATCTGCAATGCCAACTTCAAAGTCTCCGTCAGTGGCTAATAAAATGCGATGGAGGCCGCCCATAATAAACCCTTT‑CGGCGCCTGTGGATAA‑GCCTGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     <  2:275908/231‑1 (MQ=255)
                                                                                                                     cGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     >  1:11136/1‑251 (MQ=255)
                                                                                                                     cGACAGAGTAACACCAGACTCCCGCTGTTTTGTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCCCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     <  1:8308/251‑1 (MQ=255)
                                                                                                                      gACAGATTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTAAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgc                                                                                                                    <  1:725015/251‑1 (MQ=255)
                                                                                                                        cAGAGTAACACCAGTCTCCCTCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTAAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTTCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgcaa                                                                                                                  <  1:456057/251‑1 (MQ=255)
                                                                                                                                             cTGTGTTTTGACCTGTTATTCAACAGATTTTGGATCGTCAATGCAAACGTAAAATCCACCGCCAGTGGCAAAAAAAAGGCGATTGATGCCCCCCTTAA‑AAACCCCTTTCGCCGCCTGCTGATAA‑GCTAGTCCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTa                                                                                               <  2:188506/249‑1 (MQ=255)
                                                                                                                                                ttttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTc                                                                                           >  2:265297/1‑250 (MQ=255)
                                                                                                                                                tttttggcccATTGATTCAATCGTTTGGGGATCGGCAATGCAAACTTTAAAGTCACCGTCAGTGGAAAATAAAATGCGATTGATGCCGCCCTTAA‑AAACCCCTTTCGTCTCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTGCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTc                                                                                           <  1:44134/242‑1 (MQ=255)
                                                                                                                                                       aCCATTGATTCAATCGATTTTGGATAGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGCTTGATGCCGCCCTTAATAAAACCTTT‑CGTAGCCTGCTGATAA‑GCCAGTTCCAGCCCTGCACCGCCATTGGTACTGCCTTCGGCATACAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCACGGAGATAGAAGGCAATGCAATACGGGGTTCGCCATAGTAGGTCACGATGGc                                                                                   >  2:369114/1‑251 (MQ=255)
                                                                                                                                                             gATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTAATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCACTGtt                                                                             >  1:792949/1‑251 (MQ=255)
                                                                                                                                                               ttCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTAt                                                                           >  1:183773/1‑251 (MQ=255)
                                                                                                                                                                 cAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATACAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATcc                                                                         >  1:74796/1‑251 (MQ=255)
                                                                                                                                                                            ggATCGTCAATGCCAACGTTAAAGTCACCGTCAGAGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGa                                                                >  1:671887/1‑249 (MQ=255)
                                                                                                                                                                                   aaaTGCCAACGTTAAAGCCACCGTAAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTg                                                       <  1:866657/250‑1 (MQ=255)
                                                                                                                                                                                          aaCGTTAAAGTCCCCGACAGTGGCTAAAAAAAGCCGATTGATGCCGCCTTTAACAAGCCCTTT‑CGTCGCCGGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTAACCTGCTCACGAAGTTCTTTGACCaata                                                <  1:87398/251‑1 (MQ=255)
                                                                                                                                                                                            cGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCaataat                                              >  1:144685/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGTCACCGTCAGAGGCTAATAAAATGAGATTGCTGCCGCCCTTAATAACCACTTT‑CGTCGCATGCTGACAC‑GACAGATCCAGCCCGGCACCGCCATTGGTAATGCCTTCGGCAACCAACGACTCAATTGCGGCATAAATTTCCGCCTTATGCCTCCGAGAGATAGAAGGCAATGCAATACGGGAGTCGCCCGCGTCGGTGAAGAAGGCACAGTTTTACTGCTCACGAAGTTCTTTGCCCAATAATTTCa                                          >  2:552245/1‑250 (MQ=255)
                                                                                                                                                                                                      aCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGa                                      >  1:204699/1‑249 (MQ=255)
                                                                                                                                                                                                               ggataataaAATGCGATTAATGCCGCCCTTAAAAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCTCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATaa                           <  1:864748/248‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaaaattGCGATTGATGCCGCCCTTAA‑AAACCCTTTTCGTCGCCTTCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAAGCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCTATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                        <  1:517094/243‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaaaaaGGCGAGTGAGGCCGCCCTTAATAAACCCTTT‑CGGCGCTTGCTGATAA‑GCCAGTACCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTGGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                        <  2:108015/247‑1 (MQ=255)
                                                                                                                                                                                                                   aataaAATGCGCTTGATGCCGCCCTTAATAAACCCTTT‑CGTCGACTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTACCGCTGGCAATGTTATCATGATCACGAAGTTCTTTGACCAATAATTTAAACGAAGACTGGGTAAgatg                       >  2:548698/1‑248 (MQ=255)
                                                                                                                                                                                                                    aaaaaaTGTGATGGATTCCCCCTTTAATAAACCCTTC‑CGTCGCCTG‑TGAAAAAGCCAGTCCCACCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGc                      <  2:584530/249‑1 (MQ=255)
                                                                                                                                                                                                                     taaAATGCGATTGATGCCGCCCTTAATAACCCCTTT‑CGTCGCCTGCGGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACg                                                                 <  1:477259/207‑1 (MQ=255)
                                                                                                                                                                                                                            cgagtGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              <  1:579770/247‑1 (MQ=255)
                                                                                                                                                                                                                            cGATTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              >  1:737599/1‑251 (MQ=255)
                                                                                                                                                                                                                             gATTGATGCCGCCCTTAATAAACCATTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              >  2:444505/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTTGATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcat            >  1:694870/1‑251 (MQ=255)
                                                                                                                                                                                                                                tcaTGCCGCCCTTAATAAACCCTTT‑CGCCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           <  1:410983/248‑1 (MQ=255)
                                                                                                                                                                                                                                 gATGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           >  2:102773/1‑249 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGa        <  1:318967/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGCCCTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGATATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           >  1:188689/1‑248 (MQ=255)
                                                                                                                                                                                                                                        cccTTAATAAACCCTTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGAAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAt  <  1:672878/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTT‑CGTCGCCTGCTGATAA‑GCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAT  >  NC_000913/2384562‑2385044

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