Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,252 C→A intergenic (‑9/+3) sgcC ← / ← sgcB putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,2520CA94.4% 157.9 / ‑0.7 54intergenic (‑9/+3)sgcC/sgcBputative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (29/22);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530014‑4530498
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
aCGTAAAAGACATGAGGGGCTGGGCAGAAACGCACAACCATCACTTGGATAGGGGGACTGGGGTTGTCGATCGCCATGACGACCACCAGTCCCATACCGACAAAACCAATCCCACCCTTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCAAATCACGATAATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGtt                                                                                                                                                                                                                                            <  2:177512/251‑1 (MQ=255)
   tccaaGACATGCAGGGCTAGGCCGAAACCCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCAGGACGATCACAAGTACCATACCGACAAAACCAATCACAACCGTTACCACTGCCTTAACCGCCTGCAGCCAGGAGAGGCGGAAGAC‑AAGGCCAATCAAGATCATGATGATAGGGACGAAACCGGTACCGCCGGGAGACCGGATGTAATAAAACATAATGCCCACTTTTATTGAGTTAAc                                                                                                                                                                                                                                         >  2:278491/4‑251 (MQ=255)
          cATGCGGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGGGGCTGTCGATCGCCATGACGATAACCAGCCCCATACCGGCAAAACCAATCCCAACCGTTCCACCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGCAGAT‑AAGGCCAATCACGATAATGAAGATAGGGACGAAAACGGTACCGCCGAGAGACAGAGTGTAATCAAACATAATGACTCATTTTTTTGAGATAACagggact                                                                                                                                                                                                                                  >  2:304728/1‑246 (MQ=255)
                            aaccccaaatcaactatttGATGGGGGGACTGAGGCTTGCGACCGCCAGGACGATCACCAGCCCCATCCCGAAAAAACCAACCCCAACCGTTACCCCGCCCTTAACCGCCTGCAGCCAGGGGATGCGGAAAAA‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACTGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAAcg                                                                                                                                                                                                                <  2:126717/235‑1 (MQ=255)
                             aCGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGGGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAAcgc                                                                                                                                                                                                               >  2:33998/1‑251 (MQ=255)
                                      catcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTACTCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGCTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTAATTTTACTGAGTGAACAGAGCCTTGATTAGTTGTTTTACCGCGTCGTCGtt                                                                                                                                                                                                      >  2:403405/1‑251 (MQ=255)
                                              tGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCtg                                                                                                                                                                                              >  2:358619/1‑251 (MQ=255)
                                              ccgtgggctctcccgatctggcgaccGCCATGACGACCACCAGCCCAATCCCGACAAACCCAATCCCACCCTTTACCCCGGCATAAATCGCCTGCAGCCAGGGGAGGCGGAATAT‑AAGGCCAATCACTATAATGAGGATAGGGACGAAAACGGAACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCtg                                                                                                                                                                                              <  1:447434/226‑1 (MQ=255)
                                                aTAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTg                                                                                                                                                                                                                                >  2:160746/1‑215 (MQ=255)
                                                aTAGGGGGACGGAGGCTGTCGATCGCCATGACGATCCCCAGTCCCATCCCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAACCACGATCATGATGATAGGGACGAAAAGGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTg                                                                                                                                                                                                                                <  1:160746/215‑1 (MQ=255)
                                                    ccaaccagacgctcttccgatctaTGACGATCACCAGTCCCAAACCGACAAAACCAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCGGGGGATGGGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCa                                                                                                                                                                                        <  2:225056/228‑1 (MQ=255)
                                                             ggCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGgcgc                                                                                                                                                                               >  1:56432/1‑251 (MQ=255)
                                                                 gacGATCGCCATGAAGATCACAAGTCCCATACCGACAAAACCAATCACAACCGTTACACCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGCT‑AAGGCCAATCACGATCATGCTGCTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACCTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTACCGAGTCGTCGTTGATccatg                                                                                                                                                                                              >  2:447434/3‑229 (MQ=255)
                                                                       cGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGt                                                                                                                                                                     >  1:415353/1‑251 (MQ=255)
                                                                           aTGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCa                                                                                                                                                                                        >  1:225056/1‑228 (MQ=255)
                                                                                 acccccAGTCCCATACCGACAAAACCAACCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGgcggcg                                                                                                                                                                                <  1:308791/226‑1 (MQ=255)
                                                                                 aTCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGgcggcg                                                                                                                                                                                >  2:308791/1‑230 (MQ=255)
                                                                                             aTACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTCAGCGTAGGGATGCCGTATTCACTGtt                                                                                                                                                >  2:327080/1‑250 (MQ=255)
                                                                                             aTACCGACAAAACCAACCCCAACCGTTCCCCCTGCCTTAACCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGtt                                                                                                                                                <  1:327080/250‑1 (MQ=255)
                                                                                                 cGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGACATCACTGTTGGTTc                                                                                                                                           >  2:269127/1‑251 (MQ=255)
                                                                                                 cGACAAAACCAACCCCAACCTTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGACCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTc                                                                                                                                           <  1:303298/251‑1 (MQ=255)
                                                                                                      aacccaacccaacccgtccccctcgCCTTAACCCCCGGCAGCCGGGGGATGCGGAAGAA‑AAGGCAAATCACGATCATGATGATAGGGCAGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATg                                                                                                                                      <  2:79284/227‑1 (MQ=255)
                                                                                                       aaCCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGCTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTCTTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGACGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGTTGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                     >  2:240233/1‑251 (MQ=255)
                                                                                                                cAACCGTTACCCCTGCCTTATTCGCCTGCAGCAAGGGGTGGCGGAAGAT‑AAGCCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACaa                                                                                                                            <  1:110313/251‑1 (MQ=255)
                                                                                                                        cccccTGCCTTAATCGCCGGCGGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTcc                                                                                                                    <  1:269127/250‑1 (MQ=255)
                                                                                                                         ccccTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTcc                                                                                                                    >  2:238689/1‑250 (MQ=255)
                                                                                                                              gacTTAACCGCCTGCAGCAAGGGGAGGCGGAAGAT‑AAGGAAAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGGTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGTCGGAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATAcc                                                                                                               <  1:122539/248‑1 (MQ=255)
                                                                                                                                  taaccGCCTCCAGCCAGGGGATGAGGAATAT‑AAGGCCAAACACGAATAAGATGAAAGGACCTAAATCGGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTa                                                                                                           <  2:389811/246‑1 (MQ=255)
                                                                                                                                       gCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGtt                                                                                                      >  1:373275/1‑250 (MQ=255)
                                                                                                                                         cTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAg                                                                                                                                                                                                                                              >  2:58943/1‑112 (MQ=255)
                                                                                                                                         cTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCAGGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAg                                                                                                                                                                                                                                              <  1:58943/112‑1 (MQ=255)
                                                                                                                                            cAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTa                                                                                                                                                                    >  1:7522/1‑183 (MQ=255)
                                                                                                                                                 aGGGGATGCG‑AAGATTAAGGCCATGGACGAGCATGATGATGGGAACGAAAACGGGACCGCCGAGAGACAGTATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTGTTAACGCGTCGGCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGAGAtct                                                                                           <  1:57657/251‑1 (MQ=255)
                                                                                                                                                       tGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTc                                                                                      >  1:295579/1‑250 (MQ=255)
                                                                                                                                                        ggggaaagcaaaGGCCAACAACGATAATGATGATAGGTACAAAAACTGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGGTAACAGAGCCTTGATGTGTTGTTTTAACGCGTCGTCGGTGACCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGGAATCACTGTTGGTTCGCATGGAGGTTACAATAAGTTCCATACCGTTACAGTTAAGAGGTATCTCATTCa                                                                                     <  2:331895/241‑1 (MQ=255)
                                                                                                                                                            aaGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAc                                                                                  >  1:180935/1‑249 (MQ=255)
                                                                                                                                                                   aGGCCAATCACGATCATGAGAATAGGGAAGAAAACGGTACCGCCGATAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATGTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTgg                                                                          <  1:400787/251‑1 (MQ=255)
                                                                                                                                                                     gCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGAAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGc                                                                        >  1:348703/1‑251 (MQ=255)
                                                                                                                                                                          tCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGt                                                                     >  1:231077/1‑249 (MQ=255)
                                                                                                                                                                            aCGAACATGATGAGAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTg                                                                 <  2:56432/251‑1 (MQ=255)
                                                                                                                                                                                  atgatgATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGa                                                            >  1:429585/1‑250 (MQ=255)
                                                                                                                                                                                   tgatgaTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGTGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACATTAAGGTCCATACCGTTAAAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGaa                                                           >  1:235016/1‑250 (MQ=255)
                                                                                                                                                                                     atgatAGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATg                                                        <  1:408389/251‑1 (MQ=255)
                                                                                                                                                                                          ggggACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTg                                                   <  1:238689/250‑1 (MQ=255)
                                                                                                                                                                                               cGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTTGCTGAAATGCCCTGTtcag                                              <  2:98071/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGa                                        >  1:122937/1‑250 (MQ=255)
                                                                                                                                                                                                            gccgAGAGACAGGATGTAATCAAACATAATGCCGGAGTTTATTGAGTTAACAGAGCCTTGAGTTTTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTTGGGATGCCGTAATCACTGTTGGTTCGCAGGGAGGTGACAATAAGGTCCAGACCGTTACAGTTAAGAGGGATCTCAGTCAGACAGCATTGGGGGGTAGTGGCTGAAATGCCCTGTGCAGTCAGGAACTCTTg                                 <  2:38881/251‑1 (MQ=255)
                                                                                                                                                                                                               gagagaCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGc                                >  2:63444/1‑249 (MQ=255)
                                                                                                                                                                                                                agagaCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGc                             >  2:106954/1‑251 (MQ=255)
                                                                                                                                                                                                                      aGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgttt                       >  1:354269/1‑249 (MQ=255)
                                                                                                                                                                                                                       ggATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTCAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAACTGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtgtg                      >  1:303822/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaTCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCa                <  1:270516/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCCGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGa         >  1:315908/1‑250 (MQ=255)
                                                                                                                                                                                                                                         tGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCg    >  2:119817/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cctcatTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAc  >  2:243263/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAT‑AAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530014‑4530498

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