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→AAAyfbKIPR002035/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

ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAT  >  NC_000913/2384562‑2385044
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTTCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGGCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTtc                                                                                                                                                                                                                                            <  2:240014/250‑1 (MQ=255)
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTtcg                                                                                                                                                                                                                                           <  1:855375/251‑1 (MQ=255)
       aCTTTCTGCTCTTCAGAGATGGTAACAATGTAGCTGTAGTTGCCGTTACCAACATCTGCAATTCGCACCATCATCGCCTCGGTGTAATTGCGATTCCCCACGCCAAACGTCGACAGAGTGACACCAGTCTCCCGCTGTTTTTTGACCATTGATTCTATCTATTTTGGAGCGTCAATGCCAGCGTTAAAGTCACCGCCAGTGGCGCAAAAAATGCGATTGATGCCGCCCTGAATAATCCCTTTCTTCGCct                                                                                                                                                                                                                                     <  2:537101/250‑1 (MQ=255)
            cTGCGCTTCAGAGGGGGTATCACTGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTAGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACAATTGATTAAATCGATTTTGGATCGTCAATGCCAAAGTTAAAGTCACCGTCAGTGGCTAATATAATGCGATTGATGACGCCCTTAATAAACCCTTTCGTCGCctgctgc                                                                                                                                                                                                                                >  2:561104/1‑249 (MQ=255)
            cTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGAt                                                                                                                                                                                                                               >  1:165564/1‑251 (MQ=255)
               cgcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATa                                                                                                                                                                                                                              >  1:108015/1‑249 (MQ=255)
               cgcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATa                                                                                                                                                                                                                              >  1:637858/1‑249 (MQ=255)
                gcTTCAGAGAGGGAATCAATGTAGCTGTAGTTGCCGTTACCACCATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACCGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCTGATAAGc                                                                                                                                                                                                                          >  2:359351/1‑251 (MQ=255)
                    cAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGt                                                                                                                                                                                                                      >  1:584530/1‑251 (MQ=255)
                        gagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTcc                                                                                                                                                                                                                   >  1:853928/1‑250 (MQ=255)
                                  aTGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCAcc                                                                                                                                                                                                        <  1:435665/251‑1 (MQ=255)
                                                     tACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCg                                                                                                                                                                                     >  1:566776/1‑251 (MQ=255)
                                                       ccAACATCGGCAATTCGCACCATCATTGCCCTGTTGTAATTGCTATTCCCCTCGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTAATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGc                                                                                                                                                                                   <  1:576928/251‑1 (MQ=255)
                                                        cAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCa                                                                                                                                                                                                    >  2:546955/1‑233 (MQ=255)
                                                        cAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCa                                                                                                                                                                                                    <  1:546955/233‑1 (MQ=255)
                                                                 cAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCgc                                                                                                                                                                          >  1:32784/1‑249 (MQ=255)
                                                                           catcatTGCCTCTTTGTAATTGCTATCCCCCCCGCCAACCGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCgg                                                                                                                                                               <  1:550782/251‑1 (MQ=255)
                                                                              catTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAATCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCAt                                                                                                                                                            >  1:800542/1‑251 (MQ=255)
                                                                                       gTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCg                                                                                                                                                   >  1:398427/1‑251 (MQ=255)
                                                                                       gTTGTAATTACTATTACACACGCCAAACGTCGACAGAGTAACACAAGACTCCAGCTGTTTTTTGACCATTGATTCAATCGATTTAGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGCAGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAACCCGGCACCGCCATTGGTACTGACTTCGACATCCAGCGGATCAATTGCGGCATTCATTTCCg                                                                                                                                                   >  2:44134/1‑251 (MQ=255)
                                                                                          gTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCt                                                                                                                                                <  1:652253/251‑1 (MQ=255)
                                                                                                  tATTCCCCCCGCCAAACGTCGACAGAGTAACCCCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCAAACGTTAAAGTCCCCGTCAGTGGCTAATAAAAGGCGATTGATGCCGCCCTTAATAAACCCTTTCGGCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGa                                                                                                                                           <  1:362430/248‑1 (MQ=255)
                                                                                                  tATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTAATTTCCGCCTTATGa                                                                                                                                           >  2:362430/1‑248 (MQ=255)
                                                                                                   atcccccACGCCAAACGCCGAAAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGAGGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTc                                                                                                                                        <  1:354794/247‑1 (MQ=255)
                                                                                                     tCCCCACGCCAACCGAAGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATCCAATCGATTTTGGAGCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAAGAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCGGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTccc                                                                                                                                      <  1:300386/250‑1 (MQ=255)
                                                                                                        ccACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCCGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCgg                                                                                                                                    >  2:864748/1‑249 (MQ=255)
                                                                                                           cGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                                >  1:235725/1‑250 (MQ=255)
                                                                                                            gCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                                >  2:318967/1‑249 (MQ=255)
                                                                                                                 aCGTCGACAGAGTAACAGCAGCCTCCCGCTGGTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTAACCGTCAGTGTCTAATAAAATGCGCTTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCAATTGGTCCTGCCTTCGGCATCCAGCGCATCAATTGCGGCATTAATTTCCGCCTTATTTCTCCCgtctatatattgt                                                                                                                         >  2:166062/1‑239 (MQ=255)
                                                                                                                  cGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgc                                                                                                                         >  1:188506/1‑250 (MQ=255)
                                                                                                                    ccgaagtatgacacccaaccccccGCTGTTTTATGCCCATTGTTTCAACCGATTGTGGTTCGTAATTCCCAACTTTAATGCCACCGTCGGTGGTTAATAAAAGGCAATTGTTGCCCCCCTTAATAACCCCCTTCCACGCCTGTTGATAAGCCAGTTCTGGCCCGGCCCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaat                                                                                                                      <  2:258678/230‑1 (MQ=255)
                                                                                                                     caacagagaaaccccagcaccccGGTGTTTTTTGCCAATTGATTCAATCGATTTTGGATCTGCAATGCCAACTTCAAAGTCTCCGTCAGTGGCTAATAAAATGCGATGGAGGCCGCCCATAATAAACCCTTTCGGCGCCTGTGGATAAGCCTGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     <  2:275908/231‑1 (MQ=255)
                                                                                                                     cGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     >  1:11136/1‑251 (MQ=255)
                                                                                                                     cGACAGAGTAACACCAGACTCCCGCTGTTTTGTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCCCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatg                                                                                                                     <  1:8308/251‑1 (MQ=255)
                                                                                                                      gACAGATTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTAAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgc                                                                                                                    <  1:725015/251‑1 (MQ=255)
                                                                                                                        cAGAGTAACACCAGTCTCCCTCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTAAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTTCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgcaa                                                                                                                  <  1:456057/251‑1 (MQ=255)
                                                                                                                                             cTGTGTTTTGACCTGTTATTCAACAGATTTTGGATCGTCAATGCAAACGTAAAATCCACCGCCAGTGGCAAAAAAAAGGCGATTGATGCCCCCCTTAAAAACCCCTTTCGCCGCCTGCTGATAAGCTAGTCCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTa                                                                                               <  2:188506/249‑1 (MQ=255)
                                                                                                                                                ttttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTc                                                                                           >  2:265297/1‑250 (MQ=255)
                                                                                                                                                tttttggcccATTGATTCAATCGTTTGGGGATCGGCAATGCAAACTTTAAAGTCACCGTCAGTGGAAAATAAAATGCGATTGATGCCGCCCTTAAAAACCCCTTTCGTCTCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTGCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTc                                                                                           <  1:44134/242‑1 (MQ=255)
                                                                                                                                                       aCCATTGATTCAATCGATTTTGGATAGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGCTTGATGCCGCCCTTAATAAAACCTTTCGTAGCCTGCTGATAAGCCAGTTCCAGCCCTGCACCGCCATTGGTACTGCCTTCGGCATACAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCACGGAGATAGAAGGCAATGCAATACGGGGTTCGCCATAGTAGGTCACGATGGc                                                                                   >  2:369114/1‑251 (MQ=255)
                                                                                                                                                             gATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTAATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCACTGtt                                                                             >  1:792949/1‑251 (MQ=255)
                                                                                                                                                               ttCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTAt                                                                           >  1:183773/1‑251 (MQ=255)
                                                                                                                                                                 cAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATACAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATcc                                                                         >  1:74796/1‑251 (MQ=255)
                                                                                                                                                                            ggATCGTCAATGCCAACGTTAAAGTCACCGTCAGAGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGa                                                                >  1:671887/1‑249 (MQ=255)
                                                                                                                                                                                   aaaTGCCAACGTTAAAGCCACCGTAAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTg                                                       <  1:866657/250‑1 (MQ=255)
                                                                                                                                                                                          aaCGTTAAAGTCCCCGACAGTGGCTAAAAAAAGCCGATTGATGCCGCCTTTAACAAGCCCTTTCGTCGCCGGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTAACCTGCTCACGAAGTTCTTTGACCaata                                                <  1:87398/251‑1 (MQ=255)
                                                                                                                                                                                            cGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCaataat                                              >  1:144685/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGTCACCGTCAGAGGCTAATAAAATGAGATTGCTGCCGCCCTTAATAACCACTTTCGTCGCATGCTGACACGACAGATCCAGCCCGGCACCGCCATTGGTAATGCCTTCGGCAACCAACGACTCAATTGCGGCATAAATTTCCGCCTTATGCCTCCGAGAGATAGAAGGCAATGCAATACGGGAGTCGCCCGCGTCGGTGAAGAAGGCACAGTTTTACTGCTCACGAAGTTCTTTGCCCAATAATTTCa                                          >  2:552245/1‑250 (MQ=255)
                                                                                                                                                                                                      aCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGa                                      >  1:204699/1‑249 (MQ=255)
                                                                                                                                                                                                               ggataataaAATGCGATTAATGCCGCCCTTAAAAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCTCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATaa                           <  1:864748/248‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaaaattGCGATTGATGCCGCCCTTAAAAACCCTTTTCGTCGCCTTCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAAGCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCTATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                        <  1:517094/243‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaaaaaGGCGAGTGAGGCCGCCCTTAATAAACCCTTTCGGCGCTTGCTGATAAGCCAGTACCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTGGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                        <  2:108015/247‑1 (MQ=255)
                                                                                                                                                                                                                   aataaAATGCGCTTGATGCCGCCCTTAATAAACCCTTTCGTCGACTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTACCGCTGGCAATGTTATCATGATCACGAAGTTCTTTGACCAATAATTTAAACGAAGACTGGGTAAgatg                       >  2:548698/1‑248 (MQ=255)
                                                                                                                                                                                                                    aaaaaaTGTGATGGATTCCCCCTTTAATAAACCCTTCCGTCGCCTGTGAAAAAGCCAGTCCCACCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGc                      <  2:584530/249‑1 (MQ=255)
                                                                                                                                                                                                                     taaAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCGGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACg                                                                 <  1:477259/207‑1 (MQ=255)
                                                                                                                                                                                                                            cgagtGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              <  1:579770/247‑1 (MQ=255)
                                                                                                                                                                                                                            cGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              >  1:737599/1‑251 (MQ=255)
                                                                                                                                                                                                                             gATTGATGCCGCCCTTAATAAACCATTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc              >  2:444505/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcat            >  1:694870/1‑251 (MQ=255)
                                                                                                                                                                                                                                tcaTGCCGCCCTTAATAAACCCTTTCGCCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           <  1:410983/248‑1 (MQ=255)
                                                                                                                                                                                                                                 gATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           >  2:102773/1‑249 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGa        <  1:318967/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGATATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatc           >  1:188689/1‑248 (MQ=255)
                                                                                                                                                                                                                                        cccTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGAAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAt  <  1:672878/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAT  >  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: