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

ACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCAGTTGAACA  >  NC_000913/1314598‑1315066
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
acacGATAGCGTACGCCCTGATCCGTGGGTGTTTATTTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGTGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCTGAATATTCTGAATCATCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTg                                                                                                                                                                                                                            <  2:740695/251‑1 (MQ=255)
    gATAGCGTACGCCTCGATCCATGGGTGTTTATGTTCTCAGCAGGATATAGTTTTTAATTCCGCACAAAAAAGACCCCGTAATATACGGGGACAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCTAGTCGACCATTTGCTTTTCCTCATTTAAAATCACTTCGAAGGTTGGAAttta                                                                                                                                                                                                                           >  2:587176/1‑246 (MQ=255)
     aTAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGGTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCGTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCt                                                                                                                                                                                                                       <  2:432837/251‑1 (MQ=255)
       aGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTg                                                                                                                                                                                                                      <  2:94211/250‑1 (MQ=255)
          gTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTACTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGACGCTGtctc                                                                                                                                                                                                                  >  2:278496/1‑248 (MQ=255)
          gTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATc                                                                                                                                                                                                                  >  2:572838/1‑251 (MQ=255)
            aagcatcgagacgtggggtgtgaggTTCTCAGCAGGATAGCGTTTGTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTACGCTTCTACGCCGTCAGTTTCAGAGCGAATTAAAAATTTGTCAGTTGTTTGCGGAATAGTCTTAATCAGCCAGTTGGCCATTTGCTGTTCCTCATTTAAAATCGCTTCAAAGGTGGGAAGTGAAGCGGTATCAc                                                                                                                                                                                                                <  2:895041/227‑1 (MQ=255)
             cGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGGCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCAcc                                                                                                                                                                                                               <  1:55579/251‑1 (MQ=255)
                  cGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAt                                                                                                                                                                                                          >  2:371020/1‑251 (MQ=255)
                  cGATCCGTGGGTGTTTATGTTCTCAGCAGGATAGCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAAGCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCa                                                                                                                                                                                                           <  2:157499/250‑1 (MQ=255)
                   gATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGAATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAtt                                                                                                                                                                                                         >  2:818791/1‑251 (MQ=255)
                    aTCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCAttt                                                                                                                                                                                                        >  1:353937/1‑251 (MQ=255)
                          gggcgtTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATCAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGACTTAAAAATTTCTCAGTTGTTTGCGGCATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTtgct                                                                                                                                                                                                  >  1:909362/5‑251 (MQ=255)
                               ttgaTGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTGCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATGGAAGCTGTATCACCGGCATGTTTtgctgctg                                                                                                                                                                                              <  2:127281/247‑1 (MQ=255)
                                    gTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGAAATGGTAGAAAGGGTATTATTGCTTCGCGTCTACGCCATCAGTTTCAGAGAGAATTAAAAATTTCTCAGTTGTTGGCGGAATATTCTGAATAAGCCAGGCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAAGTGAAGCTGTATCACCGGCATTTGTTGCTGCTGCTAAt                                                                                                                                                                                         <  2:769867/250‑1 (MQ=255)
                                      tctcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAg                                                                                                                                                                                       <  1:818791/250‑1 (MQ=255)
                                       cgcagccagatttCGTTTTTAATTTCGCACAAAAACGATCTCGTAATAGTAGGGGTCATTAAGGACATTGTATTAAGTGTATTATGTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGGTTGCGTAATATTCTTAATCAGCCAGTCGGCTATTTGCTTTTCCGCATTTAATATCGCTTCGAGGGTTTGAATTGAAGCTGTATCACCGGCATTTTGTGCTGCTGCTAATAGTg                                                                                                                                                                                     <  2:415895/239‑1 (MQ=255)
                                           gcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                                                 >  1:838506/1‑251 (MQ=255)
                                                  atCGTTTTTAATTCCGCAAAAAAACGTCCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGGATTATTTCTTCGCTTCTACGACAGCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCTGAATATTCTGAATCAGCCAGTCGGCTATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAg                                                                                                                                                                          <  2:86497/251‑1 (MQ=255)
                                                    cGTTTTTACTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGTTATAAAGGGTATTATTTCTTAGTTTCTACGCCATCCGATTAAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTACTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACAGGCATTTTTTGCTGTTGCTAATAGTGAgtggtacaaggc                                                                                                                                                                        >  2:759937/1‑240 (MQ=255)
                                                     gTTTTTAATTCCGCACAAAAACGACCCCGTAATATAAGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGATTCTACGCCATCGGTTTCAGAGCGAATTAAAAATTTCTCAGTAGTTTGCGGCATATTCTGAATCAGCCAGTGGGCAATTTGCTTTTCCTCATTTATAATCGCTTCGAAGGTTGGAATTGAAACTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCg                                                                                                                                                                       >  1:885895/1‑251 (MQ=255)
                                                        tttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAt                                                                                                                                                                     <  1:93675/250‑1 (MQ=255)
                                                        tttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAtt                                                                                                                                                                    >  1:257806/1‑251 (MQ=255)
                                                           aaTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTc                                                                                                                                                                  >  1:538638/1‑250 (MQ=255)
                                                                         aCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGAc                                                                                                                                                   >  1:228875/1‑251 (MQ=255)
                                                                          cgacagcGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGAAGTATTTCTTCGCTTCTACGCCAGCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGACTATTCTGAATCTGCCAGTCGACAATTTGCTTTTCATCATTTAAAATCTCTTCGAAGGTTGGACTTGTAGCTGTATCACCGGCATTTTTTGCTGCTGATACTAGTGAGGTGTAACTGGAGATTTCAAATTGCTCGAAgccc                                                                                                                                                  >  2:179354/7‑248 (MQ=255)
                                                                               ccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCTCCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATcc                                                                                                                                             <  1:466633/251‑1 (MQ=255)
                                                                                   aaTATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCt                                                                                                                                          >  1:139738/1‑250 (MQ=255)
                                                                                   aaTATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTa                                                                                                                                         <  1:129189/251‑1 (MQ=255)
                                                                                   aaTATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATAGCTTCGAAGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGcc                                                                                                                                          >  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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
ACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATGTCATTACGATCAAGAATAGTTTCCAGTTGAACA  >  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: