Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,073,967 T→G 6.5% E15D (GAA→GAC rutA ← pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,073,9670TG 6.5% 210.8 / 7.5 62E15D (GAA→GACrutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base G (2/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

TGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTGTTGCCAATCTACAAGAGGGGAGAGCGCATGACGCAAGGCGC  >  NC_000913/1073699‑1074255
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                
tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAg                                                                                                                                                                                                                                                                               >  1:465357/1‑288 (MQ=255)
tGGCGGAAATCGGAAGGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGGGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAg                                                                                                                                                                                                                                                                               <  2:394013/288‑1 (MQ=255)
      aaaTCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACg                                                                                                                                                                                                                                                                        <  2:48898/289‑1 (MQ=255)
      aaaTCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACg                                                                                                                                                                                                                                                                        <  2:465357/289‑1 (MQ=255)
      aaaTCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACg                                                                                                                                                                                                                                                                        <  2:188333/289‑1 (MQ=255)
        aTCTGAATGCGCGAGGTCACGCCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACggg                                                                                                                                                                                                                                                                      <  1:444372/289‑1 (MQ=255)
             aaTGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGACTACGCCAATCTTCATCAACAAACTCTCTTCGTCTCGTATAGGGAACGTCAGACgggcgcgt                                                                                                                                                                                                                                                                 >  2:444372/1‑284 (MQ=255)
                  gcgGGGTCACGGCCCCCAGCCCCGCCATGAAGGTGAACGACTCAAGGTTGTGATCCCAGAACGCAGGTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             <  2:348044/288‑1 (MQ=255)
                  gcgAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             <  2:8099/288‑1 (MQ=255)
                  gcgAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             <  1:673334/288‑1 (MQ=255)
                  gcgAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACGCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             <  2:680484/288‑1 (MQ=255)
                   cgAGGTACCGGCCGCCAGCCCCGCCATCAAGGTGAACGACCAAAGGTTGTGAGCCCGGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCGTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGTATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATc                                                                                                                                                                                                                                                            <  2:38495/288‑1 (MQ=255)
                         cACGTCCGCCAGCCCCGCCAACAGTGGGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTGGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTGGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCata                                                                                                                                                                                                                                                     <  2:67504/289‑1 (MQ=255)
                                 ccAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGc                                                                                                                                                                                                                                                        >  2:460192/1‑278 (MQ=255)
                                 ccAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGc                                                                                                                                                                                                                                                        <  1:460192/278‑1 (MQ=255)
                                        cGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGcc                                                                                                                                                                                                                                         >  1:147916/1‑286 (MQ=255)
                                        cGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGcc                                                                                                                                                                                                                                         <  2:147916/286‑1 (MQ=255)
                                                     aaCGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGttt                                                                                                                                                                                                                         >  1:329883/1‑289 (MQ=255)
                                                        gACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCa                                                                                                                                                                                                                      <  1:130022/289‑1 (MQ=255)
                                                                    tgATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGgcg                                                                                                                                                                                                                                                                   <  2:296168/232‑1 (MQ=255)
                                                                    tgATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGgcg                                                                                                                                                                                                                                                                   >  1:296168/1‑232 (MQ=255)
                                                                    tgATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAACCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTAATCATCAACCTCACTTCGTCCCGTAAAGTGAAACTCAGACGGGGCGCTGCATACTGCATATACTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGcaa                                                                                                                                                                                                          >  2:616039/1‑287 (MQ=255)
                                                                         ccAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGtttt                                                                                                                                                                                                     >  1:20648/1‑289 (MQ=255)
                                                                                  aGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             >  2:458978/1‑224 (MQ=255)
                                                                                  aGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAt                                                                                                                                                                                                                                                             <  1:458978/224‑1 (MQ=255)
                                                                                              aaaGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAAt                                                                                                                                                                                <  2:688208/289‑1 (MQ=255)
                                                                                              aaaGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAAt                                                                                                                                                                                <  2:20648/289‑1 (MQ=255)
                                                                                                  ccACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGt                                                                                                                                                                            <  2:329883/289‑1 (MQ=255)
                                                                                                                    gACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGtt                                                                                                                                                                <  1:343382/283‑1 (MQ=255)
                                                                                                                    gACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGtt                                                                                                                                                                >  2:343382/1‑283 (MQ=255)
                                                                                                                         ggCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAACAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAatat                                                                                                                                                     >  2:478046/1‑289 (MQ=255)
                                                                                                                         ggCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAata                                                                                                                                                      >  1:508287/1‑288 (MQ=255)
                                                                                                                          gCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAAGTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAatat                                                                                                                                                     <  2:705243/288‑1 (MQ=255)
                                                                                                                               atcgacatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCAAATAGGTACGAATACGTCAATTTTTATCATCAACCTCTCTTCGTCTCCTAAAGTGAAAGTCAGATGGCGCGCTGCATCCTGCATATCCTTTTCAGCCGCGGATTGGCTTGTTTGCACAGCGGATGCACATTTTTAAAAAGTTACTGTTGTTAATCTGTTAAAATTACTTTATCACAAAGATCTGTTa                                                                                                                                               >  2:556095/1‑289 (MQ=255)
                                                                                                                                         tAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGa                                                                                                                                      <  1:120033/288‑1 (MQ=255)
                                                                                                                                            tGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTaa                                                                                                                                  >  2:102039/1‑289 (MQ=255)
                                                                                                                                             gCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTaa                                                                                                                                  <  1:36838/288‑1 (MQ=255)
                                                                                                                                               tCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAAc                                                                                                                                <  2:508287/288‑1 (MQ=255)
                                                                                                                                                   cTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCa                                                                                                                                                                                                                                                              <  1:687324/158‑1 (MQ=255)
                                                                                                                                                   cTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCa                                                                                                                                                                                                                                                              >  2:687324/1‑158 (MQ=255)
                                                                                                                                                             aaTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTaa                                                                                                                                                                                                  >  2:637150/1‑208 (MQ=255)
                                                                                                                                                             aaTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTaa                                                                                                                                                                                                  <  1:637150/208‑1 (MQ=255)
                                                                                                                                                              aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTaaaaa                                                                                                                                           >  1:296120/1‑262 (MQ=255)
                                                                                                                                                              aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTaaaaa                                                                                                                                           <  2:296120/262‑1 (MQ=255)
                                                                                                                                                                   attATTCAGTTCAAGTGGTGGCAGATACTGAGGCTCGTGGGTAGAAATGAGCCAGCCGTTTGTGGCAATAGGTACGAATACGCCAATTTTCATCATCAAGCTCGCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATACTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATa                                                                                                            <  2:332139/287‑1 (MQ=255)
                                                                                                                                                                            ttCAAAGGTCGGCATGTACTGCGGTGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTGTCATCATCACCCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATTCTGCATATCCTTTTCAGCCGCGAATTGGCTTGTTTGCAAAGCGGATGCCAGCTATTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTAGCTAAAATAAATGGTAAAAAGTGGACTAAACGGTCAAAACAGATGCACATAAAACATgaa                                                                                                   >  2:11558/1‑286 (MQ=255)
                                                                                                                                                                                 aGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGTTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCtgtg                                                                                             >  2:288788/1‑289 (MQ=255)
                                                                                                                                                                                 aGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAAAGGCTAAAACAGTTGCACATAAAACATGCATCtttg                                                                                             >  2:78538/1‑289 (MQ=255)
                                                                                                                                                                                      ggCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGt                                                                                                                                                                                                                                                                                    >  2:270971/1‑101 (MQ=255)
                                                                                                                                                                                      ggCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGt                                                                                                                                                                                                                                                                                    <  1:270971/101‑1 (MQ=255)
                                                                                                                                                                                                  ggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAg                                                                             <  1:663903/288‑1 (MQ=255)
                                                                                                                                                                                                                aaTGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATgaga                                                                                   <  2:463074/268‑1 (MQ=255)
                                                                                                                                                                                                                aaTGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATgaga                                                                                   >  1:463074/1‑268 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCcgggc                                                             <  1:141199/285‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCcgggc                                                             >  2:141199/1‑285 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCGTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCtg                                                                                               <  1:399868/251‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCGTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCtg                                                                                               <  1:399885/251‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCGTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCtg                                                                                               >  2:399885/1‑251 (MQ=255)
                                                                                                                                                                                                                     gccagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCGTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCtg                                                                                               >  2:399868/1‑251 (MQ=255)
                                                                                                                                                                                                                                        agttACGAATACTCCAATTTTCATCATCAACCTCTCTTCTTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTTTTGcc                                      <  1:11558/286‑1 (MQ=255)
                                                                                                                                                                                                                                           tACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTGTTGCCa                                     <  1:102039/287‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTGTTGCCAATCTACAagag                          >  1:166663/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                            ttCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCAGCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTGTTGCCAATCTACAAGAGGGGAGAGCGCATGACGCAAGgcgc  <  1:635325/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                
TGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGTTAATGTTATTAATCTGTTAACATTACGTTATCTAAAATATCTGGTAAAAAGTGGACTAAACGGTCAAAACAGTTGCACATAAAACATGCATCTGTGCGCGATGAGAGTGCAGAAGGTCGAGGCCGGGCGGGGGTTTTGCTATCCTGTTGCCAATCTACAAGAGGGGAGAGCGCATGACGCAAGGCGC  >  NC_000913/1073699‑1074255

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