Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,589 T→G 6.8% G204G (GGT→GGG)  waaC → ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,5890TG 6.8% 163.0 / 2.1 68G204G (GGT→GGG) waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base G (5/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAC  >  NC_000913/3796364‑3796807
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
tCGTAAGCATCATATTGCAAAACCGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTAGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACAATGGCCGGAAGAACCACGGCGGGGATTCATTGGttt                                                                                                                                                                                                                           >  1:110341/1‑228 (MQ=255)
ccGTAAGCATCATATTGCAAAACAGCAGATCGCCGTAGAACGCGCGCACGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGGCCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTa                                                                                                                                                                                                                          <  1:635910/228‑1 (MQ=255)
  gTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTACCTACAGATGCTGGCGCATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGCGAATTGATTGGTTTACt                                                                                                                                                                                                                        >  1:519626/1‑229 (MQ=255)
          catATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGAGGATAAAAACCGGCCGGAAGAACACTGGGGGCAATTAAGTGGTTTAATGGGGGAtt                                                                                                                                                                                                                >  2:682423/1‑229 (MQ=255)
          catATTGAAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGAtt                                                                                                                                                                                                                <  1:682424/229‑1 (MQ=255)
                 cAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCGGATTCCGGac                                                                                                                                                                                                          >  1:698190/1‑227 (MQ=255)
                  aaaacgtgagcaCGCCGTAAGACGTCACCGGGAAATGGTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCTTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATa                                                                                                                                                                                                        <  1:492477/221‑1 (MQ=255)
                  aaaaCAGCAGCACGCCGTAGAACGCACCCGGGAACTGGTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAAt                                                                                                                                                                                                         <  1:247929/228‑1 (MQ=255)
                   aaaCAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATa                                                                                                                                                                                                        <  1:357090/228‑1 (MQ=255)
                         cagcaCGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATGATAAACACTCGCCGGACTAACACTTGGGGGAATTGCTTGGGTTACTAGCTGATTCAGTAATaagaaata                                                                                                                                                                                                 >  2:246727/1‑222 (MQ=255)
                              cGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTATCTGGCCTATTCAGGAAAAC‑GGATTAAACtt                                                                                                                                                                                           >  1:325915/1‑229 (MQ=255)
                                ccGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTc                                                                                                                                                                                          >  1:99273/1‑228 (MQ=255)
                                   tAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGTGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACGGGCTGATTCAGGAGTAC‑GGATTAAACTTCCGt                                                                                                                                                                                       >  1:159435/1‑228 (MQ=255)
                                          ccccccggAACTGGTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGgcgcg                                                                                                                                                                               <  1:491263/222‑1 (MQ=255)
                                          ccccccgcaattttttGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGgcgcg                                                                                                                                                                               <  2:99273/216‑1 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGCCGACCCGTGATGATAAACACTGGCCGGAAGAACACAGGGCAGAATGTATTGGTTTACTGGCTGATTCACGAACAC‑GGAATAAACATCCGTGGGgcgcgc                                                                                                                                                                              >  2:743503/1‑229 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCAGGAAGAACACTGGCGGGATTTGTTTGGGTTACTGGATGATTCAGGAAAAC‑GGTTTAAACATACGTGGGgcgcgc                                                                                                                                                                              >  2:313762/1‑229 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTAATGATAAACACTGGCCGGAAGAAAACTGTCGAGAATTGATTGGTTTACTGGCTGATTTAGGAAAAC‑GGTATAAAATTCCTGggggggggc                                                                                                                                                                              >  2:491263/1‑224 (MQ=255)
                                            cccggcaacttgtttccAAAAGTTTTGGCTATACCAAACCGCAAAACCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGcgcc                                                                                                                                                                             <  1:726817/214‑1 (MQ=255)
                                             ccGGGATCTTTTTTCCAAAAGTTTGGGTTTTAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGcgccg                                                                                                                                                                            <  2:325915/229‑1 (MQ=255)
                                             ccGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGTCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGGTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTTCGTGGGGCGcgccc                                                                                                                                                                            >  1:755660/1‑228 (MQ=255)
                                                         tGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGGGAATTGATTGGTTTACTTGCTGATTCAGGAATAC‑CGCTTAAACTTCCGTGGGGCGCGCCGCATGAGgaagaa                                                                                                                                                                >  1:691446/1‑229 (MQ=255)
                                                                           taGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGCAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCCGGAATAA‑GGATTAAACTTCCGGTGGGCGCGCCGCATGAGGAGGAACGGGGGAACCGACTgggg                                                                                                                                              >  1:703085/1‑227 (MQ=255)
                                                                                 accgcaaaccCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCCGGGGGCGCGCCGCATGAGGAAGAACGGGCGCAACGACTGGCGGAAggc                                                                                                                                        >  1:478214/1‑228 (MQ=255)
                                                                                 accgcaaaccCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGGGATGCTAAACACAGGCCGGAAGAACACTGGCGAACATTGATTGGGTTCCTGGCTGATTCAGGAAAAC‑GGGTTAAAACTCCTTGGGGCGCGCCGAATGAGGAAGAAGGCGCGAAACGAACGGCGAAAGGa                                                                                                                                        >  2:188885/1‑229 (MQ=255)
                                                                                      aaaccCAGGGCGATTATGCATTTGCACAGCATTGTATGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGc                                                                                                                                   <  1:188885/229‑1 (MQ=255)
                                                                                            cgggCGATTTTGTTATCTCCCAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACCCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACCGGCGAAACGCCTGGCGGAAGGATTTGCTTATGt                                                                                                                             <  2:358600/228‑1 (MQ=255)
                                                                                             gggCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGGTTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGAGCCCCCTGAGGAAGAACGGGCGAAACGACAGGCGGCAGGGTTTGCTTATGtt                                                                                                                            >  2:579489/1‑229 (MQ=255)
                                                                                             gggCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTTAGGGATAA‑GAGTTAAACTTCCGTGGGGCGCGCCGCCTGCGGGAGAAAGGGCGAAAAGAACGGCGGGAgggattgttttggtt                                                                                                                            >  2:232103/1‑216 (MQ=255)
                                                                                                  aTTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGCTCGATTCAGGAAAAC‑GGATTAAACTTCCGTGGGGGCCGCCGCATGAGGAAGAACGGGGGAAACGAATGGGGGAAGGATTTGCTTATGTTGAGGt                                                                                                                       >  1:42767/1‑229 (MQ=255)
                                                                                                    tATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGCAACGACTGGCGGAAGGATTTGCTTATGTTTAAggaa                                                                                                                     >  1:291195/1‑226 (MQ=255)
                                                                                                           tCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGc                                                                                                                                                                                                                     >  2:361805/1‑127 (MQ=255)
                                                                                                           tCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGc                                                                                                                                                                                                                     <  1:361805/127‑1 (MQ=255)
                                                                                                            cGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGCTTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTCGGGCGGGCCGCATGAGGAAGAACGGGGGAAACGACTGGCGGACGGAGTTGGTTTGGTTGGAGAAATGCCGaaa                                                                                                             >  2:255053/1‑228 (MQ=255)
                                                                                                             ccccacttttttCTG‑CGAACCTGCCTAGAAATGCGCGGGAATATGCCGTATTTCTTCATGCGACGGCCCGTGATGATAAACCCTGGCCGGAAGAACCCTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAt                                                                                                           <  2:42767/221‑1 (MQ=255)
                                                                                                              tccaGCATTTTCGGAAGAACCTGCCTACAAGTTCGCGCGAAGATGCCGTATTTCTTCATGCGACGACCCGTGATGTTAAACCCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAt                                                                                                           <  2:495953/227‑1 (MQ=255)
                                                                                                                caGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACCGGCGAGAATTGATTGGCTTACTGGCTGATTCAGGAATAC‑GGGTTAAACTTCCGTGGGGCGCGCCGCATTAGGAATAACGGGCGAAAAGACGTGCGGAAGGGGTGGGTTGTGTTGAAgaaatggcggaaaggg                                                                                                         >  2:202857/1‑214 (MQ=255)
                                                                                                                              gacgcgtccTAAAGGTGTTGGCGAATAGGCCGTTTTTCGTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGt                                                                                            <  1:514093/221‑1 (MQ=255)
                                                                                                                                  cggCCTAAAGATGCTGGCGAATATGCCGTATGTCTTCATGCGACGACCCGTGAGGATAAACGCTGGCAGGAAGAACACTGGCGAGAATTGATTGGGTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGccc                                                                                       <  1:441865/227‑1 (MQ=255)
                                                                                                                                                 ggCGAATATGCCGTATTTCTTCATGCGCCGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGCGTTGATTGGTTTACTGGCTGATACAGGAATAC‑GGATTAAACTTACTACGGGCGCGCCGCGTGAGGAAGAACGGCCGAAACGACTGTCGGAATGCATTGATTATGTTGGAGGATGGGCGGAGATGATAGTGGAAGGCGtgtcacccggacggggcgggg                                                                        >  2:664230/1‑209 (MQ=255)
                                                                                                                                                  gCGGATTTGCCGTATTTCTTCATGCGACGGCCAGGGATGATAAACAATGGCCGGAATAACACTGGCGAGAATGGATTGGTTTACTGGCTGATTCAGGAACAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGc                                                                       <  1:750148/229‑1 (MQ=255)
                                                                                                                                                             ggaTTGCTTCATGCGACGACCCGTGATGATAAACGCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAGCTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTa                                                             <  1:202857/226‑1 (MQ=255)
                                                                                                                                                                tttCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGCCTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGACGATGAGGCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTCAATTTTTAgggt                                                         >  1:35896/1‑229 (MQ=255)
                                                                                                                                                                 ttcttcATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTc                                                        <  1:46914/229‑1 (MQ=255)
                                                                                                                                                                    ttccctgcacggcCCGTGGTGATAAACACTGGCCGTAAGAACACCGGCGAGAATTGATTGGTTTACTGGCTGGTTCAGGAATCG‑TGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGt                                                     <  2:110341/221‑1 (MQ=255)
                                                                                                                                                                     tcAGGCGACGACCCGTGATTATAAACCCTGGCCGGAAGAACACTGGCGAGGATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCgg                                                      <  1:642938/227‑1 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGGGTTGCCCTCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGAt                                                 >  1:585650/1‑229 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAgtg                                                          <  2:205964/220‑1 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAgtg                                                          >  1:205964/1‑220 (MQ=255)
                                                                                                                                                                          cgactaCCCGTGATGATAAACACTGGCCGGAAGAACGCTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGGATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATAc                                               <  1:653605/229‑1 (MQ=255)
                                                                                                                                                                                   tgtatgatAATCACTGGCCGGAAGACCACTGGCGAGAATTGATTGGTTTGCTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAg                                      <  1:316090/226‑1 (MQ=255)
                                                                                                                                                                                           aaaCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTTAAGGATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCCGGTGCTGGCCGGGGCGAAATTTGTATGGTAGGGGGATAAGGGGGTagaacaattaa                              >  2:117985/1‑219 (MQ=255)
                                                                                                                                                                                           aaaCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAAGGTCGGGGGGTACGGGGTTAAGCCATTTaa                              >  1:250686/1‑229 (MQ=255)
                                                                                                                                                                                              actgtgCCGGAAAAACACTGGCGAGAAGGGAGTGGTTTACTGGCTGCTGCAGG‑ATACGGGATTAAACTTCCGGGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAATTATTGGCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTTTAGAGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  2:755659/224‑1 (MQ=255)
                                                                                                                                                                                               ccTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:626704/227‑1 (MQ=255)
                                                                                                                                                                                               ccTGGCCGGAAGAAAACGGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTGCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:60213/227‑1 (MQ=255)
                                                                                                                                                                                                 tGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGGTTTGCTTATGTTGAAGAATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGAGTGCTGGCCGGGGGTAAATTTGTAGGGTCGGTGGATAAGGGGTTTAGCCATTTAAaggcgg                        >  1:235583/1‑229 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACCCTGGCTAGGATTGGTTGGTTTACTGGCTGATTCAGGAATCC‑GGATTAAACTTCCGTGGGGCGAGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCa                      <  1:313762/229‑1 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGGGTTGTCCGCTGGCCGGCCGGGGCTAAATTTGTAGTGTCGGTGGTGAAGGGGGTAAGGCATTTAACgggggc                       >  1:744795/1‑228 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCTCGTGCGGGCCGGGGCTAACTTTGTAGTGTCGTGGGAGTCGGGGGTAAGCCCTTTAACggcggcg                      >  1:345438/1‑228 (MQ=255)
                                                                                                                                                                                                        aagaaCATTGGGGAGAATTGATTGGTTTACGGGCTGATTCAGGAAGAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGa                 <  2:291195/229‑1 (MQ=255)
                                                                                                                                                                                                         agaaCCCTGGGGAGAATTGATTGTTTTACTGGCTGATTCAGGAATCC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGAt                <  2:703084/229‑1 (MQ=255)
                                                                                                                                                                                                                  gCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           >  1:690520/1‑209 (MQ=255)
                                                                                                                                                                                                                  gCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  2:690519/209‑1 (MQ=255)
                                                                                                                                                                                                                     agaATTGATTGGGTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAGGGGTCTGGAAGGCGTTGCCCCCGGGTTGGCCGGGGCCAAATTTGTGAGGTAGTGGGATTCGGGGGTACGCCCTTTATCGGCGGGACTgagaagagacaaaaaa    >  2:313537/1‑214 (MQ=255)
                                                                                                                                                                                                                      gtattgattgTGTTGCTGGGGGTTTCAGGAATAC‑GGAATAACCTTTCGTGGGGCGCGCCGAATGATTAAGAACGGGCGAAACGTATGGCGGAAGGATTTGCATATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCa   <  2:29139/227‑1 (MQ=255)
                                                                                                                                                                                                                       gattgattgGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCACATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc  <  1:172329/228‑1 (MQ=255)
                                                                                                                                                                                                                       atttgtttgGTTTACTGGCGAATTCAGGAATAC‑GGATTAAACTTCCGGGGGGCGCGCCGAATGAGGAAGAACGGGAGAAACGGCTGGCGGAAGGATTTGCTTATGTTGAAGTATGGCCGAAGATGAGGCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc  <  2:519626/227‑1 (MQ=255)
                                                                                                                                                                                                                            attgGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:575908/199‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
TCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATAC‑GGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAC  >  NC_000913/3796364‑3796807

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