Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,158,246 T→C G385G (GGT→GGC ptsG → PTS system, glucose‑specific IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,2460TC85.7% 667.2 / 91.5 280G385G (GGT→GGCptsGPTS system, glucose‑specific IIBC component
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (133/107);  total (152/128)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

TGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTAC  >  NC_000913/1158154‑1158338
                                                                                            |                                                                                            
  gCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCg                                                                                            >  1:1876092/1‑93 (MQ=255)
     gTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGa                                                                                         >  2:3031143/1‑93 (MQ=255)
     gTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGa                                                                                         >  1:283523/1‑93 (MQ=255)
      ttCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGaa                                                                                        <  1:643288/93‑1 (MQ=255)
      ttCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGaa                                                                                        <  1:1556968/93‑1 (MQ=255)
        cccGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGa                                                                                      <  1:6935132/93‑1 (MQ=255)
         ccGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGAc                                                                                     >  2:2588821/1‑93 (MQ=255)
         ccGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGAc                                                                                     >  2:1333395/1‑93 (MQ=255)
          cGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGAGAAGAcg                                                                                    >  2:6894020/1‑93 (MQ=255)
           gATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGAcgc                                                                                   <  1:7046871/93‑1 (MQ=255)
           gATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGAcgc                                                                                   <  1:7476330/93‑1 (MQ=255)
            aTCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGAcgcg                                                                                  >  2:5425081/1‑93 (MQ=255)
             tcgtcgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGa                                                                                 >  1:1689329/1‑93 (MQ=255)
             tcgtcgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGa                                                                                 >  2:1448207/1‑93 (MQ=255)
             tcgtcgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGa                                                                                 >  1:913426/1‑93 (MQ=255)
             tcgtcgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTAATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGa                                                                                 >  1:7301246/1‑93 (MQ=255)
               gtcgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACt                                                                               >  2:386507/1‑93 (MQ=255)
                 cgGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGa                                                                             >  2:5164121/1‑93 (MQ=255)
                   gTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAg                                                                           >  1:6311865/1‑93 (MQ=255)
                    tATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGa                                                                          <  1:5558020/93‑1 (MQ=255)
                    tATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGa                                                                          <  2:2772442/93‑1 (MQ=255)
                    tATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGa                                                                          >  1:4444441/1‑93 (MQ=255)
                    tATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGa                                                                          <  2:4986518/93‑1 (MQ=255)
                     aTCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGAt                                                                         >  2:1613148/1‑93 (MQ=255)
                     aTCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGAt                                                                         <  1:4218488/93‑1 (MQ=255)
                     aTCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGAt                                                                         <  2:7503092/93‑1 (MQ=255)
                     aTCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGAt                                                                         >  1:1103265/1‑93 (MQ=255)
                      tCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATg                                                                        >  2:2038941/1‑93 (MQ=255)
                      tCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATg                                                                        >  1:3413213/1‑93 (MQ=255)
                       cGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGc                                                                       <  1:1197249/93‑1 (MQ=255)
                        ggTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCa                                                                      >  1:5047064/1‑93 (MQ=255)
                        ggTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCa                                                                      >  2:7170505/1‑93 (MQ=255)
                         gTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCaa                                                                     >  2:3310223/1‑93 (MQ=255)
                         gTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCaa                                                                     >  2:591661/1‑93 (MQ=255)
                          ttATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCaaa                                                                    <  2:2212767/93‑1 (MQ=255)
                           tATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCaaaa                                                                   >  2:913184/1‑93 (MQ=255)
                            aTGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAg                                                                  <  1:6695630/93‑1 (MQ=255)
                             tGCGATTGTTTACTACACCATCTTCCGCGTGCTGCTTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGc                                                                 >  1:1048632/1‑93 (MQ=255)
                             tGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGc                                                                 <  1:6822424/93‑1 (MQ=255)
                             tGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGc                                                                 >  2:6484441/1‑93 (MQ=255)
                                gATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGAc                                                              >  1:5686461/1‑93 (MQ=255)
                                 aTTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACa                                                             >  1:4788462/1‑93 (MQ=255)
                                  ttGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAg                                                            >  1:4179100/1‑93 (MQ=255)
                                  ttGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAg                                                            >  2:510444/1‑93 (MQ=255)
                                   tGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAgg                                                           <  1:29595/93‑1 (MQ=255)
                                    gTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGt                                                          >  2:1957577/1‑93 (MQ=255)
                                    gTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGt                                                          >  1:6784165/1‑93 (MQ=255)
                                    gTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGt                                                          <  1:5304846/93‑1 (MQ=255)
                                     tttACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTa                                                         >  1:3614331/1‑93 (MQ=255)
                                      ttACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTAc                                                        <  1:5504716/93‑1 (MQ=255)
                                       tactacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTAcc                                                       <  1:966472/93‑1 (MQ=255)
                                        actacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCa                                                      >  1:5869079/1‑93 (MQ=255)
                                        actacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCa                                                      <  2:6280636/93‑1 (MQ=255)
                                        actacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCa                                                      >  1:2202563/1‑93 (MQ=255)
                                         ctacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAg                                                     <  2:1482615/93‑1 (MQ=255)
                                          tacACCATCTTCCGCGTGCTGATTAAAGCACTGGGTCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGc                                                    >  2:2400913/1‑93 (MQ=255)
                                          tacACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGc                                                    >  2:3101908/1‑93 (MQ=255)
                                              ccATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACGGGTACCAGCGaaa                                                <  2:4354683/93‑1 (MQ=255)
                                              ccATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGaaa                                                <  2:1092617/93‑1 (MQ=255)
                                              ccATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGaaa                                                <  2:4876625/93‑1 (MQ=255)
                                              ccATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGaaa                                                <  1:5642227/93‑1 (MQ=255)
                                              ccATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGaaa                                                <  2:6944242/93‑1 (MQ=255)
                                               cATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAAt                                               <  2:4084136/93‑1 (MQ=255)
                                                   ttCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCa                                           <  1:6152749/93‑1 (MQ=255)
                                                   ttCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCa                                           >  1:3960943/1‑93 (MQ=255)
                                                    tCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCAc                                          <  1:3159757/93‑1 (MQ=255)
                                                    tCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCAc                                          <  2:4721193/93‑1 (MQ=255)
                                                       gcgTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCgg                                       >  2:5203092/1‑93 (MQ=255)
                                                       gcgTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCgg                                       >  1:1835758/1‑93 (MQ=255)
                                                          tGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGctc                                    >  1:2112669/1‑93 (MQ=255)
                                                          tGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGctc                                    <  1:3028136/93‑1 (MQ=255)
                                                          tGCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAATGGCACCGGctct                                   >  1:5838152/1‑93 (MQ=255)
                                                           gCTGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGctct                                   >  2:535083/1‑93 (MQ=255)
                                                             tGATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTgg                                 >  2:2880513/1‑93 (MQ=255)
                                                              gATTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGt                                >  1:3827531/1‑93 (MQ=255)
                                                               aTTAAAGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGtt                               >  1:2624838/1‑93 (MQ=255)
                                                                   aaGCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTtgctg                           >  1:4688472/1‑93 (MQ=255)
                                                                     gCACTGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCa                         >  2:5555735/1‑93 (MQ=255)
                                                                        cTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCAttt                      <  2:1079532/93‑1 (MQ=255)
                                                                         tGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtg                     <  2:6403245/93‑1 (MQ=255)
                                                                         tGGATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtg                     <  2:3710673/93‑1 (MQ=255)
                                                                          ggATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtgg                    <  2:6294336/93‑1 (MQ=255)
                                                                           gATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtggt                   <  2:3961118/93‑1 (MQ=255)
                                                                           gATCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtggt                   <  2:6757018/93‑1 (MQ=255)
                                                                            aTCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtggtg                  <  2:6131547/93‑1 (MQ=255)
                                                                             tCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtggtgg                 >  1:6004144/1‑93 (MQ=255)
                                                                             tCTGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTtggtgg                 >  1:4076343/1‑93 (MQ=255)
                                                                               tGAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTa               >  2:1738699/1‑93 (MQ=255)
                                                                                gAAAACGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTaa              <  2:6104563/93‑1 (MQ=255)
                                                                                  aaaCGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAg            >  2:3558806/1‑93 (MQ=255)
                                                                                  aaaCGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAg            >  2:7359291/1‑93 (MQ=255)
                                                                                    aCGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGaa          >  2:2643761/1‑93 (MQ=255)
                                                                                    aCGCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGaa          >  2:2999867/1‑93 (MQ=255)
                                                                                      gCCGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGaaaa        >  2:420325/1‑93 (MQ=255)
                                                                                       ccGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAAc       <  2:3246855/93‑1 (MQ=255)
                                                                                        cGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACa      <  2:4645549/93‑1 (MQ=255)
                                                                                        cGGGCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACa      <  2:7080534/93‑1 (MQ=255)
                                                                                         gggCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACAt     >  2:1995039/1‑93 (MQ=255)
                                                                                         gggCCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACAt     >  2:1570491/1‑93 (MQ=255)
                                                                                           gTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTa   <  1:4925685/93‑1 (MQ=255)
                                                                                            |                                                                                            
TGGCTGTTCCCGATCGTCGGTATCGGTTATGCGATTGTTTACTACACCATCTTCCGCGTGCTGATTAAAGCACTGGATCTGAAAACGCCGGGTCGTGAAGACGCGACTGAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTAC  >  NC_000913/1158154‑1158338

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