Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,700 A→G 16.2% 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,7000AG16.2% 410.5 / 67.5 179intergenic (‑195/+33)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (73/77);  new base G (13/16);  total (86/93)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 181 total aligned reads displayed.

AACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCC  >  NC_000913/3322601‑3322798
                                                                                                   |                                                                                                  
aaCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACa                                                                                                    >  2:2880887/1‑100 (MQ=255)
  cccatgCGAAGGTCATTCGCAAGGGCAATCTTATCGCATTAAATCAACTGCGAAAAACAACAGCAAAACCGAAGTGTTCGCTGCCAAAAAATGCGCCAAt                                                                                                  >  1:6290/6‑100 (MQ=255)
  ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAt                                                                                                  >  2:1510960/1‑100 (MQ=255)
   cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATc                                                                                                 >  1:3700131/1‑100 (MQ=255)
    aaaGCGAATAACCTTCGCAAGGGCAATTTTAACGCATTGTATCAACTTCGCAAAAAAACACCAAAACAGATATGTTAGCTGAGAAAAAAAGAGACact                                                                                                  >  1:7000/1‑96 (MQ=255)
      aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                              <  1:1467433/100‑1 (MQ=255)
       gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                             <  1:193798/100‑1 (MQ=255)
            tGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGACAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttt                                                                                        >  1:2061928/1‑100 (MQ=255)
            tGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttt                                                                                        <  1:439049/100‑1 (MQ=255)
             gACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                       >  1:768117/1‑100 (MQ=255)
              aacTTCTCAAGTGCAATTTTATCGAATTGTAACAAAAGCAGAAAAAAACAGCAAATCCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttt                                                                                      <  1:3926266/98‑1 (MQ=19)
              aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttt                                                                                      <  1:1181112/100‑1 (MQ=38)
               ccTTCGCAAGGGCAATTTTATCGCATTTTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                     <  2:534092/100‑1 (MQ=255)
               ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                     >  1:1597692/1‑100 (MQ=255)
                 ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGc                                                                                   <  2:47162/100‑1 (MQ=255)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                  <  1:3993386/100‑5 (MQ=38)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                  <  2:1740063/100‑5 (MQ=37)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGcc                                                                                  <  1:2341314/100‑1 (MQ=255)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGcc                                                                                  <  2:3177519/100‑1 (MQ=255)
                   cGCAAGGGCAATTGTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAGAGGCGGCGATTGTCGCCttttttctca                                                                                 <  2:2190650/100‑6 (MQ=37)
                    gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                                <  2:838357/100‑7 (MQ=37)
                    gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                                <  2:924849/100‑7 (MQ=37)
                    gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                <  2:182220/100‑1 (MQ=255)
                     cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCTATGTGTTCGCTGAGAAAAAAGTCGGCGATTGTCGCCttttttctcagc                                                                               <  1:560051/100‑8 (MQ=37)
                     cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                               <  1:2186143/100‑1 (MQ=255)
                      aaggcaaaTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                <  1:1602032/92‑1 (MQ=255)
                       ngggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTa                                                                             >  2:1866781/2‑100 (MQ=255)
                       aGGACAATTTTATCGCATTGTATCAACTGCGAAAAAAAACAGCAAAACCGATGTGTTCGCTTAGAAAAAAGGCGACAATCGCCGCTTTTTTAGCCAGtt                                                                              <  1:2719558/99‑1 (MQ=255)
                        gggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTAt                                                                            >  1:3533444/1‑100 (MQ=255)
                        gggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTAt                                                                            >  1:3352196/1‑100 (MQ=255)
                         ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATc                                                                           >  1:2476297/1‑100 (MQ=255)
                         ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATc                                                                           >  2:4303328/1‑100 (MQ=255)
                          gCAATTTTATCGCATTGTATCAACTGCGGAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                                         <  2:3495742/100‑1 (MQ=255)
                           cAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGGAAAAAGCCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                                         >  1:3099451/1‑100 (MQ=255)
                            aaTTTTATCGCATTGTATGAACTGAGGAAAAAAACAGCAAAACCGATTTTTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                        <  1:1284477/100‑1 (MQ=255)
                              ttttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACg                                                                      <  2:4049186/100‑1 (MQ=255)
                               tttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGc                                                                     <  2:2531753/100‑1 (MQ=255)
                               tttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGc                                                                     >  2:1173554/1‑100 (MQ=255)
                                ttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGca                                                                    >  1:2521065/1‑99 (MQ=255)
                                  aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCttt                                                                  >  2:323348/1‑100 (MQ=255)
                                   tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                 <  2:1193518/100‑1 (MQ=255)
                                   tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                 >  1:4136633/1‑100 (MQ=255)
                                    cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                 <  1:3555793/99‑1 (MQ=255)
                                     gCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTATTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATTTAACGCTTTaaa                                                               >  1:2437037/1‑100 (MQ=255)
                                          gTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCt                                                          >  2:1429563/1‑100 (MQ=255)
                                           tATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtt                                                         <  2:2590480/100‑1 (MQ=255)
                                            aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCttt                                                        >  1:3024971/1‑100 (MQ=255)
                                            aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCttt                                                        >  2:1214287/1‑100 (MQ=255)
                                              cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCtttt                                                       <  1:2230685/99‑1 (MQ=255)
                                               aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                     >  1:1014423/1‑100 (MQ=255)
                                                aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                    <  1:2145159/100‑1 (MQ=255)
                                                aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                    <  2:2083269/100‑1 (MQ=255)
                                                  tGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGc                                                  >  1:1361471/1‑100 (MQ=255)
                                                   gCGGCAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGAAAATAGCAGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTATTACTGc                                                  >  1:657658/1‑99 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGc                                                  <  1:3190837/99‑1 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                 >  2:1592797/1‑100 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                 >  1:1838212/1‑100 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                 >  2:390243/1‑100 (MQ=255)
                                                        aaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTACATCGGc                                            <  1:612787/100‑1 (MQ=255)
                                                           aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                         >  2:1675286/1‑100 (MQ=255)
                                                                gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                    >  2:3726520/1‑100 (MQ=255)
                                                                gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                    >  1:230118/1‑100 (MQ=255)
                                                                gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAGCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGCAgt                                    <  2:4111471/100‑1 (MQ=255)
                                                                 cacaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTCGCCAGTTATCTAACCCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                      >  1:2625330/4‑97 (MQ=255)
                                                                 cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                   >  1:1925029/1‑100 (MQ=255)
                                                                  aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtgt                                  >  1:1513399/1‑100 (MQ=255)
                                                                   aaaCCGATGTGTTCGCTGAGACAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATACGgtgt                                  >  1:1396036/1‑99 (MQ=255)
                                                                   aaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGtgtg                                 >  1:1899315/1‑100 (MQ=255)
                                                                   aaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGcagtaaa                                 <  2:4020954/100‑7 (MQ=255)
                                                                    aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                                <  2:552979/100‑1 (MQ=255)
                                                                      ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaa                              >  1:3842735/1‑100 (MQ=255)
                                                                       cGATGTGTTCGCTGAGAAAAAATGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaaa                             <  2:4376670/100‑1 (MQ=255)
                                                                       cGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaaa                             <  2:3892626/100‑1 (MQ=255)
                                                                        gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAt                            <  1:1377365/100‑1 (MQ=255)
                                                                            tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCac                        <  2:674776/100‑2 (MQ=255)
                                                                              ttCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTc                      <  2:247712/100‑1 (MQ=255)
                                                                                cGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCACATTCAGTCAc                    >  1:2736173/1‑100 (MQ=255)
                                                                                cGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAc                    >  1:2237857/1‑100 (MQ=255)
                                                                                 actggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                   >  1:257838/8‑100 (MQ=25)
                                                                                 actggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                   >  2:1358933/8‑100 (MQ=37)
                                                                                  ctggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCt                  >  2:2601297/7‑100 (MQ=37)
                                                                                  cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCt                  >  2:2044946/1‑100 (MQ=255)
                                                                                    ggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                >  1:2679797/5‑100 (MQ=37)
                                                                                    gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                >  1:920695/1‑100 (MQ=255)
                                                                                    gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAACTg                 >  1:1756773/1‑99 (MQ=255)
                                                                                     gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg               >  2:6497/4‑100 (MQ=38)
                                                                                     gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg               >  2:3763798/4‑100 (MQ=38)
                                                                                     gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg               >  2:1136439/4‑100 (MQ=38)
                                                                                     gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                >  1:4152293/4‑99 (MQ=38)
                                                                                      gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc              >  1:2950645/1‑100 (MQ=255)
                                                                                      gaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc              <  2:1638289/100‑1 (MQ=255)
                                                                                       aaaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCt             <  2:150505/100‑1 (MQ=255)
                                                                                        aaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt            >  2:1569576/1‑100 (MQ=38)
                                                                                         aaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTc           >  1:3170674/1‑100 (MQ=255)
                                                                                          aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGGTTc           >  1:1695181/1‑99 (MQ=255)
                                                                                            aGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTc        >  1:3844646/1‑100 (MQ=255)
                                                                                              gCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCtt      <  1:1128187/100‑1 (MQ=255)
                                                                                               cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc     >  1:424477/1‑100 (MQ=255)
                                                                                                gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg    <  2:1958935/100‑1 (MQ=255)
                                                                                                 aCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc   <  1:1670710/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCC  >  NC_000913/3322601‑3322798

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: