Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,705 C→T 16.7% intergenic (‑200/+28) 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,7050CT16.7% 466.6 / 73.0 191intergenic (‑200/+28)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (80/79);  new base T (13/19);  total (93/98)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 208 total aligned reads displayed.

AAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCAC  >  NC_000913/3322606‑3322804
                                                                                                   |                                                                                                   
aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                                     >  1:995320/1‑100 (MQ=255)
aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                                     >  2:1244399/1‑100 (MQ=255)
  gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                                   <  2:544910/100‑1 (MQ=255)
    gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                                 >  2:1611445/1‑100 (MQ=255)
    gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                                 >  2:717745/1‑100 (MQ=255)
    gAATGACCTTCGCAAGGGCAATTTAATCGCATTGTATCAATTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                                 <  2:1989808/100‑1 (MQ=255)
    gAAGGACCAGCGCAAGAGCAATGTCATAGGATTGTAGCAAAGAGAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                                 <  2:1486954/100‑1 (MQ=255)
          ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCtttttt                                                                                           <  2:253052/100‑2 (MQ=38)
          ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                           <  1:3426731/100‑1 (MQ=255)
          ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                           <  1:3866544/100‑1 (MQ=255)
          ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                           <  1:1992629/100‑1 (MQ=255)
           cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttn                                                                                          <  1:3209830/100‑3 (MQ=38)
           cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAg                                                                                          >  2:3362563/1‑100 (MQ=255)
            ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttct                                                                                         <  1:782722/100‑4 (MQ=38)
            ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGc                                                                                         <  1:707096/100‑1 (MQ=255)
             tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                        <  1:733673/100‑5 (MQ=37)
              cgtaaGGGCAATTTTATTGCATTGTATCAATTGCGAAAAAAAACAGCAGCACCTCTGTGTTCGCTGAGAAAAAAGGCTGCGATTGTCGTCttttttctca                                                                                       <  2:1447269/97‑6 (MQ=255)
                cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                                      <  2:1650035/99‑7 (MQ=37)
                cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                                     >  1:635895/1‑100 (MQ=255)
                cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                                     >  2:1341860/1‑100 (MQ=255)
                cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                                     <  2:59529/100‑1 (MQ=255)
                  aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTa                                                                                   <  1:1379741/100‑1 (MQ=255)
                        aTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                             >  1:1929060/1‑100 (MQ=255)
                        aTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                             >  2:3317861/1‑100 (MQ=255)
                          tttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGc                                                                           >  1:3158694/1‑100 (MQ=255)
                            tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCtt                                                                         <  2:241532/100‑1 (MQ=255)
                             aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCttt                                                                        <  1:1061356/100‑1 (MQ=255)
                              tcttccgatctCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                       <  2:2762637/90‑1 (MQ=255)
                              tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                       >  2:2292613/1‑100 (MQ=255)
                              tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGGTATCTAACGCTTTa                                                                       >  2:2044022/1‑100 (MQ=255)
                               cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaa                                                                      <  1:1623736/100‑1 (MQ=255)
                                gCATTGTATCAACTGCGGAAAAAAGCAGCAAAACCGAGGTGTTCGCTGCGAAAAAAGGCGTCAATCGCCGCCTTTTTAGCCAGTTATCTAACTCTTTa                                                                       <  2:530797/98‑1 (MQ=255)
                                gCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaaa                                                                     >  1:4089046/1‑100 (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)
                                         cAACTGCGGCAAAAAACAGCAAAACCGATGTGTTCGCTTAGAAAAAAGGCGACAATCGCCGCCTTTTTATCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                            >  2:3232151/1‑100 (MQ=255)
                                         cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTaagatcggaaga                                                            >  1:2762644/1‑90 (MQ=255)
                                         cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                            <  2:3969958/100‑1 (MQ=255)
                                          aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                           <  2:658057/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                                                       >  2:1720979/1‑100 (MQ=255)
                                              gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                       >  2:744304/1‑100 (MQ=255)
                                              gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                       >  1:3249475/1‑100 (MQ=255)
                                                 gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATc                                                     >  2:3111109/1‑99 (MQ=255)
                                                 gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                                    >  2:3024938/1‑100 (MQ=255)
                                                 gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                                    >  2:2686397/1‑100 (MQ=255)
                                                 gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                                    >  1:2648068/1‑100 (MQ=255)
                                                  aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                   <  1:143962/100‑1 (MQ=255)
                                                  aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                   <  1:4069882/100‑1 (MQ=255)
                                                   aaaaaaCCGCAAAACGGATGTGTGCGCTGAGAAAAAAGGCGGCAATCGCCGCCTTTTTAGCCAGTGATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                                  <  2:2892073/100‑1 (MQ=255)
                                                        aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                             <  1:4034086/100‑1 (MQ=255)
                                                          aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                           <  2:1558997/100‑1 (MQ=255)
                                                           gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGCAgt                                          <  1:508507/100‑1 (MQ=255)
                                                            cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                         <  1:744304/100‑1 (MQ=255)
                                                             aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtgt                                        <  1:443367/100‑1 (MQ=255)
                                                               agccGATGTGTTCGCTGAGAAAAATGGCGACAATCGCCGCCTGTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                                      <  2:3996274/98‑1 (MQ=255)
                                                               aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                                      >  1:178460/1‑100 (MQ=255)
                                                                aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                                     <  1:3078144/100‑1 (MQ=255)
                                                                aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACCGCATCGCCGATGcagtaaaag                                     <  1:3053015/100‑9 (MQ=255)
                                                                 ccGATGTGTTCGCTGAGAAAAAAGGCGCCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                                     <  1:2085835/99‑1 (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)
                                                                      gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTCTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                               <  1:2943905/100‑1 (MQ=255)
                                                                       tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                              <  2:2987487/100‑1 (MQ=255)
                                                                       tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                              <  2:3944125/100‑1 (MQ=255)
                                                                       tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                              <  2:4150400/100‑1 (MQ=255)
                                                                       tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                              <  2:691033/100‑1 (MQ=255)
                                                                          tCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCa                           <  1:4046630/100‑1 (MQ=255)
                                                                            gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                         >  1:3678830/1‑100 (MQ=255)
                                                                            gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                         <  1:2301765/100‑1 (MQ=255)
                                                                             cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATATAACGCGTTAAACGGCTTCTACTGCAACCGCGATGCGGTGCGCAAATTCAGTcgcat                        >  2:3708940/1‑96 (MQ=255)
                                                                               ggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                      >  2:2124277/5‑100 (MQ=37)
                                                                               gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                      <  2:1914451/100‑1 (MQ=255)
                                                                                gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                     >  1:930302/4‑100 (MQ=38)
                                                                                agaAAAAAGGCGACAATCGCCGTCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                      >  2:432989/1‑99 (MQ=255)
                                                                                agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                     >  2:1524262/1‑100 (MQ=255)
                                                                                 ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTCACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc                    >  2:838655/3‑100 (MQ=37)
                                                                                 ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTACCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc                    >  2:2239496/3‑100 (MQ=37)
                                                                                 caaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc                    >  1:2851612/2‑100 (MQ=38)
                                                                                  aaaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCt                   <  1:3011553/100‑1 (MQ=255)
                                                                                   aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                  >  2:1322730/1‑100 (MQ=255)
                                                                                   aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                  >  1:3859914/1‑100 (MQ=255)
                                                                                   aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                  >  1:2577689/1‑100 (MQ=255)
                                                                                     aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg                <  1:3058816/100‑1 (MQ=255)
                                                                                     aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg                <  1:1166159/100‑1 (MQ=255)
                                                                                      aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCCCCTGCGCTTCGt               >  2:1094413/1‑100 (MQ=255)
                                                                                        ggCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt             >  1:1224186/1‑100 (MQ=255)
                                                                                          cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc           >  1:916604/1‑100 (MQ=255)
                                                                                          cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc           >  1:2405202/1‑100 (MQ=255)
                                                                                           gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg          >  2:802217/1‑100 (MQ=255)
                                                                                            aCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc         <  1:2146341/100‑1 (MQ=255)
                                                                                             cAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTACACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc         >  2:3625009/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
AAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCAC  >  NC_000913/3322606‑3322804

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: