Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,902,722 G→A G168G (GGC→GGT bglF ← PTS system beta‑glucosides‑specific enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,902,7220GA98.9% 621.1 / ‑3.7 179G168G (GGC→GGTbglFPTS system beta‑glucosides‑specific enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (92/85);  total (93/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 188 total aligned reads displayed.

GCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTAC  >  NC_000913/3902633‑3902810
                                                                                         |                                                                                        
gcATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccaccg                                                                                       >  2:219731/1‑93 (MQ=255)
 cATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGa                                                                                      >  2:2431323/1‑93 (MQ=255)
   tCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAg                                                                                    >  2:3125154/1‑93 (MQ=255)
   tCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAg                                                                                    >  1:2676371/1‑93 (MQ=255)
    ccGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGc                                                                                   <  1:2258212/93‑1 (MQ=255)
     cGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCg                                                                                  >  2:1676551/1‑93 (MQ=255)
     cGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCg                                                                                  >  1:4025697/1‑93 (MQ=255)
      gCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGTATTACCACCGAAGCGt                                                                                 >  1:3121294/1‑93 (MQ=255)
      gCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGt                                                                                 >  1:371345/1‑93 (MQ=255)
      catttttGCCCGTGCTCGAAAGCAGTCAGAATGAATGGATGCCCTAACGCGCCACCAATCTCCATGGCAGTAAATGGATTACCACCGAAGCGt                                                                                 <  2:902742/91‑1 (MQ=255)
       cTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGtt                                                                                <  2:2622142/93‑1 (MQ=255)
       cTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGtt                                                                                <  1:2570148/93‑1 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                               <  1:933227/93‑1 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                               <  1:4162470/93‑1 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                               >  1:4310732/1‑93 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                               >  1:3431295/1‑93 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGctt                                                                               >  2:2047629/1‑93 (MQ=255)
        tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCCCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                               >  1:4713475/1‑93 (MQ=255)
            tGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTcccc                                                                           <  1:2706677/93‑1 (MQ=255)
             gCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCcg                                                                          <  1:510345/93‑1 (MQ=255)
               ccGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCcgcg                                                                        >  1:3395899/1‑93 (MQ=255)
                 tttCTCGAACGCAGTCAGAATTAATGCATGCACTAACGCTCCACCAATCTCCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGgt                                                                      <  2:4321991/92‑1 (MQ=255)
                  ttCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGgtg                                                                     >  1:1006215/1‑93 (MQ=255)
                   tctcGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGgtgt                                                                    >  2:2048184/1‑93 (MQ=255)
                   tctcGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGgtgt                                                                    <  2:1806297/93‑1 (MQ=255)
                      cGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATc                                                                 <  2:2008276/93‑1 (MQ=255)
                      cGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATc                                                                 >  1:443079/1‑93 (MQ=255)
                       gAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATcc                                                                <  1:2344290/93‑1 (MQ=255)
                       gAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATcc                                                                <  1:3933658/93‑1 (MQ=255)
                        aaaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATccc                                                               <  2:1441374/93‑1 (MQ=255)
                          aGCAGTCAGAATTAATGGATGCACTACCGCTCCACCAATCACCATGGCAGCAAATGGGTTACCACCGAAGCGTTTCACCGCGGGGTCTCTCAg                                                             >  2:984840/1‑93 (MQ=255)
                          aGCAGGCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAg                                                             >  1:4441155/1‑93 (MQ=255)
                           gCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAgg                                                            <  1:4377414/93‑1 (MQ=255)
                            cAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAggt                                                           >  1:459352/1‑92 (MQ=255)
                            cAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGa                                                           <  1:477582/93‑1 (MQ=255)
                            cAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGa                                                           >  1:459353/1‑93 (MQ=255)
                             aGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGAt                                                          <  1:4101794/93‑1 (MQ=255)
                               tCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTa                                                        >  2:1409384/1‑93 (MQ=255)
                                cAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTAt                                                       <  1:3455199/93‑1 (MQ=255)
                                 aGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc                                                      >  2:4001633/1‑93 (MQ=255)
                                 aGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc                                                      >  1:3729019/1‑93 (MQ=255)
                                 aGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc                                                      <  1:136004/93‑1 (MQ=255)
                                 aGAATTAACGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc                                                      <  2:741866/93‑1 (MQ=255)
                                  gAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCg                                                     >  2:437765/1‑93 (MQ=255)
                                  gAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCg                                                     >  2:1009183/1‑93 (MQ=255)
                                   aaTTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCgg                                                    >  1:4651427/1‑93 (MQ=255)
                                    aTTAATGAATGCACTAACTCGCCTCCAATCTCCATGGCAGTAAATTGATTACCACCGAAGCGTTTTCCCGCGGTGTATCCTAGGATTTTCtgg                                                   <  2:3325916/93‑1 (MQ=255)
                                         tGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaagaa                                              >  2:1730657/1‑93 (MQ=255)
                                         tGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaagaa                                              <  2:3875666/93‑1 (MQ=255)
                                            aTGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCa                                           >  1:669944/1‑93 (MQ=255)
                                             tGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCaa                                          >  1:3330250/1‑93 (MQ=255)
                                              gCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCaaa                                         >  2:4103173/1‑93 (MQ=255)
                                              gCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCaaa                                         >  2:2760167/1‑93 (MQ=255)
                                                 cTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGCACCAAAACa                                      >  1:155198/1‑93 (MQ=255)
                                                 cTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACa                                      >  2:3901829/1‑93 (MQ=255)
                                                  tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGCACCAAAACaa                                     >  1:2337432/1‑93 (MQ=255)
                                                  tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACaa                                     >  1:4230129/1‑93 (MQ=255)
                                                  tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACaa                                     <  2:1363894/93‑1 (MQ=255)
                                                   aaCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAg                                    >  2:230763/1‑93 (MQ=255)
                                                   aaCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAg                                    <  2:2916559/93‑1 (MQ=255)
                                                   aaCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAg                                    <  1:2592052/93‑1 (MQ=255)
                                                    aCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAgg                                   <  1:2062733/93‑1 (MQ=255)
                                                     cGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGc                                  <  2:3496628/93‑1 (MQ=255)
                                                      gCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCa                                 >  1:1037358/1‑93 (MQ=255)
                                                      gCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCa                                 <  2:1118154/93‑1 (MQ=255)
                                                          ccccaATCACCAAGGAAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTTTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATctc                             <  2:1508833/91‑3 (MQ=255)
                                                            ccaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTg                           <  2:4681334/93‑1 (MQ=255)
                                                            ccaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTg                           <  1:3218095/93‑1 (MQ=255)
                                                             caATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTgg                          <  1:3077734/93‑1 (MQ=255)
                                                             caATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTgg                          <  1:3795311/93‑1 (MQ=255)
                                                              aaTCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGACACAACAAAACGACTCACTGgc                         >  1:2224377/1‑93 (MQ=255)
                                                              aaTCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgc                         <  2:4093537/93‑1 (MQ=255)
                                                               aTCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgcg                        >  2:3947598/1‑93 (MQ=255)
                                                                tCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgcgc                       >  2:1382189/1‑93 (MQ=255)
                                                                 cACCATGGCAGTAAATGAAATACCACCGAAGCGTTTCCCCGCGATGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgcgca                      <  2:4209441/93‑2 (MQ=255)
                                                                  aCCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTg                     >  1:1681116/1‑93 (MQ=255)
                                                                   ccATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGa                    <  1:952857/93‑1 (MQ=255)
                                                                    caaggCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                   <  2:1195360/90‑1 (MQ=255)
                                                                    cATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                   <  1:167963/93‑1 (MQ=255)
                                                                    cATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                   <  2:3765698/93‑1 (MQ=255)
                                                                    cATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                   <  1:3748360/93‑1 (MQ=255)
                                                                    cATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                   >  1:904595/1‑93 (MQ=255)
                                                                     aTGTCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAAt                  <  2:3228095/93‑1 (MQ=255)
                                                                        gCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATaaa               <  2:4460140/93‑1 (MQ=255)
                                                                        gCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATaaa               >  2:2194419/1‑93 (MQ=255)
                                                                          aGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAAt             <  2:2232559/93‑1 (MQ=255)
                                                                           gTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAAtt            >  1:2925491/1‑93 (MQ=255)
                                                                             aaaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaa          <  1:1132561/93‑1 (MQ=255)
                                                                              aaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGTGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaaa         >  2:4402550/1‑93 (MQ=255)
                                                                              aaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaaa         >  1:2121733/1‑93 (MQ=255)
                                                                              aaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaaa         >  2:4538955/1‑93 (MQ=255)
                                                                                tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATa       >  2:3927021/1‑93 (MQ=255)
                                                                                tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATa       >  1:2051993/1‑93 (MQ=255)
                                                                                tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATa       <  2:544519/93‑1 (MQ=255)
                                                                                tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATa       >  1:946768/1‑93 (MQ=255)
                                                                                tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATa       <  2:1602993/93‑1 (MQ=255)
                                                                                 ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa      <  1:3008963/93‑1 (MQ=255)
                                                                                 ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa      <  1:391330/93‑1 (MQ=255)
                                                                                     tACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTAc  <  2:2807305/93‑1 (MQ=255)
                                                                                     tACCACCGAAGCGTTTCCCCCCGGTGTTTCCCAGGATTATCGGGAAGGACCAAAACACGGCATCCATGGCGCTCAATAAAATTAAATAACTAc  >  2:1927188/1‑93 (MQ=255)
                                                                                         |                                                                                        
GCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTAC  >  NC_000913/3902633‑3902810

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