Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,175,991 T→A K111* (AAG→TAG)  gatZ ← putative tagatose‑1,6‑bisphosphate aldolase 2 chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,175,9910TA97.2% 211.3 / ‑1.4 71K111* (AAG→TAG) gatZputative tagatose‑1,6‑bisphosphate aldolase 2 chaperone
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (34/35);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

GCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCG  >  NC_000913/2175744‑2176233
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gcATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATAc                                                                                                                                                                                                                                                 >  1:525085/1‑251 (MQ=255)
 cATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCCGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATAcc                                                                                                                                                                                                                                                >  2:226514/1‑251 (MQ=255)
  atattCAACATGGGTGATGTGTCCTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCCCACTCGGGCAATCTGTCGCCACTCGTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                               <  2:669138/248‑1 (MQ=255)
  aTCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATAcc                                                                                                                                                                                                                                                >  1:461618/1‑250 (MQ=255)
      tCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                               <  1:803417/247‑1 (MQ=255)
      tCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCAAGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTc                                                                                                                                                                                                                                           >  1:762292/1‑251 (MQ=255)
         acaggggTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGAc                                                                                                                                                                                                                                        <  2:777872/247‑1 (MQ=255)
                tGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGAtttt                                                                                                                                                                                                                                 >  1:785828/1‑251 (MQ=255)
                 gATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGAttttt                                                                                                                                                                                                                                <  2:133441/251‑1 (MQ=255)
                          tgactgaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATcgccg                                                                                                                                                                                                                       >  1:230928/1‑251 (MQ=255)
                           gtccgaagtgcgcatgtctccacGCCCGGAACCGGAACTTCGGTGCCAAAGACATAGCGCAGTTGCTGACTCTGGCAAGCTTGCGCCACACTTACCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATAGTCCTATACCAGCTCGACGGCTTTGTCCATcgcc                                                                                                                                                                                                                        <  2:362434/227‑1 (MQ=255)
                           gactgaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATcgcc                                                                                                                                                                                                                        <  1:793411/249‑1 (MQ=255)
                             ctgaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCAACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGccgcat                                                                                                                                                                                                                    >  1:415189/1‑251 (MQ=255)
                                     cgcTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTAGACGGCTTTTTCCATCGCCGCATCCGCAttt                                                                                                                                                                                                            >  2:653977/1‑251 (MQ=255)
                                         ggCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTc                                                                                                                                                                                                          <  1:492991/249‑1 (MQ=255)
                                            cTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGCCGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTtgct                                                                                                                                                                                                     >  2:330309/1‑251 (MQ=255)
                                                ccGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGACGCATCCGCATTTTCTTGCTGCca                                                                                                                                                                                                 <  1:550821/251‑1 (MQ=255)
                                                          ccggaacTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGc                                                                                                                                                                                       >  1:266601/1‑251 (MQ=255)
                                                                     gTGCCAATGCCATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATc                                                                                                                                                                            <  2:572410/251‑1 (MQ=255)
                                                                     gTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATc                                                                                                                                                                            >  1:428822/1‑251 (MQ=255)
                                                                     gTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCAGCCGCATTTTCTTGCTGCCAGCAGTTTGGACCCAGATGATc                                                                                                                                                                            <  1:736749/251‑1 (MQ=255)
                                                                        ccAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGAGCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATcgcc                                                                                                                                                                         <  2:266601/251‑1 (MQ=255)
                                                                          aaTGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATcgccgc                                                                                                                                                                       >  1:167870/1‑251 (MQ=255)
                                                                                 tAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGaat                                                                                                                                                                >  1:299471/1‑251 (MQ=255)
                                                                                      cAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGc                                                                                                                                                           <  1:270197/251‑1 (MQ=255)
                                                                                         ttGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGtt                                                                                                                                                        >  1:566885/1‑251 (MQ=255)
                                                                                          tGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGtt                                                                                                                                                        <  1:573091/250‑1 (MQ=255)
                                                                                          tGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGtt                                                                                                                                                        <  1:167166/250‑1 (MQ=255)
                                                                                            cTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTcgc                                                                                                                                                     >  1:393782/1‑251 (MQ=255)
                                                                                             tCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTcgcg                                                                                                                                                    >  1:817185/1‑251 (MQ=255)
                                                                                              cACGCTGGCAAACTGGCGCCACACCTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCTTTTCTGGTGCGAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGTCAGCAGTTTGGCCCCCGATGATCGCCGCCGAGAATAATGCGTTcgc                                                                                                                                                     <  2:244159/249‑1 (MQ=255)
                                                                                               aCGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTg                                                                                                                                                  <  2:230928/251‑1 (MQ=255)
                                                                                                     gtaatttttCGCCACATGTTCCGCAGTCAAGCAAAGCCCAGCAGCTCCTTTCGTAAACGTTGTTGGTGCTTATGGGATTGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACGGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTGTCCATCGCCGCATCCGCATGTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCaaa                                                                                                                                              <  2:512554/241‑1 (MQ=255)
                                                                                                            gTCGCCCCACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTCCCGCAACCGTGTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATGTTTCCATCGCCGCATCCTCATTTTCTTGCTGCCAGCAGTTTGGCTCCAGATGATCGCCGCCGAGAATAATGTGTTCGCGTGCaaatcccac                                                                                                                                        <  2:560938/248‑7 (MQ=255)
                                                                                                                ccACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTCACGGTATAGGATCCCCCGCGCAGGACATTGACGCATCAAGATGCATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCGTCCGCATATTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGTGTGAAAACCCAACTTTATcc                                                                                                                                 >  2:755894/1‑250 (MQ=255)
                                                                                                                  acacTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGc                                                                                                                               >  2:694162/1‑251 (MQ=255)
                                                                                                                   cacTTTCCGCAGCAAAGCAAAGCACAGCATCTCGTTCCGCAACCGTTTCTTGTGCTAACGGTATGGGAACCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGTATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCa                                                                                                                              <  1:328234/251‑1 (MQ=255)
                                                                                                                       ttCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCg                                                                                                                          >  1:142314/1‑251 (MQ=255)
                                                                                                                               caaagcaaagCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACa                                                                                                                  <  1:155518/251‑1 (MQ=255)
                                                                                                                                aaagcaaagcACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGCTGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACa                                                                                                                  >  2:15014/1‑250 (MQ=255)
                                                                                                                                       atcacaGCAGCTCGTTCCGCACCCGTTTCTTGGGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATTCCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTcgcg                                                                                                          <  2:461618/249‑1 (MQ=255)
                                                                                                                                       agcaCAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTATCGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTAACGGCAATCGTAAAAACAAATTcgcg                                                                                                          <  2:773765/251‑1 (MQ=255)
                                                                                                                                            agcagcTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGt                                                                                                     <  1:158388/251‑1 (MQ=255)
                                                                                                                                               agcTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTg                                                                                                  <  1:599371/251‑1 (MQ=255)
                                                                                                                                                    gTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGCCGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCCCTCGTAAAAACAACTTCGCGAAAATCTGCCggg                                                                                             >  2:858517/1‑250 (MQ=255)
                                                                                                                                                    gTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGgt                                                                                             <  1:192138/251‑1 (MQ=255)
                                                                                                                                                    gTTCCGCAACCGTTGCTGTTACTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTGCTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCACATCCGCCTTTTCTTGCCGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGTGCATAATTAGTTCGCGTGCAAACACAACATTATCGGCAGTCGTAAACACAAATGCGCTAAAGACTGCAGgt                                                                                             >  2:485089/1‑251 (MQ=255)
                                                                                                                                                         gCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCAt                                                                                        >  1:220675/1‑251 (MQ=255)
                                                                                                                                                            aCCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTcc                                                                                     <  1:226514/251‑1 (MQ=255)
                                                                                                                                                                            aCGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAat                                                                     <  1:653977/251‑1 (MQ=255)
                                                                                                                                                                                  tGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCACCTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCACattgatt                                                                >  1:540156/1‑250 (MQ=255)
                                                                                                                                                                                   gggAACCCCCGCGCAGGACATGAACGCAGCAAGATGAATTTTACTGAAGCCGGCCCGAACATATTCCTATCCTCGCTCGACGGATTTTTCCATCGCCGCAGCCGCATTTTCTTGCTGCCAGCAGTGTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCATCTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAattgagtg                                                               <  2:63831/250‑4 (MQ=255)
                                                                                                                                                                                      cccccccGCGCAGGTGAGTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTACTCCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGTTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAACCCCATCTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCGGCCGTTGTCATTCCGGTATAACCGCCAAATTGATTGACCt                                                           <  2:334886/249‑1 (MQ=255)
                                                                                                                                                                                       tCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGTACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGGGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAATACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCATATTGATTGACCTg                                                          >  1:429313/1‑251 (MQ=255)
                                                                                                                                                                                        cccccGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGGTTCATTCCGGGATAACCGCCAATTTTATTGACCTgg                                                         >  2:685057/1‑251 (MQ=255)
                                                                                                                                                                                          cccGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGtt                                                       <  1:489613/251‑1 (MQ=255)
                                                                                                                                                                                                cAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATATCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATAGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCCACTTTATCGGCAATCGTAAAAACCAATTCGCGAAAGTCTGCCGGTGTACTTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGt                                                 >  1:571937/1‑251 (MQ=255)
                                                                                                                                                                                                 gtgACAGTGACGCAACAAGATTAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTCCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCATTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGCCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGtt                                                <  2:815623/249‑1 (MQ=255)
                                                                                                                                                                                                  gtaCATTGACTCATCAAGATGAATTTGACTGAAGCCGGTCCGAACATATCCCTAAACCAGCTCGCCGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCTATCGTAAAAACAAATGCGCGAAAGTCTGCCGGTTTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTg                                               <  1:292965/249‑1 (MQ=255)
                                                                                                                                                                                                     cATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCtt                                            <  1:534354/251‑1 (MQ=255)
                                                                                                                                                                                                      aTTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAGACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTc                                           >  1:334104/1‑251 (MQ=255)
                                                                                                                                                                                                      aTTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTc                                           >  2:140374/1‑251 (MQ=255)
                                                                                                                                                                                                        tGACGCCTCAAGATGAATTTTACTGAAGCCTGCACGAACATATTCCTATACCAGCTCAACTGATTTTTCCATCCCCGCATCCGCATTTTCTTGCTGCCAGCCGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCTTGCAAACCCACCTTTATCGGCAATCGTAAAACCAAATTCGAGCATGTCTGCCGGTGTCATTCCGGTTTAACCCCCAAATTTTTAGACCTGTTTTGCCGTTACTTcgt                                         >  2:843522/1‑249 (MQ=255)
                                                                                                                                                                                                           cGCATCAGGATGCATTATACTGAAGCCGGCACGCACATACTCCTATACCAGCTCGACGGATTTTAACATCGCCGCCTCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGACGTGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTCAAAACAACTTCGAGCAAGTCTGCCCGTGTCATTCAGGTCTAACCGCACAATTGTTTGACCTGGATTGACGTTGcaccaata                                       >  2:872259/1‑243 (MQ=255)
                                                                                                                                                                                                                       aaTTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGcgcg                          >  1:795644/1‑251 (MQ=255)
                                                                                                                                                                                                                            tACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTg                     >  1:196784/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTg                     <  1:15014/250‑1 (MQ=255)
                                                                                                                                                                                                                               tGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTg                  >  1:645168/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTATTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATc             >  2:277048/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cACGAACATATTCCTATACCAGCTCGACGTATTTTCCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCTGTTTGGCCACATATGATCGCCGCCGAGAATAATGCGGTCGCTTGCAAACCCAACTTTATCGGCAATCGTAAAATCAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCTCTTTGCGCGTGCTGTTGCGATCAAAt         <  1:381935/251‑1 (MQ=255)
                                                                                                                                                                                                                                               atatTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAACCCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCCGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAgcg  >  2:378742/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCG  >  NC_000913/2175744‑2176233

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