Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,533 A→G 11.1% L301P (CTG→CCG)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,5330AG11.1% 189.9 / 12.3 63L301P (CTG→CCG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base G (4/3);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCG  >  NC_000913/2967247‑2967808
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                   
aTCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAgc                                                                                                                                                                                                                                                                                   >  1:490876/1‑289 (MQ=255)
 tCGACCACCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGGACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGAACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAgc                                                                                                                                                                                                                                                                                   <  2:258236/288‑1 (MQ=255)
    aCCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGAGCCGGCGCAc                                                                                                                                                                                                                                                                               >  1:129978/1‑289 (MQ=255)
         cAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTaaa                                                                                                                                                                                                                                                                          >  1:579722/1‑289 (MQ=255)
                  gcgATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGcc                                                                                                                                                                                                                                                                 <  2:490876/289‑1 (MQ=255)
                    gATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCg                                                                                                                                                                                                                                                               <  1:83218/289‑1 (MQ=255)
                           aTCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCtt                                                                                                                                                                                                                                                          >  1:302549/1‑287 (MQ=255)
                           aTCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCtt                                                                                                                                                                                                                                                          <  2:302549/287‑1 (MQ=255)
                             acaCGCTCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCCCGCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTgag                                                                                                                                                                                                                                                      <  2:129978/287‑1 (MQ=255)
                                   tCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGtt                                                                                                                                                                                                                                                <  2:141628/289‑1 (MQ=255)
                                                   cagcaCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTa                                                                                                                                                                                                                                               <  2:490143/274‑1 (MQ=255)
                                                   cagcaCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTa                                                                                                                                                                                                                                               >  1:490143/1‑274 (MQ=255)
                                                           aaGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTc                                                                                                                                                                                                                        >  1:248338/1‑289 (MQ=255)
                                                                  aaTATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGGTCGGCCCTTGATTAGCGTCATCAAGCTGAAACAGCACTCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCaact                                                                                                                                                                                                                 >  2:84837/1‑287 (MQ=255)
                                                                            cccATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGCATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATg                                                                                                                                                                                                        >  1:57842/1‑288 (MQ=255)
                                                                                           aTCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGct                                                                                                                                                                                                                                                             >  1:99649/1‑220 (MQ=255)
                                                                                           aTCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGct                                                                                                                                                                                                                                                             <  2:99649/220‑1 (MQ=255)
                                                                                           aTCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         >  2:174065/1‑288 (MQ=255)
                                                                                            tCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGAACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCa                                                                                                                                                                                       <  2:371533/289‑1 (MQ=255)
                                                                                              cccAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAGGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGct                                                                                                                                                                                                                                                             <  1:527444/217‑1 (MQ=255)
                                                                                              cccAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAGGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGct                                                                                                                                                                                                                                                             >  2:527444/1‑217 (MQ=255)
                                                                                              cccAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCg                                                                                                                                                                                                                                                               >  2:383838/1‑215 (MQ=255)
                                                                                              cccAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCg                                                                                                                                                                                                                                                               <  1:383838/215‑1 (MQ=255)
                                                                                              cccAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCAc                                                                                                                                                                                      >  1:575353/1‑288 (MQ=255)
                                                                                                cAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATg                                                                                                                                                                                                        >  2:343771/1‑268 (MQ=255)
                                                                                                cAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATg                                                                                                                                                                                                        <  1:343771/268‑1 (MQ=255)
                                                                                                cAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAAACGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCActc                                                                                                                                                                                    >  2:476083/1‑288 (MQ=255)
                                                                                                    gcACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTg                                                                                                                                                                                  <  1:433697/286‑1 (MQ=255)
                                                                                                      aCGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAcc                                                                                                                                                                             <  2:436086/289‑1 (MQ=255)
                                                                                                      aCGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAAACGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAcc                                                                                                                                                                             <  1:476083/289‑1 (MQ=255)
                                                                                                      aCGTACATTGAGCGAAGAATGGAAGTTGCGTGGGCCATTTCGAACGACAAAACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGAACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAcc                                                                                                                                                                             <  2:248338/289‑1 (MQ=255)
                                                                                                                    cagcaTGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         <  2:341686/263‑1 (MQ=255)
                                                                                                                    cagcaTGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         >  1:341686/1‑263 (MQ=255)
                                                                                                                        aTGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAAcc                                                                                                                                                              <  1:278426/286‑1 (MQ=255)
                                                                                                                                 ggCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCACACAAt                                                                                                                                                  >  1:224295/1‑289 (MQ=255)
                                                                                                                                           tctcGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCa                                                                                                                                        >  1:153746/1‑289 (MQ=255)
                                                                                                                                                    acaccaCCGACTAAGCGATCCTGGGGCAGCTATGCAAGGTTTGTCGCTGACAGTGCATCTGCGACCAGAATGAAACGTTGCGGCCAGGCGTTCGGCCCTTGAGTAGCGTCATCAAGATGAAACAGCAATCGCTTACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCggg                                                                                                                                  <  2:57842/286‑1 (MQ=255)
                                                                                                                                                     caacaCCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCGACTCTGCCACCGAGCCGTTATCAACCTCGGCAAACCATTAGCGACGCAGCCGGgttg                                                                                                                               >  1:243238/1‑286 (MQ=255)
                                                                                                                                                           cGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc                                                                                                                                                                         >  2:316357/1‑240 (MQ=255)
                                                                                                                                                           cGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc                                                                                                                                                                         <  1:316357/240‑1 (MQ=255)
                                                                                                                                                                     aTCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         <  2:192963/214‑1 (MQ=255)
                                                                                                                                                                     aTCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         >  1:192963/1‑214 (MQ=255)
                                                                                                                                                                                gCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         <  2:469623/203‑1 (MQ=255)
                                                                                                                                                                                gCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                                                                                                                                                                                         >  1:469623/1‑203 (MQ=255)
                                                                                                                                                                                gCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGGGTCCGAAAGCAGGCGCGAGTAACGATAGaa                                                                                                   >  2:178574/1‑289 (MQ=255)
                                                                                                                                                                                gCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTAGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAg                                                                                                                                                                                                                                                                                    <  1:452096/112‑1 (MQ=255)
                                                                                                                                                                                gCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTAGATTAGCGTCATCAAGATGAAACAGAAATCGCTGACCCAg                                                                                                                                                                                                                                                                                    >  2:452096/1‑112 (MQ=255)
                                                                                                                                                                                       aaGAGTTGTCGCTGACATTTAATCTGCCACCAGAATGAAACGGTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAgc                                                                                             <  2:579722/288‑1 (MQ=255)
                                                                                                                                                                                        aGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAgccg                                                                                           <  2:153746/289‑1 (MQ=255)
                                                                                                                                                                                                                        aaTGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGCAAAGATCGAAAATACCCGCCGTTTCTTTTTGTTCACCGGCGGTAAAGCGt                                                           >  2:584294/1‑289 (MQ=255)
                                                                                                                                                                                                                                  tCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGcc                                                                                                                                                                                                                                                                 <  1:55164/81‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGcc                                                                                                                                                                                                                                                                 >  2:55164/1‑81 (MQ=255)
                                                                                                                                                                                                                                  tCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTGGTGCACCGGCGGCAAAGCGTTTGCTGTAgc                                                 >  2:340478/1‑289 (MQ=255)
                                                                                                                                                                                                                                                 ggCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCaa                                                                >  2:296348/1‑259 (MQ=255)
                                                                                                                                                                                                                                                 ggCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCaa                                                                <  1:296348/259‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttGATTATCGTCATCAAGATGACACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCCCTTAGCGCGGCAAACTGTTCGGCAACTTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTGGGCACACAAGTCGCGACGCAGCCGGGTCTCCGAAAGCCGGTTCGAGTAACGATCGAAAATAGCCGCCGTTTCTTTTTGTGCCCCGCCGGCAAACCTTTTGCTGTCGCGGCTACACTTGTTTCACGct                             >  2:358748/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               gTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTcccg                     >  1:247573/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                  atcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAgcgc                  >  1:31754/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                       gATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGCTTTCCGcttctctcctcgacccgg             >  2:17304/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            aaCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCCGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTc       >  2:365180/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 cAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACcg                                                                       <  1:416868/220‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACcg                                                                       >  2:416868/1‑220 (MQ=255)
                                                                                                                                                                                                                                                                                  aaTCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTcg  <  2:31754/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                   
ATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCG  >  NC_000913/2967247‑2967808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: