Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,463,185 A→G G190G (GGT→GGC treB ← PTS trehalose‑specific system enzyme IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,463,1850AG85.7% 713.4 / 99.4 272G190G (GGT→GGCtreBPTS trehalose‑specific system enzyme IIBC
Reads supporting (aligned to +/- strand):  ref base A (21/18);  new base G (100/133);  total (121/151)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 275 total aligned reads displayed.

ATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCA  >  NC_000913/4463096‑4463274
                                                                                         |                                                                                         
aTGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCa                                                                                        <  2:1035773/93‑1 (MQ=255)
 tGCCAAAGTCCCACACTTCCGGTAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCaa                                                                                       <  1:350921/93‑1 (MQ=255)
 tGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCaa                                                                                       >  2:4348909/1‑93 (MQ=255)
 tGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCaa                                                                                       >  2:3872330/1‑93 (MQ=255)
  gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAg                                                                                      <  2:5410769/93‑1 (MQ=255)
  gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAg                                                                                      >  1:3870521/1‑93 (MQ=255)
  gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAg                                                                                      <  2:1116083/93‑1 (MQ=255)
  gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAg                                                                                      <  2:647797/93‑1 (MQ=255)
   ccAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAgg                                                                                     >  1:2130109/1‑93 (MQ=255)
   ccAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAgg                                                                                     <  1:4741192/93‑1 (MQ=255)
    cAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGa                                                                                    >  1:85245/1‑93 (MQ=255)
     aaaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGAt                                                                                   <  1:4453121/93‑1 (MQ=255)
      aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATc                                                                                  <  2:5185386/93‑1 (MQ=255)
      aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATc                                                                                  >  1:4254606/1‑93 (MQ=255)
      aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATc                                                                                  >  2:5546415/1‑93 (MQ=255)
        ctCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCgg                                                                                <  1:194737/92‑1 (MQ=255)
         tCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGc                                                                               <  1:6456633/93‑1 (MQ=255)
          cccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCg                                                                              <  1:2887413/93‑1 (MQ=255)
           ccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGt                                                                             >  1:746603/1‑93 (MQ=255)
           ccACACGTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGt                                                                             <  2:1423601/93‑1 (MQ=255)
             acacTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTgc                                                                           >  2:5926852/1‑93 (MQ=255)
              cacTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTgcc                                                                          <  1:6347343/93‑1 (MQ=255)
               acTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTgccg                                                                         <  2:463586/93‑1 (MQ=255)
                 ttCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTgccgcc                                                                       <  1:3604101/93‑1 (MQ=255)
                   ccGGTAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCa                                                                     <  2:3506118/93‑1 (MQ=255)
                    cGGTAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAt                                                                    >  2:3289044/1‑93 (MQ=255)
                    cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAt                                                                    <  1:2652934/93‑1 (MQ=255)
                    cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAt                                                                    <  1:2002163/93‑1 (MQ=255)
                    cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAt                                                                    <  2:1006759/93‑1 (MQ=255)
                     ggCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAtt                                                                   <  2:2913317/93‑1 (MQ=255)
                      gCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAttt                                                                  >  2:2354099/1‑93 (MQ=255)
                       cAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAtttt                                                                 <  2:2005071/93‑1 (MQ=255)
                        aGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAttttt                                                                >  2:4167530/1‑93 (MQ=255)
                         gctgctGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAtttttt                                                               >  2:1823497/1‑93 (MQ=255)
                         gctgctGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCAtttttt                                                               <  2:6166749/93‑1 (MQ=255)
                         gctgctGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtttttt                                                               <  2:3035349/93‑1 (MQ=255)
                            gctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCa                                                            >  1:3551199/1‑93 (MQ=255)
                             ctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCAc                                                           >  2:426992/1‑93 (MQ=255)
                                cccGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGc                                                        <  2:589417/93‑1 (MQ=255)
                                  cGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTg                                                      >  2:3066205/1‑93 (MQ=255)
                                  cGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTg                                                      <  1:1131832/93‑1 (MQ=255)
                                  cGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTg                                                      <  1:125392/93‑1 (MQ=255)
                                  cGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTg                                                      <  1:5474688/93‑1 (MQ=255)
                                   gAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                     >  1:3435309/1‑93 (MQ=255)
                                     gCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGAcc                                                   <  2:451626/93‑1 (MQ=255)
                                     gCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGAcc                                                   <  2:2440581/93‑1 (MQ=255)
                                     gCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGAcc                                                   <  2:5152912/93‑1 (MQ=255)
                                     gCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGAcc                                                   >  2:4006862/1‑93 (MQ=255)
                                       aGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAg                                                 >  1:3093002/1‑93 (MQ=255)
                                        gATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGc                                                <  1:5736442/93‑1 (MQ=255)
                                        gATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGc                                                >  1:479151/1‑93 (MQ=255)
                                          tAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCaa                                              >  1:4903427/1‑93 (MQ=255)
                                           aaGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCaaa                                             <  1:3212379/93‑1 (MQ=255)
                                            aGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAAt                                            <  1:5796776/93‑1 (MQ=255)
                                            aGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAAt                                            >  1:4324890/1‑93 (MQ=255)
                                              cGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATac                                          >  2:5653497/1‑93 (MQ=255)
                                               gTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATacc                                         >  1:613905/1‑93 (MQ=255)
                                               gTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACACCTGACCAGCAAAAacc                                         >  2:3707896/1‑93 (MQ=255)
                                                ttCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccg                                        >  1:40277/1‑93 (MQ=255)
                                                 tcatcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccga                                       >  2:699488/1‑93 (MQ=255)
                                                  tatcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgac                                      <  2:5501894/92‑1 (MQ=255)
                                                  catcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgac                                      >  2:3171920/1‑93 (MQ=255)
                                                  catcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgac                                      <  2:173378/93‑1 (MQ=255)
                                                  catcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgac                                      <  2:265791/93‑1 (MQ=255)
                                                  catcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgac                                      <  1:5382170/93‑1 (MQ=255)
                                                   atcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgacc                                     <  2:6510676/93‑1 (MQ=255)
                                                   atcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgacc                                     >  1:1552698/1‑93 (MQ=255)
                                                       gCTGTGGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCa                                 >  2:857022/1‑93 (MQ=255)
                                                         tgtgGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                               >  1:5887711/1‑93 (MQ=255)
                                                          gtgGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGt                              <  1:4067019/93‑1 (MQ=255)
                                                          gtgGAGAAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGt                              <  2:2513821/93‑1 (MQ=255)
                                                             gagaAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTaga                           <  2:3865833/93‑1 (MQ=255)
                                                              agaAACCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaa                          <  2:3392835/93‑1 (MQ=255)
                                                              agaAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaa                          <  1:1463986/93‑1 (MQ=255)
                                                               gaAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaag                         >  1:3026207/1‑93 (MQ=255)
                                                                aaaCCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaaga                        <  2:5608566/93‑1 (MQ=255)
                                                                aaaCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaaga                        <  1:1196750/93‑1 (MQ=255)
                                                                 aaCCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaa                       <  1:253622/93‑1 (MQ=255)
                                                                  aCCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaaa                      >  1:1767010/1‑92 (MQ=255)
                                                                  aCCAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGGagaagaag                      >  2:156346/1‑93 (MQ=255)
                                                                   ccAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaaga                     <  1:2439273/93‑1 (MQ=255)
                                                                   ccAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaaga                     >  2:4491805/1‑93 (MQ=255)
                                                                    cAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGAt                    <  2:3219430/93‑1 (MQ=255)
                                                                    cAGTGTCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGAt                    <  2:6470158/93‑1 (MQ=255)
                                                                        gtCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCt                <  1:231744/93‑1 (MQ=255)
                                                                         tCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCtt               <  2:1644147/93‑1 (MQ=255)
                                                                         tCACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCtt               <  1:1423493/93‑1 (MQ=255)
                                                                          cACGCCAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTc              <  2:2428465/93‑1 (MQ=255)
                                                                              ccAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCg          >  2:4407351/1‑93 (MQ=255)
                                                                              ccAAGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCg          <  2:2134821/93‑1 (MQ=255)
                                                                              ccAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCg          >  2:6612762/1‑93 (MQ=255)
                                                                                 aGCACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATc       <  1:1054349/93‑1 (MQ=255)
                                                                                 aGCACGATGCCAAGGATCGGAGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATc       >  2:6372467/1‑93 (MQ=255)
                                                                                 aGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATc       <  2:750868/93‑1 (MQ=255)
                                                                                   cACGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCaa     <  2:934724/93‑1 (MQ=255)
                                                                                    aCGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAAc    <  2:2325961/93‑1 (MQ=255)
                                                                                     cGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAGGATCGCTTCACCGATCAAcc   <  1:3708126/93‑1 (MQ=255)
                                                                                     cGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAAcc   <  1:5067303/93‑1 (MQ=255)
                                                                                     cGATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAAcc   >  2:3009964/1‑93 (MQ=255)
                                                                                      gATGCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACca  <  2:5952364/93‑1 (MQ=255)
                                                                                         |                                                                                         
ATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCA  >  NC_000913/4463096‑4463274

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