Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,518 T→C E11E (GAA→GAG araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,5180TC92.2% 150.1 / ‑1.7 52E11E (GAA→GAGaraDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (23/24);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGG  >  NC_000913/66275‑66758
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
gccgccAATGGACGGGAATGCCTGATAGAGCAGTAGGTGAGTTTGCGGATCGGAGGAGGGCGTTTTCGCACCTTCAACCACTTCACCGGGTGCGAGGCTAACCGCGACCATATCGTCATCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACTAGGTTATGTTTTGGCAGCGCCATGTGGGCTTCTAAtt                                                                                                                                                                                                                                           <  2:426025/251‑2 (MQ=255)
    ccAAAGGACGGGAATGCATGATTGAGCAGTCTGTGAGGTGGCGTATCGGGTGTGGGCGTTTTCGCACCTTCAACCACGTCACCGGGTTCGATGCTAACCACGACCATCGCGTCAGCGGTCATGATTCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGTGCCGACGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCTCCAGGTTGGCCTCTAATatct                                                                                                                                                                                                                                        <  2:235404/250‑4 (MQ=255)
      aaTGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCGCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGc                                                                                                                                                                                                                                     >  1:109344/1‑251 (MQ=255)
        tGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGt                                                                                                                                                                                                                                   >  1:463043/1‑251 (MQ=255)
           acgggAATGCCTGATAGAGCCGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGCCCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTCATACCTGGCGTTTg                                                                                                                                                                                                                                >  2:445770/3‑251 (MQ=255)
           acgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                                                                                                                                                                                                                <  2:310626/249‑1 (MQ=255)
             gggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                                                                                                                                                                                                                >  1:282904/1‑249 (MQ=255)
                       gATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGGCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCGTGTGAGCGTGACCAGGTTAGGTTGTGGCAGCGCCAGGTTGGCTTCTAATACCTGGTGTTTGAGATCTTCTa                                                                                                                                                                                                                      <  1:407734/249‑1 (MQ=255)
                       gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCATGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTaa                                                                                                                                                                                                                     <  2:418595/250‑1 (MQ=255)
                          agagCAGTCGGTGAGTTGGCGTATCGGAGTAGGGCTCTTTCGCACCTTCAACCACTTCCCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGCCGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCATCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACAt                                                                                                                                                                                                                  <  2:129173/250‑1 (MQ=255)
                            agcagtCGGTGAGGTGGCGTATCGGAGGCGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGACAGCGGTCATGATGCTGTACTCGACGCCGGAAGGTTTGATAACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTGATAACTGACGTTTGAGATCTTCTAAAATGtt                                                                                                                                                                                                               >  2:399782/1‑251 (MQ=255)
                                   ggTGAGTTGGCGTAGCGGAGGAGTGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCt                                                                                                                                                                                                        <  2:199514/251‑1 (MQ=255)
                                       aGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAATGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATCCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGt                                                                                                                                                                                                    >  1:244836/1‑251 (MQ=255)
                                               tatCGGAGGCGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTACCCACGACAATATCGTCAGCGGTCATGCTGCTGTAATCGTCGCAGTAAGGTTTGATCACAAAGTAGACGCGCTCGCGATCAACGGCACTGACGTTGCCCCATGCGAGCGCGACCAGGTTATATTTTGGCAGCGCCAGGTTGACCGCTAATACCTTGGTTTTGGGATCTTCTACAATGTTAACTACTGCGTGACGGagt                                                                                                                                                                                            >  2:334826/3‑249 (MQ=255)
                                                    gaggagGTCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATAGCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTGGCCCCATGTGAGCGTGACCAGTTTATGTTTTGGCATCGCCAGGTTGGCCTCTAATGCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCtt                                                                                                                                                                                       <  1:440028/251‑1 (MQ=255)
                                                       gaggGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCCCGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGCCGTTGCCCCATGTGAGCGAGACCAGGTTATGTCTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTATCTCCTTCTTGCCGGTTGCGCTTTGc                                                                                                                                                                                    >  2:282039/1‑251 (MQ=255)
                                                        aggGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCt                                                                                                                                                                                   <  1:400719/251‑1 (MQ=255)
                                                                 tcgcaCCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATAGACGCCGGAAGGATTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGcc                                                                                                                                                                          >  1:443901/5‑251 (MQ=255)
                                                                             ccACTTCACCGGTTTCGATGCTAACCCCGACCATATCGGCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGTCTACAAAATCGCa                                                                                                                                                              <  2:109344/251‑1 (MQ=255)
                                                                             acACTTCCCCGGTTTCGATGCTAACAACGACCATATCGTCAGCGGTCATGATGCTGGAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGGGACCGGGTTATGTTTTGGCAGCTCCAGGTTGGCCTCTAATACCTGGCGTTTGCGATCTCCTAACATGTTGACTCCTTCGTGCCGGATGCGCATTGCTTATCCGGCCTACAAATTCGc                                                                                                                                                               >  2:328971/2‑250 (MQ=255)
                                                                                  taaCCGGTTGCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGCGCGTGACCCGGTAATGTTTTGGCAGCGCCAGGCTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAAACTGTTGACTCCTTCGTCCCGGATGCGCCTTGCTTATCCGGCCTACATCATCGCAGCgtg                                                                                                                                                         >  2:237211/3‑251 (MQ=255)
                                                                                             gATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCTgggg                                                                                                                                              >  2:119512/1‑248 (MQ=255)
                                                                                             gATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTg                                                                                                                                                 <  1:119512/248‑1 (MQ=255)
                                                                                                    aCCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCATCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAcgcg                                                                                                                                       <  1:90375/251‑1 (MQ=255)
                                                                                                        cGACCATATCGGCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACgcgctag                                                                                                                                   <  2:282904/251‑4 (MQ=255)
                                                                                                          aCCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAATGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCg                                                                                                                                 <  2:385792/251‑1 (MQ=255)
                                                                                                                  gTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGGTAAGACGCGCCAGCGTaaagttg                                                                                                                         >  2:308699‑M1/1‑244 (MQ=255)
                                                                                                                           cATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGAGCTCGCGATCAACGGCGCTGACGATGCCCCATGTGAGCGAGACCAGGATATGTTTTGGCAGCGCCAGGTTGGCCTCTAATAACTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTAAAAAATCGACGCGGTTAGGCCTGATAAAACGCGCCAGCGTaaagtgggcactgtg                                                                                                                 >  1:324942‑M1/1‑235 (MQ=255)
                                                                                                                              gatgctgTAATCGACGCCGGAAGGTTTGCTCACAAAGACGCGGCGCTCGCGATCAACGGCGCTGACGTTGCCACATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagc                                                                                                             >  1:272939‑M1/4‑232 (MQ=255)
                                                                                                                                 gcagtaaTCGACGCCGGAAGGTTTGATCACAAAGGCGCCGCGCTCGCGAGCAACGGCGCTGACGTGGCCCCATGTGATCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATCCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCTAGCGTaaagttggcgctgttaagcca                                                                                                           <  2:260598‑M1/247‑22 (MQ=255)
                                                                                                                                    gtatacgacgcCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCTATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagat                                                                                                        <  2:429037‑M1/245‑25 (MQ=255)
                                                                                                                                    gTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagat                                                                                                        >  1:293026‑M1/1‑226 (MQ=255)
                                                                                                                                              ccGGAAGGATTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGAAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctatcg                                                                                              >  1:294845‑M1/1‑216 (MQ=255)
                                                                                                                                                   aGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGAGCTTTGCTTATCCg                                                                                        >  1:97731/1‑251 (MQ=255)
                                                                                                                                                    ggTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGCGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTGTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctc                                                                                       <  1:123860‑M1/251‑42 (MQ=255)
                                                                                                                                                    ggTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctc                                                                                       >  1:111821‑M1/1‑210 (MQ=255)
                                                                                                                                                    ggTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctc                                                                                       <  2:111821‑M1/251‑42 (MQ=255)
                                                                                                                                                    ggTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctc                                                                                       <  1:194296‑M1/251‑42 (MQ=255)
                                                                                                                                                         gATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttc                                                                                   >  1:292930‑M1/1‑205 (MQ=255)
                                                                                                                                                            cacaAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatg                                                                                <  1:420630‑M1/250‑49 (MQ=255)
                                                                                                                                                                    cgccgcGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatc                                                                       <  1:325874‑M1/251‑58 (MQ=255)
                                                                                                                                                                    cgccgcCCTCGCGATCACCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTTGCAGCGCCAGGTTGGCCTCTAATTCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatc                                                                       <  2:293026‑M1/251‑58 (MQ=255)
                                                                                                                                                                          gcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagct                                                                 <  1:421409‑M1/251‑64 (MQ=255)
                                                                                                                                                                             cgcgATTAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGATTCCTTCGTTCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgat                                                              <  1:54415‑M1/251‑67 (MQ=255)
                                                                                                                                                                              gcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatc                                                             >  1:38230‑M1/1‑184 (MQ=255)
                                                                                                                                                                                             gACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccag                                               <  1:155069‑M1/250‑82 (MQ=255)
                                                                                                                                                                                              aCGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccagcg                                             >  1:83949‑M1/1‑168 (MQ=255)
                                                                                                                                                                                                   gCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccagcgtgacc                                        <  1:308699‑M1/251‑89 (MQ=255)
                                                                                                                                                                                                       cATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTCTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttgtcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccagcgtgacctggg                                    >  2:429835‑M1/1‑159 (MQ=255)
                                                                                                                                                                                                                    ccAGGTTATGCTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagcccgatatcgctctcgctctctccttcatgccactgatcgtagctgatctctttgatttccagcgtgacctggtgacttgccagaat                       >  2:199667‑M1/1‑146 (MQ=255)
                                                                                                                                                                                                                                     aGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccagcgtgacctggtgacttgccagaatctgctgcatgatcccgg      <  2:110375‑M1/251‑123 (MQ=255)
                                                                                                                                                                                                                                         ccAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGGTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTaaagttggcgctgttaagccagatatcgctctcgatctctccttcatgccactgatcgtagctgatctctttgatttccagcgtgacctggtgacttgccagaatctgctgcatgatcccggcaat  <  2:157992‑M1/251‑127 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
GCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGG  >  NC_000913/66275‑66758

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