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→GGGwaaCADP‑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

TCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAC  >  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)
                              cGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTATCTGGCCTATTCAGGAAAACGGATTAAACtt                                                                                                                                                                                           >  1:325915/1‑229 (MQ=255)
                                ccGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTc                                                                                                                                                                                          >  1:99273/1‑228 (MQ=255)
                                   tAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGTGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACGGGCTGATTCAGGAGTACGGATTAAACTTCCGt                                                                                                                                                                                       >  1:159435/1‑228 (MQ=255)
                                          ccccccggAACTGGTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGgcgcg                                                                                                                                                                               <  1:491263/222‑1 (MQ=255)
                                          ccccccgcaattttttGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGgcgcg                                                                                                                                                                               <  2:99273/216‑1 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGCCGACCCGTGATGATAAACACTGGCCGGAAGAACACAGGGCAGAATGTATTGGTTTACTGGCTGATTCACGAACACGGAATAAACATCCGTGGGgcgcgc                                                                                                                                                                              >  2:743503/1‑229 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCAGGAAGAACACTGGCGGGATTTGTTTGGGTTACTGGATGATTCAGGAAAACGGTTTAAACATACGTGGGgcgcgc                                                                                                                                                                              >  2:313762/1‑229 (MQ=255)
                                           aCCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTAATGATAAACACTGGCCGGAAGAAAACTGTCGAGAATTGATTGGTTTACTGGCTGATTTAGGAAAACGGTATAAAATTCCTGggggggggc                                                                                                                                                                              >  2:491263/1‑224 (MQ=255)
                                            cccggcaacttgtttccAAAAGTTTTGGCTATACCAAACCGCAAAACCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGcgcc                                                                                                                                                                             <  1:726817/214‑1 (MQ=255)
                                             ccGGGATCTTTTTTCCAAAAGTTTGGGTTTTAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGcgccg                                                                                                                                                                            <  2:325915/229‑1 (MQ=255)
                                             ccGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGTCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGGTTACTGGCTGATTCAGGAATACGGATTAAACTTTCGTGGGGCGcgccc                                                                                                                                                                            >  1:755660/1‑228 (MQ=255)
                                                         tGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGGGAATTGATTGGTTTACTTGCTGATTCAGGAATACCGCTTAAACTTCCGTGGGGCGCGCCGCATGAGgaagaa                                                                                                                                                                >  1:691446/1‑229 (MQ=255)
                                                                           taGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGCAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCCGGAATAAGGATTAAACTTCCGGTGGGCGCGCCGCATGAGGAGGAACGGGGGAACCGACTgggg                                                                                                                                              >  1:703085/1‑227 (MQ=255)
                                                                                 accgcaaaccCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCCGGGGGCGCGCCGCATGAGGAAGAACGGGCGCAACGACTGGCGGAAggc                                                                                                                                        >  1:478214/1‑228 (MQ=255)
                                                                                 accgcaaaccCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGGGATGCTAAACACAGGCCGGAAGAACACTGGCGAACATTGATTGGGTTCCTGGCTGATTCAGGAAAACGGGTTAAAACTCCTTGGGGCGCGCCGAATGAGGAAGAAGGCGCGAAACGAACGGCGAAAGGa                                                                                                                                        >  2:188885/1‑229 (MQ=255)
                                                                                      aaaccCAGGGCGATTATGCATTTGCACAGCATTGTATGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGc                                                                                                                                   <  1:188885/229‑1 (MQ=255)
                                                                                            cgggCGATTTTGTTATCTCCCAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACCCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACCGGCGAAACGCCTGGCGGAAGGATTTGCTTATGt                                                                                                                             <  2:358600/228‑1 (MQ=255)
                                                                                             gggCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGGTTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGAGCCCCCTGAGGAAGAACGGGCGAAACGACAGGCGGCAGGGTTTGCTTATGtt                                                                                                                            >  2:579489/1‑229 (MQ=255)
                                                                                             gggCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTTAGGGATAAGAGTTAAACTTCCGTGGGGCGCGCCGCCTGCGGGAGAAAGGGCGAAAAGAACGGCGGGAgggattgttttggtt                                                                                                                            >  2:232103/1‑216 (MQ=255)
                                                                                                  aTTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGCTCGATTCAGGAAAACGGATTAAACTTCCGTGGGGGCCGCCGCATGAGGAAGAACGGGGGAAACGAATGGGGGAAGGATTTGCTTATGTTGAGGt                                                                                                                       >  1:42767/1‑229 (MQ=255)
                                                                                                    tATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGCAACGACTGGCGGAAGGATTTGCTTATGTTTAAggaa                                                                                                                     >  1:291195/1‑226 (MQ=255)
                                                                                                           tCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGc                                                                                                                                                                                                                     >  2:361805/1‑127 (MQ=255)
                                                                                                           tCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGc                                                                                                                                                                                                                     <  1:361805/127‑1 (MQ=255)
                                                                                                            cGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGCTTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTCGGGCGGGCCGCATGAGGAAGAACGGGGGAAACGACTGGCGGACGGAGTTGGTTTGGTTGGAGAAATGCCGaaa                                                                                                             >  2:255053/1‑228 (MQ=255)
                                                                                                             ccccacttttttCTGCGAACCTGCCTAGAAATGCGCGGGAATATGCCGTATTTCTTCATGCGACGGCCCGTGATGATAAACCCTGGCCGGAAGAACCCTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAt                                                                                                           <  2:42767/221‑1 (MQ=255)
                                                                                                              tccaGCATTTTCGGAAGAACCTGCCTACAAGTTCGCGCGAAGATGCCGTATTTCTTCATGCGACGACCCGTGATGTTAAACCCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAt                                                                                                           <  2:495953/227‑1 (MQ=255)
                                                                                                                caGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACCGGCGAGAATTGATTGGCTTACTGGCTGATTCAGGAATACGGGTTAAACTTCCGTGGGGCGCGCCGCATTAGGAATAACGGGCGAAAAGACGTGCGGAAGGGGTGGGTTGTGTTGAAgaaatggcggaaaggg                                                                                                         >  2:202857/1‑214 (MQ=255)
                                                                                                                              gacgcgtccTAAAGGTGTTGGCGAATAGGCCGTTTTTCGTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGt                                                                                            <  1:514093/221‑1 (MQ=255)
                                                                                                                                  cggCCTAAAGATGCTGGCGAATATGCCGTATGTCTTCATGCGACGACCCGTGAGGATAAACGCTGGCAGGAAGAACACTGGCGAGAATTGATTGGGTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGccc                                                                                       <  1:441865/227‑1 (MQ=255)
                                                                                                                                                 ggCGAATATGCCGTATTTCTTCATGCGCCGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGCGTTGATTGGTTTACTGGCTGATACAGGAATACGGATTAAACTTACTACGGGCGCGCCGCGTGAGGAAGAACGGCCGAAACGACTGTCGGAATGCATTGATTATGTTGGAGGATGGGCGGAGATGATAGTGGAAGGCGtgtcacccggacggggcgggg                                                                        >  2:664230/1‑209 (MQ=255)
                                                                                                                                                  gCGGATTTGCCGTATTTCTTCATGCGACGGCCAGGGATGATAAACAATGGCCGGAATAACACTGGCGAGAATGGATTGGTTTACTGGCTGATTCAGGAACACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGc                                                                       <  1:750148/229‑1 (MQ=255)
                                                                                                                                                             ggaTTGCTTCATGCGACGACCCGTGATGATAAACGCTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAGCTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTa                                                             <  1:202857/226‑1 (MQ=255)
                                                                                                                                                                tttCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGCCTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGACGATGAGGCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTCAATTTTTAgggt                                                         >  1:35896/1‑229 (MQ=255)
                                                                                                                                                                 ttcttcATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTc                                                        <  1:46914/229‑1 (MQ=255)
                                                                                                                                                                    ttccctgcacggcCCGTGGTGATAAACACTGGCCGTAAGAACACCGGCGAGAATTGATTGGTTTACTGGCTGGTTCAGGAATCGTGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGt                                                     <  2:110341/221‑1 (MQ=255)
                                                                                                                                                                     tcAGGCGACGACCCGTGATTATAAACCCTGGCCGGAAGAACACTGGCGAGGATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCgg                                                      <  1:642938/227‑1 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGGGTTGCCCTCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGAt                                                 >  1:585650/1‑229 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAgtg                                                          <  2:205964/220‑1 (MQ=255)
                                                                                                                                                                        tGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAgtg                                                          >  1:205964/1‑220 (MQ=255)
                                                                                                                                                                          cgactaCCCGTGATGATAAACACTGGCCGGAAGAACGCTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGGATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATAc                                               <  1:653605/229‑1 (MQ=255)
                                                                                                                                                                                   tgtatgatAATCACTGGCCGGAAGACCACTGGCGAGAATTGATTGGTTTGCTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAg                                      <  1:316090/226‑1 (MQ=255)
                                                                                                                                                                                           aaaCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTTAAGGATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCCGGTGCTGGCCGGGGCGAAATTTGTATGGTAGGGGGATAAGGGGGTagaacaattaa                              >  2:117985/1‑219 (MQ=255)
                                                                                                                                                                                           aaaCACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAAGGTCGGGGGGTACGGGGTTAAGCCATTTaa                              >  1:250686/1‑229 (MQ=255)
                                                                                                                                                                                              actgtgCCGGAAAAACACTGGCGAGAAGGGAGTGGTTTACTGGCTGCTGCAGGATACGGGATTAAACTTCCGGGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAATTATTGGCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTTTAGAGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  2:755659/224‑1 (MQ=255)
                                                                                                                                                                                               ccTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:626704/227‑1 (MQ=255)
                                                                                                                                                                                               ccTGGCCGGAAGAAAACGGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTGCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:60213/227‑1 (MQ=255)
                                                                                                                                                                                                 tGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGGTTTGCTTATGTTGAAGAATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGAGTGCTGGCCGGGGGTAAATTTGTAGGGTCGGTGGATAAGGGGTTTAGCCATTTAAaggcgg                        >  1:235583/1‑229 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACCCTGGCTAGGATTGGTTGGTTTACTGGCTGATTCAGGAATCCGGATTAAACTTCCGTGGGGCGAGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCa                      <  1:313762/229‑1 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGGGTTGTCCGCTGGCCGGCCGGGGCTAAATTTGTAGTGTCGGTGGTGAAGGGGGTAAGGCATTTAACgggggc                       >  1:744795/1‑228 (MQ=255)
                                                                                                                                                                                                   gCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCTCGTGCGGGCCGGGGCTAACTTTGTAGTGTCGTGGGAGTCGGGGGTAAGCCCTTTAACggcggcg                      >  1:345438/1‑228 (MQ=255)
                                                                                                                                                                                                        aagaaCATTGGGGAGAATTGATTGGTTTACGGGCTGATTCAGGAAGACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGa                 <  2:291195/229‑1 (MQ=255)
                                                                                                                                                                                                         agaaCCCTGGGGAGAATTGATTGTTTTACTGGCTGATTCAGGAATCCGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGAt                <  2:703084/229‑1 (MQ=255)
                                                                                                                                                                                                                  gCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           >  1:690520/1‑209 (MQ=255)
                                                                                                                                                                                                                  gCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  2:690519/209‑1 (MQ=255)
                                                                                                                                                                                                                     agaATTGATTGGGTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGAGGGGTCTGGAAGGCGTTGCCCCCGGGTTGGCCGGGGCCAAATTTGTGAGGTAGTGGGATTCGGGGGTACGCCCTTTATCGGCGGGACTgagaagagacaaaaaa    >  2:313537/1‑214 (MQ=255)
                                                                                                                                                                                                                      gtattgattgTGTTGCTGGGGGTTTCAGGAATACGGAATAACCTTTCGTGGGGCGCGCCGAATGATTAAGAACGGGCGAAACGTATGGCGGAAGGATTTGCATATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCa   <  2:29139/227‑1 (MQ=255)
                                                                                                                                                                                                                       gattgattgGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCACATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc  <  1:172329/228‑1 (MQ=255)
                                                                                                                                                                                                                       atttgtttgGTTTACTGGCGAATTCAGGAATACGGATTAAACTTCCGGGGGGCGCGCCGAATGAGGAAGAACGGGAGAAACGGCTGGCGGAAGGATTTGCTTATGTTGAAGTATGGCCGAAGATGAGGCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc  <  2:519626/227‑1 (MQ=255)
                                                                                                                                                                                                                            attgGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAAcgg                           <  1:575908/199‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
TCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAC  >  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: