Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,055,078 C→A 8.1% K67N (AAG→AAT glnG ← fused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,055,0780CA 8.1% 138.6 / 2.7 46K67N (AAG→AATglnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/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 = 7.41e-01

TCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAA  >  NC_000913/4054783‑4055377
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
tCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGaa                                                                                                                                                                                                                                                                                                        >  2:4935/1‑300 (MQ=255)
     aCGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAgcg                                                                                                                                                                                                                                                                                                   <  2:71409/301‑1 (MQ=255)
        cTAATAGATGTACGCGAAAGCCTACCGATAATACGGATCACGTCCTGCATGGCTGGCGCTTCGCCGATGATACCGTTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGCTCAAGCAGATTAAACAGCGCCa                                                                                                                                                                                                                                                                                                <  2:625503/301‑1 (MQ=255)
         tcttccgatctcgcgAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGCTCAAGCAGATTAAACAGCGCCat                                                                                                                                                                                                                                                                                               <  2:237372/290‑2 (MQ=255)
          aaTAGAAGAACGCGAAAGCCGACCGATAATACGGAACAGTGCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGc                                                                                                                                                                                                                                                                                              <  1:128532/301‑1 (MQ=255)
              aaagaaCGCGAAAGCTGACGGTTAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGt                                                                                                                                                                                                                                                                                          <  1:4935/300‑1 (MQ=255)
               tcagacgtgtgctcttcctatttataCGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTc                                                                                                                                                                                                                                                                                         <  2:584561/278‑1 (MQ=255)
                tgaaCTCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGTCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTcc                                                                                                                                                                                                                                                                                        <  2:278847/300‑1 (MQ=255)
                 gaaCGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGAATTGGATGACGCTGTTTAATATACTTGAGCAGCGCCAGCCCGTccc                                                                                                                                                                                                                                                                                       >  1:278544/1‑300 (MQ=255)
                       gAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTTCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCg                                                                                                                                                                                                                                                                                 <  2:361156/301‑1 (MQ=255)
                            ccgaccgaTAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGAGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATACAACGCCCCTTGTTGATAGGCGCTGACTGACGCATCCAGATCGGAATGTACGGTCATAATGAATACCGGAAGCAATGGAAGGCGCCGTATAATCTGATTGaaccacacaccccagtccattaccggccat                                                                                                                                                                                                                                                                            >  1:374848/1‑272 (MQ=255)
                                      ataCGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCCCTGACGGCAGCATCCAGCTCGGCATGTGCGGTCATACTGATGACCGGACGCATTAGATGGCGCTGATTAATCTGCTTGAGCCGCGCCAGCCCGTCAACTccgaaaagagccgtaccag                                                                                                                                                                                                                                                                  >  1:583837/1‑284 (MQ=255)
                                          ggAACACTTCCTTCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAag                                                                                                                                                                                                                                                              <  2:626645/301‑1 (MQ=255)
                                              cacGTCCTGCTTGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAagcagc                                                                                                                                                                                                                                                          <  2:333025/301‑1 (MQ=255)
                                                    cTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGCAGCATTCGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCACGTCATACGGATAGCTAAAAGACGCCCATcc                                                                                                                                                                                                                                                    >  2:359986/1‑301 (MQ=255)
                                                     tGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAATCATTGGATTGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATACg                                                                                                                                                                                                                                                   >  2:306580/1‑301 (MQ=255)
                                                         gggctggcGGTTCGCCGTTGTTATCGGGCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCGAATCCAGCGCCACTGCTTCGTCGATATCAAACGGTGTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGt                                                                                                                                                                                                                                               <  1:640255/300‑1 (MQ=255)
                                                             ggtgggttggtcggtgtgttCGGTGTTTGGGGCTTTAAGCTGACAATTATGCGGCTGCTGCGTTTATTGGTAATGACTTAGAGCTCGATCAACCAGCGCCACTGCTTCGTCGATATCAAACGTTTTGGGCAGATAATCAATCGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAAGGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCTGCGTTTTg                                                                                                                                                                                                                                           <  1:375315/282‑1 (MQ=255)
                                                                              tatCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAg                                                                                                                                                                                                                          <  2:89193/301‑1 (MQ=255)
                                                                                                                                      tGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGATCTCAAACGTCGTACAGGTTAAAACTGACCCagacagcgaaa                                                                                                                                                                  >  2:367106/1‑298 (MQ=255)
                                                                                                                                           gATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTAGAACAGGTTAAACATGACCCAGCGAACGAACGTTCa                                                                                                                                                             >  2:343771/1‑301 (MQ=255)
                                                                                                                                             aggcgcggTCAACCGGCTCAACTGCTTCGTCGTTATCCAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGTGAGGCGCTGACTGCAGCAGCCAGATCGGAATGTGCGGTCATAAGGATGACCGGAAGCAGTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAg                                                                                                                                                           <  1:457547/299‑1 (MQ=255)
                                                                                                                                                       aaCCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGCAAGCAGCCCATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTACCCCAGCGACCGCCCGTACCACCACCCAACgg                                                                                                                                                 >  1:524533/1‑301 (MQ=255)
                                                                                                                                                         ccAGTGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCGCGTTCAAGCACCCAACGGAt                                                                                                                                               <  2:71673/301‑1 (MQ=255)
                                                                                                                                                              gcCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGACAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGAtaaaaa                                                                                                                                          >  2:157190/1‑296 (MQ=255)
                                                                                                                                                                   gggttggtgggtatCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCg                                                                                                                                     <  1:306176/290‑1 (MQ=255)
                                                                                                                                                                         gtcgATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCCCAAACGTCGTACAGGTTAAACCTGCCCAAGCGAGCGCCCCTTCAAGCACCCAACGGATGAAACTATCGTCaacc                                                                                                                               >  1:548384/1‑298 (MQ=255)
                                                                                                                                                                                            ggCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCACTc                                                                                                            >  2:566554/1‑301 (MQ=255)
                                                                                                                                                                                                aaTAATCAATCGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTg                                                                                                        <  1:565840/300‑1 (MQ=255)
                                                                                                                                                                                                   aaTCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCAAGACTATCACTCGTTTCaa                                                                                                     >  2:503728/1‑300 (MQ=255)
                                                                                                                                                                                                                 ggttgAGAGGCGCTGACGGAAGCATCCAGTTCGGAATGTGCGGTCATACTGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCGTGAGCAGCGCCAGCCCGTCCGTTCCCGGCATACGGATATCTGAAATCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTa                                                                                       <  1:240214/300‑1 (MQ=255)
                                                                                                                                                                                                                   ttgATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCCAACGTCGTCCAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACACAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCACGTTGCATAAACGTCACCTTTAtt                                                                                     >  1:488253/1‑301 (MQ=255)
                                                                                                                                                                                                                                 cGGCAGCTTCTGGGTCGTAATGTGCGGTCTGATTGATGTCCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGGGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAgg                                                                       <  1:343308/301‑1 (MQ=255)
                                                                                                                                                                                                                                          ccagtgcGGAATGTGTGGTTAGAATGATGACCGGAAGCATTGGATGGCGCTGTTTAAGCTGCTGTGGCAGCGACAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgag                                                              <  1:157039/295‑1 (MQ=255)
                                                                                                                                                                                                                                             acgacgctcttccgatcttAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAc                                                           <  1:36079/283‑1 (MQ=255)
                                                                                                                                                                                                                                                    gcttccGGTCATCATTATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAAAACTCGGTaa                                                    >  2:610823/6‑300 (MQ=255)
                                                                                                                                                                                                                                                    gcttccGGTCATCATTATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGATAAAACGAGAACTCGGTaa                                                    >  2:174605/6‑300 (MQ=255)
                                                                                                                                                                                                                                                      ttccGGTCATCATTATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGAGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGAGAAACCTGCACCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATAGTCAGCGACTACCCAGACTAACACTCGTGGAATAAAAGTCACATTTAATGACTGATAGGAAGGAAAACCGAACACTCGGTAAGc                                                  >  2:542586/4‑301 (MQ=255)
                                                                                                                                                                                                                                                       tccGGTCATCATTATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGACCCCTTTATTTCCTGATAGGCAGGTAAACCGAGAACCCGGCATGcc                                                 >  1:362264/3‑301 (MQ=255)
                                                                                                                                                                                                                                                               tAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAcagaacagaagagcaaaca                                         >  2:36130/1‑283 (MQ=255)
                                                                                                                                                                                                                                                               tAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACACAGACTATCCCTCGTTGCCTCAAAGTCACCTTTATTTCCTGAAAGGCAGGTAAACAGACAACACGGTATGCACTGACCaa                                         >  1:492365/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             aGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAt                           <  2:525205/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ctgggggtGGCGCTGTTGAATCTGCTTGTGCCGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAATCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAttt                         <  1:8149/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAtta                         >  2:8149/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                     ctGGCGCTGCTCAAGCAGATTAAACAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATACACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCAGGCCAACTGGAAAATTCAAttttcactt                   >  1:50187/2‑296 (MQ=255)
                                                                                                                                                                                                                                                                                                      gCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTAAATTTTGCCTGAATGCTGATCAATCaaa  >  2:299494/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
TCAGCACGCTAATAGAAGAACGCGAAAGCCGACCGATAATACGGAACACGTCCTGCATGGCTGGCGCTTCGCCGATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTGACGGCAGCATCCAGATCGGAATGTGCGGTCATAATGATGACCGGAAGCATTGGATGGCGCTGTTTAATCTGCTTGAGCAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAA  >  NC_000913/4054783‑4055377

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