Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,700 A→G 15.3% intergenic (‑195/+33) secG ← / ← glmM preprotein translocase membrane subunit/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,7000AG15.3% 500.0 / 72.1 208intergenic (‑195/+33)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (88/88);  new base G (13/19);  total (101/107)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

ACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGC  >  NC_000913/3322602‑3322797
                                                                                                  |                                                                                                 
 ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAACACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACaa                                                                                                  >  2:462417/1‑99 (MQ=255)
 ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAt                                                                                                 >  1:322873/1‑100 (MQ=255)
 ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAt                                                                                                 >  2:4053803/1‑100 (MQ=255)
   aaaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATTTGTTCACTGAAAAAAAAGGCGACAATcg                                                                                               <  2:629034/100‑1 (MQ=255)
    aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              >  1:2940855/1‑100 (MQ=255)
    aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              <  1:706988/100‑1 (MQ=255)
    aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              >  1:995320/1‑100 (MQ=255)
      gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                            <  2:544910/100‑1 (MQ=255)
        gAATGACCTTCGCAAGGGCAATTTAATCGCATTGTATCAATTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                          <  2:1989808/100‑1 (MQ=255)
        gAAGGACCAGCGCAAGAGCAATGTCATAGGATTGTAGCAAAGAGAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                          <  2:1486954/100‑1 (MQ=255)
           tGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttt                                                                                       <  2:3013635/100‑1 (MQ=255)
            gACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                      >  1:3053982/1‑100 (MQ=255)
             aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttttt                                                                                     <  1:3052876/100‑1 (MQ=255)
              ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCtttttt                                                                                    <  2:253052/100‑2 (MQ=38)
              ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                    <  1:3426731/100‑1 (MQ=255)
              ccTACGCAAGGGCAATTTTATCGCATAGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttt                                                                                     <  2:2871658/99‑1 (MQ=38)
               cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttn                                                                                   <  1:3209830/100‑3 (MQ=38)
               cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttc                                                                                   <  2:713101/100‑3 (MQ=38)
               cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAg                                                                                   >  2:3362563/1‑100 (MQ=255)
                ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttct                                                                                  <  1:782722/100‑4 (MQ=38)
                ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttct                                                                                  <  1:4182178/100‑4 (MQ=38)
                 tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                 <  1:733647/100‑5 (MQ=37)
                 tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                 <  1:3907446/100‑5 (MQ=37)
                 tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGcc                                                                                 >  2:2980649/1‑100 (MQ=255)
                  cgtaaGGGCAATTTTATTGCATTGTATCAATTGCGAAAAAAAACAGCAGCACCTCTGTGTTCGCTGAGAAAAAAGGCTGCGATTGTCGTCttttttctca                                                                                <  2:1447269/97‑6 (MQ=255)
                   gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                               <  2:67483/100‑7 (MQ=37)
                   gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                               <  1:2090942/100‑1 (MQ=255)
                    cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                              <  2:59529/100‑1 (MQ=255)
                    cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                              >  2:1341860/1‑100 (MQ=255)
                      aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTa                                                                            <  1:1379741/100‑1 (MQ=255)
                       gggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTAt                                                                           >  1:1958461/1‑100 (MQ=255)
                            aTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                      >  1:403981/1‑100 (MQ=255)
                             ttttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACg                                                                     <  1:3070189/100‑1 (MQ=255)
                                 aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCttt                                                                 <  1:1061356/100‑1 (MQ=255)
                                  tcttccgatctCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                <  2:2762637/90‑1 (MQ=255)
                                  tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGGTATCTAACGCTTTa                                                                >  2:2044022/1‑100 (MQ=255)
                                   cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaa                                                               <  1:1623736/100‑1 (MQ=255)
                                    gCATTGTATCAACTGCGGAAAAAAGCAGCAAAACCGAGGTGTTCGCTGCGAAAAAAGGCGTCAATCGCCGCCTTTTTAGCCAGTTATCTAACTCTTTa                                                                <  2:530797/98‑1 (MQ=255)
                                     cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                             <  1:2905425/100‑1 (MQ=255)
                                     cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTAATCTAACGCTTTAAAc                                                             >  2:2452907/1‑100 (MQ=255)
                                      aTTGTATCAACTGCGGAAAAAAACAGAAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTGAGCCAGTTATCTAACGCTTTAAACg                                                            >  1:1594967/1‑100 (MQ=255)
                                        tGTATCAATTGCGGAAAAAAACAGCATAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                          <  1:2491842/100‑1 (MQ=255)
                                            tCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtttt                                                      <  1:4060684/100‑1 (MQ=255)
                                             cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaagatcggaaga                                                     >  1:2762644/1‑90 (MQ=255)
                                             cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                     <  2:3969958/100‑1 (MQ=255)
                                              aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                    >  1:2235304/1‑100 (MQ=255)
                                               aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                   <  1:2346953/100‑1 (MQ=255)
                                                cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:3943065/100‑1 (MQ=255)
                                                cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:1543436/100‑1 (MQ=255)
                                                 tGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGTCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:1046149/99‑1 (MQ=255)
                                                  gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                >  1:3249475/1‑100 (MQ=255)
                                                  gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                >  2:1720979/1‑100 (MQ=255)
                                                     gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                             >  1:2648068/1‑100 (MQ=255)
                                                     gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                             >  2:2686397/1‑100 (MQ=255)
                                                      aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                            <  1:143962/100‑1 (MQ=255)
                                                          aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                        >  2:730186/1‑100 (MQ=255)
                                                            aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                      <  1:4086584/100‑1 (MQ=255)
                                                              aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                    <  2:1558997/100‑1 (MQ=255)
                                                              aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                    <  1:2980660/100‑1 (MQ=255)
                                                               gCAAAACCGATGTGTTCGCTGAGAAACAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCCTTAAACGGCTTTTACTGCAGCGGCCATGCGgt                                   >  2:3163191/1‑100 (MQ=255)
                                                               gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                    >  1:3819920/1‑99 (MQ=255)
                                                               gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                   >  1:470468/1‑100 (MQ=255)
                                                                cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                  <  1:744304/100‑1 (MQ=255)
                                                                cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCCGCGATGCGgtg                                  >  2:2320817/1‑100 (MQ=255)
                                                                   agccGATGTGTTCGCTGAGAAAAATGGCGACAATCGCCGCCTGTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                               <  2:3996274/98‑1 (MQ=255)
                                                                   aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                               >  1:178460/1‑100 (MQ=255)
                                                                     ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGcagtaaaagc                             <  2:2150550/100‑10 (MQ=255)
                                                                     ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGCAgtaaaag                              <  2:3898662/99‑6 (MQ=255)
                                                                         tgtgTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                         <  1:3362579/100‑1 (MQ=255)
                                                                          gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGTGGGGTGCAAATTCa                        >  1:2281554/1‑100 (MQ=255)
                                                                          gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                        <  2:4089017/100‑1 (MQ=255)
                                                                           tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                       <  2:3944125/100‑1 (MQ=255)
                                                                           tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                       <  2:4150400/100‑1 (MQ=255)
                                                                            gtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCAGTGTGCAAATTCAGt                      >  2:2263673/1‑100 (MQ=255)
                                                                                gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                  <  2:834112/100‑1 (MQ=255)
                                                                                gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                  <  1:2301765/100‑1 (MQ=255)
                                                                                gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                  >  1:3678830/1‑100 (MQ=255)
                                                                                 cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATATAACGCGTTAAACGGCTTCTACTGCAACCGCGATGCGGTGCGCAAATTCAGTcgcat                 >  2:3708940/1‑96 (MQ=255)
                                                                                  tGAGAAAAAAGGCGACAATCGCCGCCTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:756195/1‑100 (MQ=255)
                                                                                   gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               <  2:1914451/100‑1 (MQ=255)
                                                                                   gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:3107444/1‑100 (MQ=255)
                                                                                    gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  1:2069632/4‑100 (MQ=38)
                                                                                    gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  1:930302/4‑100 (MQ=38)
                                                                                    gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  1:577326/4‑100 (MQ=38)
                                                                                    agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  2:1524262/1‑100 (MQ=255)
                                                                                     gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:2107283/1‑100 (MQ=255)
                                                                                     gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:3292901/1‑100 (MQ=255)
                                                                                     ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:162695/3‑100 (MQ=38)
                                                                                     caaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  1:2851612/2‑100 (MQ=38)
                                                                                       aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt           >  2:1322730/1‑100 (MQ=255)
                                                                                         aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg         <  1:1166159/100‑1 (MQ=255)
                                                                                          aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCCCCTGCGCTTCGt        >  2:1094413/1‑100 (MQ=255)
                                                                                          aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGt        <  2:1667309/100‑1 (MQ=255)
                                                                                            ggCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt      >  1:1224186/1‑100 (MQ=255)
                                                                                              cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc    >  1:2405202/1‑100 (MQ=255)
                                                                                              cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc    >  1:916604/1‑100 (MQ=255)
                                                                                               gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg   <  1:3861168/100‑1 (MQ=255)
                                                                                                aCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc  <  1:3085995/100‑1 (MQ=255)
                                                                                                aCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc  <  1:2146341/100‑1 (MQ=255)
                                                                                                 cAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTACACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc  >  2:3625009/1‑99 (MQ=255)
                                                                                                  |                                                                                                 
ACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGC  >  NC_000913/3322602‑3322797

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