Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,564,531 T→C H684R (CAT→CGT)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,564,5310TC100.0% 1482.8 / NA 352H684R (CAT→CGT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (205/147);  total (205/147)

Only 100 of 355 total aligned reads displayed.

CTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCC  >  NC_000913/3564433‑3564628
                                                                                                  |                                                                                                 
cTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGa                                                                                                 >  1:6874940/1‑101 (MQ=255)
 tttCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGAt                                                                                                >  1:6874943/1‑101 (MQ=255)
  ttCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATc                                                                                               <  1:6874390/101‑1 (MQ=255)
  ttCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATc                                                                                               >  2:6874949/1‑101 (MQ=255)
    ctcGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCa                                                                                             >  1:6874950/1‑101 (MQ=255)
      cGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGc                                                                                           <  1:6874333/101‑1 (MQ=255)
       gtagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCa                                                                                          >  1:6874956/1‑101 (MQ=255)
         agtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCAtc                                                                                        >  1:1790423/1‑101 (MQ=255)
         agtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCAtc                                                                                        >  2:6874958/1‑101 (MQ=255)
           taTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCAtctc                                                                                      <  2:6874303/101‑1 (MQ=255)
           taTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCAtctc                                                                                      >  1:6874964/1‑101 (MQ=255)
           taTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCAtctc                                                                                      >  2:6874965/1‑101 (MQ=255)
             ttCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGa                                                                                    >  1:6874970/1‑101 (MQ=255)
              tCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGAc                                                                                   >  2:6874972/1‑101 (MQ=255)
               cACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACa                                                                                  <  1:6874317/101‑1 (MQ=255)
                aCGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACAt                                                                                 >  2:6874977/1‑101 (MQ=255)
                aCGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACAt                                                                                 <  2:6874325/101‑1 (MQ=255)
                 cgcgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACAtt                                                                                >  2:6874981/1‑101 (MQ=255)
                   cgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCg                                                                              <  1:6874720/101‑1 (MQ=255)
                   cgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCg                                                                              <  1:6874637/101‑1 (MQ=255)
                      tttGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCAc                                                                           <  1:6874421/101‑1 (MQ=255)
                       ttGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACAt                                                                                 <  1:6874241/94‑1 (MQ=255)
                       ttGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCAcc                                                                          <  2:6874571/101‑1 (MQ=255)
                       ttGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCAcc                                                                          >  2:6874984/1‑101 (MQ=255)
                         gTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGt                                                                        >  1:6874987/1‑101 (MQ=255)
                            gCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCa                                                                     <  2:6874605/101‑1 (MQ=255)
                              ccTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAAc                                                                   >  1:6874992/1‑101 (MQ=255)
                              ccTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAAc                                                                   <  2:6874381/101‑1 (MQ=255)
                               cTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACg                                                                  >  2:6874993/1‑101 (MQ=255)
                                tGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGt                                                                 >  2:6874996/1‑101 (MQ=255)
                                   cgacgCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTAcc                                                              <  2:6874393/101‑1 (MQ=255)
                                    gacgCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCg                                                             >  1:6875001/1‑101 (MQ=255)
                                     acgCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGa                                                            >  2:6875005/1‑101 (MQ=255)
                                       gCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATa                                                          >  1:6875007/1‑101 (MQ=255)
                                       gCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATa                                                          <  1:6874502/101‑1 (MQ=255)
                                         aGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGt                                                        >  2:6875013/1‑101 (MQ=255)
                                          gTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTc                                                       <  1:6874356/101‑1 (MQ=255)
                                          gTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTc                                                       <  1:6874581/101‑1 (MQ=255)
                                           ttcttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCa                                                      >  1:6875016/1‑101 (MQ=255)
                                            tcttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAg                                                     <  2:6874515/101‑1 (MQ=255)
                                             cttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgc                                                    <  2:6874536/101‑1 (MQ=255)
                                              ttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcg                                                   <  2:6874188/101‑1 (MQ=255)
                                                ccACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCa                                                 >  2:6875020/1‑101 (MQ=255)
                                                ccACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCa                                                 <  1:6874537/101‑1 (MQ=255)
                                                ccACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCa                                                 >  1:6875022/1‑101 (MQ=255)
                                                  aCTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAcc                                               <  2:6874810/101‑1 (MQ=255)
                                                  aCTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCAcc                                               >  1:6875025/1‑101 (MQ=255)
                                                   cttcttCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCg                                              >  1:6875029/1‑101 (MQ=255)
                                                    ttcttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGt                                             <  1:6874626/101‑1 (MQ=255)
                                                     tcttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGtt                                            <  2:6874495/101‑1 (MQ=255)
                                                      cttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTa                                           >  2:6875033/1‑101 (MQ=255)
                                                       ttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa                                          >  2:6875037/1‑101 (MQ=255)
                                                        tcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGt                                             >  1:6875039/1‑97 (MQ=255)
                                                        tcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAg                                         <  1:6874445/101‑1 (MQ=255)
                                                           gCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgcgc                                      >  2:6875045/1‑101 (MQ=255)
                                                             tgtgTTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaa                                    >  2:6875048/1‑101 (MQ=255)
                                                              gtgtTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaaa                                   <  2:6874464/101‑1 (MQ=255)
                                                              gtgtTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaaa                                   >  2:6875053/1‑101 (MQ=255)
                                                               tgtTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCg                                                             >  1:6875056/1‑74 (MQ=255)
                                                                gtTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgc                                        >  1:6875060/1‑94 (MQ=255)
                                                                gtTACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACt                                 >  2:6875058/1‑101 (MQ=255)
                                                                 ttACCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACtt                                <  2:6874413/101‑1 (MQ=255)
                                                                   aCCAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCa                              >  2:6875066/1‑101 (MQ=255)
                                                                    ccAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCAt                             >  2:6875069/1‑101 (MQ=255)
                                                                    ccAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCAt                             >  1:6875073/1‑101 (MQ=255)
                                                                     cAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATg                            >  2:6875075/1‑101 (MQ=255)
                                                                     cAAAAATAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATg                            <  1:6874567/101‑1 (MQ=255)
                                                                      aaaaaTAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGt                           <  2:6874589/101‑1 (MQ=255)
                                                                       aaaaTAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGtt                          >  1:6875079/1‑101 (MQ=255)
                                                                         aaTAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCAt                             >  2:6875083/1‑96 (MQ=255)
                                                                          aTAAAGATATTGTCAGTACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACt                       >  1:6875086/1‑101 (MQ=255)
                                                                          aTAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACt                       >  2:6875090/1‑101 (MQ=255)
                                                                           tAAAGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTg                      >  1:6875093/1‑101 (MQ=255)
                                                                            aaaGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTgg                     <  2:6874152/101‑1 (MQ=255)
                                                                             aaGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGt                    >  1:6875097/1‑101 (MQ=255)
                                                                             aaGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGt                    <  1:6874519/101‑1 (MQ=255)
                                                                              aGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTg                   >  2:6875101/1‑101 (MQ=255)
                                                                              aGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTg                   <  2:6874352/101‑1 (MQ=255)
                                                                              aGATATTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTg                   <  2:6874503/101‑1 (MQ=255)
                                                                                atatTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGcc                 >  2:6875106/1‑101 (MQ=255)
                                                                                  atTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCgg               <  2:6874627/101‑1 (MQ=255)
                                                                                  atTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCgg               >  1:6875112/1‑101 (MQ=255)
                                                                                  atTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCgg               >  1:6875108/1‑101 (MQ=255)
                                                                                  atTGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCgg               <  1:6874765/101‑1 (MQ=255)
                                                                                   ttGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGa              <  2:6874424/101‑1 (MQ=255)
                                                                                    tGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGaa             >  2:6875117/1‑101 (MQ=255)
                                                                                    tGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGaa             >  1:6875118/1‑101 (MQ=255)
                                                                                    tGTCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGaa             <  1:6874695/101‑1 (MQ=255)
                                                                                      tCAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACATACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGc           >  2:6875122/1‑101 (MQ=255)
                                                                                       cAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCt          >  2:6875129/1‑101 (MQ=255)
                                                                                       cAGCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCt          >  1:6875126/1‑101 (MQ=255)
                                                                                         gCACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTc        >  1:6875134/1‑101 (MQ=255)
                                                                                          cACCGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTcc       >  2:6875135/1‑101 (MQ=255)
                                                                                            ccGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGtt                          <  2:6874484/80‑1 (MQ=255)
                                                                                            ccGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGt     >  2:6875140/1‑101 (MQ=255)
                                                                                            ccGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGt     <  2:6874697/101‑1 (MQ=255)
                                                                                            ccGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGt     <  1:6874092/101‑1 (MQ=255)
                                                                                             cGACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTc    >  1:6875143/1‑101 (MQ=255)
                                                                                              gACACGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTcc   <  2:6874453/101‑1 (MQ=255)
                                                                                               acacGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTccc  <  1:6874276/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCC  >  NC_000913/3564433‑3564628

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