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% 403.7 / 60.5 171intergenic (‑195/+33)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (70/74);  new base G (10/16);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 174 total aligned reads displayed.

AACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCC  >  NC_000913/3322601‑3322798
                                                                                                   |                                                                                                   
aaCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACa                                                                                                     >  1:2519233/1‑100 (MQ=255)
 aCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACaa                                                                                                    >  1:1198643/1‑100 (MQ=255)
   cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATc                                                                                                  >  1:2044454/1‑100 (MQ=255)
    aaaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATAAACTGCGGAAAAAAAAAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcg                                                                                                 >  2:1300652/1‑100 (MQ=255)
     aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                                <  1:935728/100‑1 (MQ=255)
        cGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccgc                                                                                             >  1:1276637/1‑100 (MQ=255)
         gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                            <  1:1975905/100‑1 (MQ=255)
              aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttt                                                                                      <  1:3814657/100‑1 (MQ=38)
              aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑ttttt                                                                                      >  1:503645/1‑100 (MQ=255)
              aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑ttttt                                                                                      <  1:3792834/100‑1 (MQ=255)
               ccTTCGAAAGACCACTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCACAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttt                                                                                      <  1:3273444/99‑1 (MQ=37)
               ccTGCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑tttttt                                                                                     <  2:2359274/100‑2 (MQ=38)
                cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttc                                                                                    <  1:724525/100‑3 (MQ=38)
                cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAg                                                                                    >  1:3259478/1‑100 (MQ=255)
                attCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAGAACAGCAAAACCGCTGTGTTCGCTGAGAAAAACGGCGACAATCGCCGCC‑TTTTTAg                                                                                    >  2:2562780/2‑100 (MQ=255)
                 ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGc                                                                                   >  2:3720730/1‑100 (MQ=255)
                 ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGc                                                                                   >  1:879816/1‑100 (MQ=255)
                 ttCGCAAAGACAATTTTATCGCATTGTATCAACTGCAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGTGATTGTCGCC‑ttttttct                                                                                   <  2:806366/100‑4 (MQ=37)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctc                                                                                  <  2:2414792/100‑5 (MQ=37)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctc                                                                                  <  2:3034287/100‑5 (MQ=37)
                  tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGcc                                                                                  >  2:1987543/1‑100 (MQ=255)
                   cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGACAAAAGGCGACAATCGCCGCC‑TTTTTAGcc                                                                                  >  2:1267056/1‑99 (MQ=255)
                   cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCa                                                                                 >  1:1467230/1‑100 (MQ=255)
                   cGCAAGGGCAATTTTAGCGCATTGTATAAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctca                                                                                 <  2:1141756/100‑6 (MQ=25)
                    gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                                <  2:3216198/100‑7 (MQ=37)
                    gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                                <  2:3628856/100‑7 (MQ=37)
                     cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGt                                                                               <  1:3128936/100‑1 (MQ=255)
                     cAAGAGCAACTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGt                                                                               <  2:3071804/100‑1 (MQ=255)
                      aaGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGtt                                                                              >  1:2928251/1‑100 (MQ=255)
                         ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATc                                                                           >  1:1331242/1‑100 (MQ=255)
                          gCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCt                                                                          <  1:2273243/100‑1 (MQ=255)
                            aaTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTaa                                                                        <  1:1146188/100‑1 (MQ=255)
                            aaTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTaa                                                                        >  1:1119582/1‑100 (MQ=255)
                               tttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGc                                                                     <  1:695585/100‑1 (MQ=255)
                                 tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCtt                                                                   >  2:2917901/1‑100 (MQ=255)
                                 tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCtt                                                                   <  1:1645852/100‑1 (MQ=255)
                                  aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                  <  1:2162045/100‑1 (MQ=255)
                                  aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                  >  1:2233005/1‑100 (MQ=255)
                                   tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTa                                                                 >  2:1576308/1‑100 (MQ=255)
                                    cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                >  2:2751786/1‑100 (MQ=255)
                                    cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                <  1:3363056/100‑1 (MQ=255)
                                    cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                <  2:3244136/100‑1 (MQ=255)
                                    cGCATTGTATCAACTGAGGAAAAAAACAGCAAATCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                <  1:3941910/100‑1 (MQ=255)
                                      cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAAc                                                              <  1:1915806/100‑1 (MQ=255)
                                        ttGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACgg                                                            <  1:2180176/100‑1 (MQ=255)
                                            aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCttt                                                        >  1:1188127/1‑100 (MQ=255)
                                            aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCttt                                                        >  1:2338033/1‑100 (MQ=255)
                                            aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTAGCTGAGAAAAAAAGCGCAAAACGCCGCC‑CTTTAAGCGAGTAATCTAAAGCTATAAACgg                                                            >  2:2490420/1‑96 (MQ=255)
                                              cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                      <  1:1974656/100‑1 (MQ=255)
                                               aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                     >  1:1698274/1‑100 (MQ=255)
                                               aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                     <  2:672689/100‑1 (MQ=255)
                                                aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                    >  2:2970514/1‑100 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                 <  1:717315/100‑1 (MQ=255)
                                                   gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                 >  1:3454465/1‑100 (MQ=255)
                                                       aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                             <  1:1001256/100‑1 (MQ=255)
                                                        aaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                            <  1:1953598/100‑1 (MQ=255)
                                                        aaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                            <  2:1900011/100‑1 (MQ=255)
                                                          aaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGa                                          <  2:237957/100‑1 (MQ=255)
                                                          aaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGa                                          >  2:2213842/1‑100 (MQ=255)
                                                            aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                        <  2:1197515/100‑1 (MQ=255)
                                                            aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                        >  2:3140789/1‑100 (MQ=255)
                                                             aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGCCAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                       <  1:2768313/100‑1 (MQ=255)
                                                              cAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                      >  1:979025/1‑100 (MQ=255)
                                                               aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                     >  1:3322041/1‑100 (MQ=255)
                                                                gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAGCTCTTTAAACGGCTTTTACTGCATCGCCGATGCAgt                                    <  2:1443275/100‑1 (MQ=255)
                                                                gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                    >  2:1813235/1‑100 (MQ=255)
                                                                 cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                   <  2:2169534/100‑1 (MQ=255)
                                                                 cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                   <  2:1100881/100‑1 (MQ=255)
                                                                    aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGcagtaaaa                                <  1:2650053/100‑8 (MQ=255)
                                                                     aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTACACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                               <  2:327096/100‑1 (MQ=255)
                                                                      ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaa                              >  1:2370281/1‑100 (MQ=255)
                                                                        gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAt                            <  2:1331233/100‑1 (MQ=255)
                                                                          tgtgTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCAAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                          <  2:2518295/100‑1 (MQ=255)
                                                                          tgtgTTCCCTGAGAAAAAAGGCGACAATCGCTGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                          <  2:1363829/100‑1 (MQ=255)
                                                                           gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                          <  2:620510/99‑1 (MQ=255)
                                                                           gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                         <  2:784515/100‑1 (MQ=255)
                                                                            tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                        >  1:3506007/1‑100 (MQ=255)
                                                                            tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                        <  2:2828777/100‑1 (MQ=255)
                                                                            tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                        <  1:583083/100‑1 (MQ=255)
                                                                             gtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGt                       <  1:130470/100‑1 (MQ=255)
                                                                                cGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTCATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAc                    <  2:3009011/100‑1 (MQ=255)
                                                                                 gCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                   >  2:2682427/1‑100 (MQ=255)
                                                                                  cTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCt                  >  1:2577284/1‑100 (MQ=255)
                                                                                   tGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTg                 >  2:3850129/1‑100 (MQ=255)
                                                                                    ggctaaaaaGGCGGCGATTGTCGCCTTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTg                 >  1:3544193/5‑100 (MQ=25)
                                                                                    ggctaaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                >  1:134089/5‑100 (MQ=37)
                                                                                     gctaaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCCATGCGGTGTGCAAATTCAGTCACCTgc                >  2:27612/4‑99 (MQ=37)
                                                                                     agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCACATTCAGTCACCTgcg               >  1:3891568/1‑100 (MQ=255)
                                                                                     agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                >  2:2048653/1‑99 (MQ=255)
                                                                                     agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGGTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg               >  1:3325252/1‑100 (MQ=255)
                                                                                       aaaaaaGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCt             >  1:3788822/1‑100 (MQ=255)
                                                                                        aaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt            >  1:996756/1‑100 (MQ=38)
                                                                                        aaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt            >  1:3084662/1‑100 (MQ=38)
                                                                                          aaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg          >  1:1000643/1‑100 (MQ=255)
                                                                                             ggCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt       >  1:3652066/1‑100 (MQ=255)
                                                                                               cGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGTATCGGCGATGCGGTGTGCAGATTCAGTAACCTGCGCTTCGtgtg      >  2:3577300/1‑96 (MQ=255)
                                                                                                gACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg    >  1:3039946/1‑100 (MQ=255)
                                                                                                gACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg    <  1:616868/100‑1 (MQ=255)
                                                                                                  cAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc  >  2:1628179/1‑100 (MQ=255)
                                                                                                  cAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc  >  1:986007/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCC  >  NC_000913/3322601‑3322798

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