Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,314,837 T→A I130F (ATC→TTC)  yciE ← DUF892 domain‑containing protein YciE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,8370TA98.3% 173.4 / NA 59I130F (ATC→TTC) yciEDUF892 domain‑containing protein YciE
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (31/26);  minor base G (0/1);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

ACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCAGTTGAACA  >  NC_000913/1314598‑1315066
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
acacGATAGCGTACGCCCTGATCCGTGGGTGTTTATTTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGTGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCTGAATATTCTGAATCATCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTg                                                                                                                                                                                                                            <  2:740695/251‑1 (MQ=255)
    gATAGCGTACGCCTCGATCCATGGGTGTTTATGTTCTCAGCAGGATATAGTTTTTAATTCCGCACAAAAAAGACCCCGTAATA‑TACGGGGACAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCTAGTCGACCATTTGCTTTTCCTCATTTAAAATCACTTCGAAGGTTGGAAttta                                                                                                                                                                                                                           >  2:587176/1‑246 (MQ=255)
     aTAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGGTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCGTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCt                                                                                                                                                                                                                       <  2:432837/251‑1 (MQ=255)
       aGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTg                                                                                                                                                                                                                      <  2:94211/250‑1 (MQ=255)
          gTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTACTTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGACGCTGtctc                                                                                                                                                                                                                  >  2:278496/1‑248 (MQ=255)
          gTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATc                                                                                                                                                                                                                  >  2:572838/1‑251 (MQ=255)
            aagcatcgagacgtggggtgtgaggTTCTCAGCAGGATAGCGTTTGTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTACGCTTCTACGCCGTCAGTTTCAGAGCGAATTAAAAATTTGTCAGTTGTTTGCGGAATAGTCTTAATCAGCCAGTTGGCCATTTGCTGTTCCTCATTTAAAATCGCTTCAAAGGTGGGAAGTGAAGCGGTATCAc                                                                                                                                                                                                                <  2:895041/227‑1 (MQ=255)
             cGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGGCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCAcc                                                                                                                                                                                                               <  1:55579/251‑1 (MQ=255)
                  cGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAt                                                                                                                                                                                                          >  2:371020/1‑251 (MQ=255)
                  cGATCCGTGGGTGTTTATGTTCTCAGCAGGATAGCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAAGCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCa                                                                                                                                                                                                           <  2:157499/250‑1 (MQ=255)
                   gATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGAATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAtt                                                                                                                                                                                                         >  2:818791/1‑251 (MQ=255)
                    aTCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAttt                                                                                                                                                                                                        >  1:353937/1‑251 (MQ=255)
                          gggcgtTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATCAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGACTTAAAAATTTCTCAGTTGTTTGCGGCATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTtgct                                                                                                                                                                                                  >  1:909362/5‑251 (MQ=255)
                               ttgaTGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTGCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATGGAAGCTGTATCACCGGCATGTTTtgctgctg                                                                                                                                                                                              <  2:127281/247‑1 (MQ=255)
                                    gTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGAAATGGTAGAAAGGGTATTATTGCTTCGCGTCTACGCCATCAGTTTCAGAGAGAATTAAAAATTTCTCAGTTGTTGGCGGAATATTCTGAATAAGCCAGGCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAAGTGAAGCTGTATCACCGGCATTTGTTGCTGCTGCTAAt                                                                                                                                                                                         <  2:769867/250‑1 (MQ=255)
                                      tctcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAg                                                                                                                                                                                       <  1:818791/250‑1 (MQ=255)
                                       cgcagccagatttCGTTTTTAATTTCGCACAAAAACGATCTCGTAATAGTA‑GGGGTCATTAAGGACATTGTATTAAGTGTATTATGTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGGTTGCGTAATATTCTTAATCAGCCAGTCGGCTATTTGCTTTTCCGCATTTAATATCGCTTCGAGGGTTTGAATTGAAGCTGTATCACCGGCATTTTGTGCTGCTGCTAATAGTg                                                                                                                                                                                     <  2:415895/239‑1 (MQ=255)
                                           gcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                                                 >  1:838506/1‑251 (MQ=255)
                                                  atCGTTTTTAATTCCGCAAAAAAACGTCCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGGATTATTTCTTCGCTTCTACGACAGCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCTGAATATTCTGAATCAGCCAGTCGGCTATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAg                                                                                                                                                                          <  2:86497/251‑1 (MQ=255)
                                                    cGTTTTTACTTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGTTATAAAGGGTATTATTTCTTAGTTTCTACGCCATCCGATTAAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTACTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACAGGCATTTTTTGCTGTTGCTAATAGTGAgtggtacaaggc                                                                                                                                                                        >  2:759937/1‑240 (MQ=255)
                                                     gTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TAAGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGATTCTACGCCATCGGTTTCAGAGCGAATTAAAAATTTCTCAGTAGTTTGCGGCATATTCTGAATCAGCCAGTGGGCAATTTGCTTTTCCTCATTTATAATCGCTTCGAAGGTTGGAATTGAAACTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCg                                                                                                                                                                       >  1:885895/1‑251 (MQ=255)
                                                        tttAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAt                                                                                                                                                                     <  1:93675/250‑1 (MQ=255)
                                                        tttAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAtt                                                                                                                                                                    >  1:257806/1‑251 (MQ=255)
                                                           aaTTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTc                                                                                                                                                                  >  1:538638/1‑250 (MQ=255)
                                                                         aCGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGAc                                                                                                                                                   >  1:228875/1‑251 (MQ=255)
                                                                          cgacagcGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGAAGTATTTCTTCGCTTCTACGCCAGCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGACTATTCTGAATCTGCCAGTCGACAATTTGCTTTTCATCATTTAAAATCTCTTCGAAGGTTGGACTTGTAGCTGTATCACCGGCATTTTTTGCTGCTGATACTAGTGAGGTGTAACTGGAGATTTCAAATTGCTCGAAgccc                                                                                                                                                  >  2:179354/7‑248 (MQ=255)
                                                                               ccGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCTCCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATcc                                                                                                                                             <  1:466633/251‑1 (MQ=255)
                                                                                   aaTA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCt                                                                                                                                          >  1:139738/1‑250 (MQ=255)
                                                                                   aaTA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTa                                                                                                                                         <  1:129189/251‑1 (MQ=255)
                                                                                   aaTA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATAGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGcc                                                                                                                                          >  1:105395/1‑249 (MQ=255)
                                                                                          cGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATag                                                                                                                                     >  1:602874/1‑249 (MQ=255)
                                                                                           ggggTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGc                                                                                                                                  >  2:904813/1‑251 (MQ=255)
                                                                                                cAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTg                                                                                                                             >  2:26370/1‑251 (MQ=255)
                                                                                                cAATAAGGAAATGGTATAAAGGTTATTATTTCTGCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAACCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGATAATAGTGAGGTGTATCAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCttt                                                                                                                              <  2:533012/250‑1 (MQ=255)
                                                                                                         cATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTtca                                                                                                                    <  2:344815/251‑1 (MQ=255)
                                                                                                         cATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTtc                                                                                                                     >  1:527841/1‑250 (MQ=255)
                                                                                                                        attattTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAAtata                                                                                                     <  1:972820/251‑1 (MQ=255)
                                                                                                                           attTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATacca                                                                                                  <  1:570388/251‑1 (MQ=255)
                                                                                                                                    cTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGACTATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGTAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACAGATTGAc                                                                                         >  2:421932/1‑251 (MQ=255)
                                                                                                                                       cTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGc                                                                                      >  2:643022/1‑251 (MQ=255)
                                                                                                                                              aTCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCAAGTCGGCCATTTGCTTTTCCTCATGTAAAATCGCTTCGAAGGGTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCaatcg                                                                               <  2:826883/251‑4 (MQ=255)
                                                                                                                                                        gagCGAATTCAAAATTTCTCCGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGATAATAGTGAGGTGTAACAGGCGCTTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCAtttt                                                                     >  2:469090/1‑251 (MQ=255)
                                                                                                                                                                   aaaTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGACGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGATTATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAAt                                                          >  2:203772/1‑251 (MQ=255)
                                                                                                                                                                     aTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCt                                                        >  1:451302/1‑251 (MQ=255)
                                                                                                                                                                        tctcAGTTGTTTGCGGAATATTATGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGCATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGCCATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCttt                                                      >  2:210148/1‑250 (MQ=255)
                                                                                                                                                                           cAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTaa                                                    <  1:371020/249‑1 (MQ=255)
                                                                                                                                                                               tgttTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACt                                               <  1:12400/250‑1 (MQ=255)
                                                                                                                                                                                   tGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCGTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGTTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCGTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCGTTAATGACTGAAc                                           <  2:909362/250‑1 (MQ=255)
                                                                                                                                                                                       gAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGaa                                      <  1:572838/251‑1 (MQ=255)
                                                                                                                                                                                                gAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTa                             <  1:241536/251‑1 (MQ=255)
                                                                                                                                                                                                 aaTCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTAc                            >  1:953493/1‑251 (MQ=255)
                                                                                                                                                                                                 aaTCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTAc                            <  2:353937/251‑1 (MQ=255)
                                                                                                                                                                                                          gTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGATGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCACAAGCGCAGCCATTTTACTCATGGACTCTTTAATGACTGAACGTGACATGTCATTACGATCAAGa                    >  1:908312/1‑250 (MQ=255)
                                                                                                                                                                                                          gTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAAGTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGaa                   <  1:643022/251‑1 (MQ=255)
                                                                                                                                                                                                             ggCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAAt                  >  1:953143/1‑249 (MQ=255)
                                                                                                                                                                                                              gCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGGTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGAAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGTATAGt               >  1:310804/1‑251 (MQ=255)
                                                                                                                                                                                                                    tGCTTTTCCTCATTTAACATCGCTTCGAAGGTTGTAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCa          <  1:243009/250‑1 (MQ=255)
                                                                                                                                                                                                                         atCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGTTCAAGAATAGTTTCCAGTTGa     >  2:706860/2‑250 (MQ=255)
                                                                                                                                                                                                                           ccTCAGTTAAAATCGCTTCGAAGGTTGGAAGTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATGTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCAGTTGAACa  <  2:607667/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
ACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCAGTTGAACA  >  NC_000913/1314598‑1315066

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