Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,101,728 (G)7→8 coding (39/984 nt) kdgT → 2‑keto‑3‑deoxy‑D‑gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,7211.G95.7% 335.2 / 8.5 94coding (32/984 nt)kdgT2‑keto‑3‑deoxy‑D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (54/36);  total (55/39)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 114 total aligned reads displayed.

CATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGG  >  NC_000913/4101622‑4101814
                                                                                                    |                                                                                               
cATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATccc                                                                                                  >  1:1021222/1‑100 (MQ=255)
 aTCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑g                                                                                               >  2:852784/1‑100 (MQ=255)
   cAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                             <  2:1653571/100‑1 (MQ=255)
   cAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                             >  2:198320/1‑100 (MQ=255)
   cAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggg                                                                                             <  1:1481817/100‑1 (MQ=255)
    aaCACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggg                                                                                            >  1:2837364/1‑100 (MQ=255)
    aaCACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑gggg                                                                                            <  1:2153661/100‑1 (MQ=255)
     acacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg                                                                                           <  1:2247367/100‑1 (MQ=255)
     acacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg                                                                                           <  2:2518555/100‑1 (MQ=255)
     acacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg                                                                                           >  1:1934583/1‑100 (MQ=255)
      cacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑g                                                                                         <  1:2054068/100‑1 (MQ=255)
      cacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑g                                                                                         <  1:1511622/100‑1 (MQ=255)
      cacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑g                                                                                         <  2:2619854/100‑1 (MQ=255)
        cTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑ggg                                                                                       >  2:691407/1‑99 (MQ=255)
         tGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAATATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑ggga                                                                                      <  2:771933/100‑3 (MQ=255)
           aTGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑gggatg                                                                                    <  2:1916040/100‑5 (MQ=255)
            tGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑gggatga                                                                                   <  1:497023/100‑6 (MQ=255)
            tGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑gggatga                                                                                   <  1:2347166/100‑6 (MQ=255)
            tGCTCCATTGAGGAATTACGCATCAGACCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg‑gggatga                                                                                   <  2:1421626/100‑6 (MQ=255)
             gCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑Ggatgat                                                                                   >  2:2952951/1‑100 (MQ=255)
              cTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGatgatg                                                                                  >  1:2543617/1‑100 (MQ=255)
                ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCt                                                                                <  2:2908815/100‑1 (MQ=255)
                ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCt                                                                                >  2:39290/1‑100 (MQ=255)
                 cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTc                                                                               <  2:2501579/100‑1 (MQ=255)
                     gAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTcc                                                                           >  2:2909616/1‑100 (MQ=255)
                     gAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTcc                                                                           >  1:425850/1‑100 (MQ=255)
                     gAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTcc                                                                           >  1:2930865/1‑100 (MQ=255)
                       ggAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCg                                                                         >  2:382071/1‑100 (MQ=255)
                          aTTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTa                                                                      >  1:862409/1‑100 (MQ=255)
                             aCGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTc                                                                   >  2:1246252/1‑100 (MQ=255)
                             aCGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTc                                                                   <  1:1011531/100‑1 (MQ=255)
                              cGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTcc                                                                  >  1:1660628/1‑100 (MQ=255)
                              cGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTcc                                                                  >  2:2644382/1‑100 (MQ=255)
                                  tCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTgg                                                              <  1:2205513/100‑1 (MQ=255)
                                   cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGgc                                                             <  1:1011878/100‑1 (MQ=255)
                                    aGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGgcg                                                            <  1:2684600/100‑1 (MQ=255)
                                     gCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTAATCCTTGgcgc                                                           >  1:697805/1‑100 (MQ=255)
                                       ccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGCTGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                         >  2:1068668/1‑100 (MQ=255)
                                        cTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACt                                                        >  1:1521632/1‑100 (MQ=255)
                                           ccgatcttGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgtg                                                     <  1:1595402/93‑1 (MQ=255)
                                           aaaaaTATGGCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgtg                                                     <  2:697802/100‑1 (MQ=255)
                                           aaaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgtg                                                     >  2:2341297/1‑100 (MQ=255)
                                           aaaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACtgtg                                                     <  1:198322/100‑1 (MQ=255)
                                            aaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGc                                                    >  2:459319/1‑100 (MQ=255)
                                            aaaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGc                                                    >  1:988176/1‑100 (MQ=255)
                                             aaaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGcc                                                   <  1:809680/100‑1 (MQ=255)
                                                   gCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                             <  1:1027902/100‑1 (MQ=255)
                                                   gCCGACAGGGGATGGAAATGCAGGTAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGcc                                                   <  2:2183861/94‑1 (MQ=255)
                                                      gACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                          <  1:2687528/100‑1 (MQ=255)
                                                      gACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtct                                          >  2:2738872/1‑100 (MQ=255)
                                                       aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtctc                                         >  2:2193227/1‑100 (MQ=255)
                                                       aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCg                                        <  2:1528960/100‑1 (MQ=255)
                                                        cAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCg                                        >  1:151992/1‑100 (MQ=255)
                                                          ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                      >  2:797061/1‑100 (MQ=255)
                                                          ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                      >  1:1945038/1‑100 (MQ=255)
                                                          ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                      >  1:1403469/1‑100 (MQ=255)
                                                          ggTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                      <  1:2568131/100‑1 (MQ=255)
                                                           gtggtgGAGTTGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTAGTACCGCTATTCCTTGGCGCACTGTGAGACACCTTCTCGc                                       >  2:1299393/1‑98 (MQ=255)
                                                           gTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGcc                                      <  2:2161212/100‑1 (MQ=255)
                                                           gTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCTCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGc                                       >  2:147176/1‑98 (MQ=255)
                                                           gTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCg                                     >  2:2241618/1‑100 (MQ=255)
                                                             naTGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTACCGCCATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                                   >  1:524011/2‑100 (MQ=255)
                                                             gATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggg                                   >  2:1622966/1‑100 (MQ=255)
                                                              aTGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCgggg                                  >  1:1339729/1‑100 (MQ=255)
                                                               tGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCggggg                                 >  2:186839/1‑100 (MQ=255)
                                                                 gAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCg                               >  2:1583868/1‑100 (MQ=255)
                                                                  aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgg                              >  1:329983/1‑100 (MQ=255)
                                                                  aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgg                              <  1:852787/100‑1 (MQ=255)
                                                                   aaTGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGAGGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCggg                             <  1:860492/100‑1 (MQ=255)
                                                                    aTGCAAATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCgggg                            <  2:2930849/100‑1 (MQ=255)
                                                                     tGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGCTAATGGTAGTCCCGCTATTTCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGa                           >  1:642960/1‑100 (MQ=255)
                                                                      gCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaa                          <  1:1207847/100‑1 (MQ=255)
                                                                       cAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaaa                         <  1:39291/100‑1 (MQ=255)
                                                                         gATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                       >  2:1302016/1‑100 (MQ=255)
                                                                         gATAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAata                       >  2:2209935/1‑100 (MQ=255)
                                                                           tAAAACGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtt                     >  1:1921402/1‑100 (MQ=255)
                                                                            aaaaCGCTCGATTGAGAAAAGCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAttt                    <  2:772664/100‑1 (MQ=255)
                                                                               aCGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg                 <  2:1834538/100‑1 (MQ=255)
                                                                               aCGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg                 <  1:862940/100‑1 (MQ=255)
                                                                               aCGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTgg                 <  1:120262/100‑1 (MQ=255)
                                                                                cGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa                <  1:1711747/100‑1 (MQ=255)
                                                                                cGCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa                <  1:2373992/100‑1 (MQ=255)
                                                                                cGCTCGATTGAGAAAATCCC‑‑GGGGGTGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCAAACCTTCTCGCCGGGGGCGGGGAAATATTTTg                  <  2:731586/98‑1 (MQ=255)
                                                                                 gCTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGAt               >  1:1295299/1‑100 (MQ=255)
                                                                                  cTCGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATc              <  1:397956/100‑1 (MQ=255)
                                                                                    cGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCan            <  1:380800/100‑2 (MQ=255)
                                                                                    cGATTGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCAt            >  1:200058/1‑100 (MQ=255)
                                                                                      aTTGAGAAAATCCC‑‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCa         >  1:422668/1‑100 (MQ=255)
                                                                                      aTTGAGAAAAGCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATGCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTc          <  1:685989/100‑1 (MQ=255)
                                                                                        tGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAc        >  1:2335105/1‑100 (MQ=255)
                                                                                        tGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAc        <  1:2811261/100‑1 (MQ=255)
                                                                                        tGAGAAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAc        >  2:2061767/1‑100 (MQ=255)
                                                                                         gagaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAcc       <  2:1682885/100‑1 (MQ=255)
                                                                                         gagaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAcc       >  1:36053/1‑100 (MQ=255)
                                                                                         gagaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAcc       <  2:1021217/100‑1 (MQ=255)
                                                                                          agaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACGCCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCa      >  2:2091358/1‑100 (MQ=255)
                                                                                          agaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCGTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCa      >  2:2565623/1‑100 (MQ=255)
                                                                                           gaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCAcc       <  2:1315850/98‑1 (MQ=255)
                                                                                           gaAAATCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCaa     >  1:1524870/1‑100 (MQ=255)
                                                                                              aaTCCCG‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACgg  <  1:2190970/100‑1 (MQ=255)
                                                                                                    |                                                                                               
CATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGG‑GGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGG  >  NC_000913/4101622‑4101814

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