Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,528 T→C Q8R (CAG→CGG)  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,5280TC100.0% 140.5 / NA 42Q8R (CAG→CGG) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/24);  total (18/24)

GGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAA  >  NC_000913/66288‑66766
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
gggAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                                                                                                                                                                        >  1:640‑M1/1‑245 (MQ=255)
            tAGAGCCGCAGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACAGGTTTCGATGCTAACCACGACCATATCGTCAGCGGGCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTAGAGAACAACGGCGCTGACGTAGACCCATGAGAGCGGGAACAGGTTGTGGATTGGCAGCGCCAGGTGGGCTTCTAATACCCGGCGcccgagacccgcctgcgt                                                                                                                                                                                                                          >  2:228873‑M1/1‑233 (MQ=255)
             agagCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtc                                                                                                                                                                                                                         <  1:302720‑M1/251‑20 (MQ=255)
               agcagcCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcga                                                                                                                                                                                                                       <  2:184497‑M1/251‑22 (MQ=255)
                       gTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcccgcgtcgacctgcagg                                                                                                                                                                                                               >  2:170181‑M1/1‑222 (MQ=255)
                       gTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcagg                                                                                                                                                                                                               <  1:170181‑M1/251‑30 (MQ=255)
                                 gtgtCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACAAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgccagcgccgcactgcaggtatggggaa                                                                                                                                                                                                      >  2:424616‑M1/1‑212 (MQ=255)
                                 gtggcggagGAGGGCTTGTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTGTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGGGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaa                                                                                                                                                                                                     <  2:148080‑M1/247‑40 (MQ=255)
                                   gtCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgcgacccgcacgcgtcgacctgcaggtatggagaaac                                                                                                                                                                                                    >  1:325556‑M1/1‑210 (MQ=255)
                                      ggaggaGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAAACACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgaccagcaggtatggagaaacagt                                                                                                                                                                                                 >  1:396420‑M1/1‑207 (MQ=255)
                                       gagtagggcttTTTCGTACCTTCAACCACGTCACCGGTTTCGATGCTAACCACGACCATGTCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGGCGCGCTCGCGATCAACGGCGCTTACGTTGCGCCTTGTGAGCGTGACCGGTTTGTGGTTTGGCAGCGCCAGGTTGGCTTCTAAGACCCGGCGcccgagtccctaacgcgtcgacctgcaggtgtggagaatgagtag                                                                                                                                                                                               <  2:335954‑M1/247‑46 (MQ=255)
                                              gCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttg                                                                                                                                                                                        >  1:304956‑M1/1‑199 (MQ=255)
                                               cTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTACCCCCGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGCCCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgc                                                                                                                                                                                       >  1:148722‑M1/1‑198 (MQ=255)
                                                   gtcGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATAGCGTCAGCGGTCATGACGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGTCGTTGCCCCATGTGAGCGTGACCAGGTTGTGGTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgat                                                                                                                                                                                    <  1:149425‑M1/249‑57 (MQ=255)
                                                     cGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataa                                                                                                                                                                                  >  1:319562‑M1/1‑192 (MQ=255)
                                                       ttccttCAACCACTTCACCGGGTTCGATGCTAACCTCGACCATATCGTCAGCTGTCATGACGCTGTAATCGACGCCGGAAGGTATGATCACAAAGACGTCGCGCTCGCGATCAACGTCGCGGTCGTTGCCGCATGTGAACGGGTCCAGGTTGTGTTTTGGCAGTGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcggcgacctgcaggtatggagaaacagtagagagttgcgataaa                                                                                                                                                                                 <  1:283603‑M1/247‑60 (MQ=255)
                                                                     tCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                                                                                                                                                                  <  2:370313‑M1/250‑75 (MQ=255)
                                                                             ttCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCTCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagactcgcacgcgtcgacctgcaggtatggagaaacagtagagaggtgcgataaaaagcttcaggtaggatccgctaa                                                                                                                                                          <  2:196466‑M1/250‑83 (MQ=255)
                                                                              tCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaat                                                                                                                                                         <  1:120021‑M1/250‑84 (MQ=255)
                                                                                      aaCCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACTGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagtggcgataaaaagcgtcaggtaggatccgctaatcttatgga                                                                                                                                                 <  2:396420‑M1/250‑92 (MQ=255)
                                                                                        ccACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataa                                                                                                                                              <  2:337556‑M1/251‑95 (MQ=255)
                                                                                              ccATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgc                                                                                                                                        <  2:8115‑M1/251‑101 (MQ=255)
                                                                                                       cAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcata                                                                                                                               <  1:330984‑M1/251‑110 (MQ=255)
                                                                                                            gTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaa                                                                                                                           >  2:118078‑M1/1‑137 (MQ=255)
                                                                                                                         aaTCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgaagccgtg                                                                                                             >  2:42596‑M1/1‑124 (MQ=255)
                                                                                                                          aTCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTGTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgggc                                                                                                            <  2:361634‑M1/251‑129 (MQ=255)
                                                                                                                             gACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaacagcgtacggtaggatccgctaatcttatggagaaaaatgctatggcatagcaaagtgtgccgcagtgcaa                                                                                                          >  2:320172‑M1/1‑120 (MQ=255)
                                                                                                                                                                cgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattataga                                                                       <  2:274688‑M1/250‑166 (MQ=255)
                                                                                                                                                                cgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattataga                                                                       <  2:68175‑M1/250‑166 (MQ=255)
                                                                                                                                                                 gcgATCCACGGAGCTGACGTTGCTCCATGTGAGCGTGACCAGGTTGTGTTTTGGCCGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagac                                                                      >  2:253048‑M1/1‑84 (MQ=255)
                                                                                                                                                                    agcaaCGGCGCTGACGTTGCCCCATGTGAGCGTGACTAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCTGCGcccgagacccgcacgcgtcggcctgcaggtatggagaaacagtagagaggtgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatgtcatagcaaagtgtgacgccgggcaaataatcgatgtggacttgtctgccgtgattatagacactt                                                                  <  1:335496‑M1/249‑171 (MQ=255)
                                                                                                                                                                            cgcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttaca                                                          <  1:118078‑M1/251‑179 (MQ=255)
                                                                                                                                                                                               gAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttg                                       >  1:315251‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                       ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgct                               >  1:162120‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                       ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgaccggcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgct                               <  1:233001‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                        cAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctt                              <  2:640‑M1/251‑207 (MQ=255)
                                                                                                                                                                                                         aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgcttt                             <  2:291298‑M1/251‑208 (MQ=255)
                                                                                                                                                                                                              gtgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgtta                         <  2:319562‑M1/250‑212 (MQ=255)
                                                                                                                                                                                                                   ttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaacgtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaat                   >  2:416‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                         gcgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaatgctttt             <  2:304956‑M1/251‑224 (MQ=255)
                                                                                                                                                                                                                                gTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgccctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatgtctttggtcccgctttgttacagaatgcttttaatacgc      >  2:149836‑M1/1‑21 (MQ=255)
                                                                                                                                                                                                                                    gCTTCTAATACCCGGCGcccgagacccgcaagcgtcgacctgcaggtatggagaaacagtagcgagttgcgataaaaagcgacaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgaagccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaatgcttttactaagcgggg  >  2:344361‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
GGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAA  >  NC_000913/66288‑66766

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