Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,703 C→T 15.6% intergenic (‑198/+30) 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,7030CT15.6% 412.2 / 56.0 166intergenic (‑198/+30)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (67/73);  new base T (10/16);  total (77/89)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 175 total aligned reads displayed.

CAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCC  >  NC_000913/3322604‑3322802
                                                                                                   |                                                                                                    
cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATc                                                                                                      >  1:2044454/1‑100 (MQ=255)
 aaaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATAAACTGCGGAAAAAAAAAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcg                                                                                                     >  2:1300652/1‑100 (MQ=255)
  aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                                    <  1:935728/100‑1 (MQ=255)
   aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                                   <  2:1910677/100‑1 (MQ=255)
      gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAGGGCGACAATCgccgcc                                                                                                >  1:347311/1‑100 (MQ=255)
      gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                                <  1:1975905/100‑1 (MQ=255)
         tGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACTATCGCCGCC‑ttt                                                                                            <  1:2593868/100‑1 (MQ=255)
          gACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑tttt                                                                                           >  1:2843148/1‑100 (MQ=255)
           aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttt                                                                                          <  2:1874222/100‑1 (MQ=38)
           aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttt                                                                                          <  1:3814657/100‑1 (MQ=38)
            ccTTCGAAAGACCACTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCACAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttt                                                                                          <  1:3273444/99‑1 (MQ=37)
            ccTGCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑tttttt                                                                                         <  2:2359274/100‑2 (MQ=38)
             cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttc                                                                                        <  2:2340518/100‑3 (MQ=38)
             cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttc                                                                                        <  1:724525/100‑3 (MQ=38)
              ttCGCAAGGGCAATTTTATCGCATTGTATTAACTGCGGAAAAAAACAGCCAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttct                                                                                       <  1:717190/100‑4 (MQ=37)
              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttct                                                                                       <  2:2409874/100‑4 (MQ=38)
              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGc                                                                                       >  1:879816/1‑100 (MQ=255)
              ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGc                                                                                       >  2:3720730/1‑100 (MQ=255)
              ttCGCAAAGACAATTTTATCGCATTGTATCAACTGCAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGTGATTGTCGCC‑ttttttct                                                                                       <  2:806366/100‑4 (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)
                cGCAAGGGCAATTTTAGCGCATTGTATAAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctca                                                                                     <  2:1141756/100‑6 (MQ=25)
                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                                    <  2:3628856/100‑7 (MQ=37)
                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                                    <  2:3216198/100‑7 (MQ=37)
                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                                    <  1:2145477/100‑7 (MQ=37)
                 gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAg                                                                                    >  1:1197524/1‑100 (MQ=255)
                  cAAGAGCAACTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGt                                                                                   <  2:3071804/100‑1 (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)
                               aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGGTGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                      <  1:1640355/100‑1 (MQ=255)
                               aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAACAGGCGACACTCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                      <  1:3144850/100‑1 (MQ=255)
                               aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                      >  1:2233005/1‑100 (MQ=255)
                               aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                      <  1:2162227/100‑1 (MQ=255)
                                tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTa                                                                     >  2:1576308/1‑100 (MQ=255)
                                 cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                    <  2:3244136/100‑1 (MQ=255)
                                 cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                    >  2:2751786/1‑100 (MQ=255)
                                 cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                    <  1:3363056/100‑1 (MQ=255)
                                 cGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAGCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTa                                                                     <  1:1930941/99‑1 (MQ=255)
                                 cGCATTGTATCAACTGAGGAAAAAAACAGCAAATCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTaa                                                                    <  1:3941910/100‑1 (MQ=255)
                                   cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                  <  1:1915806/100‑1 (MQ=255)
                                     ttGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACgg                                                                <  2:2174415/100‑1 (MQ=255)
                                       gTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCt                                                              >  1:418609/1‑100 (MQ=255)
                                         aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCttt                                                            >  1:2338033/1‑100 (MQ=255)
                                            aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                         >  1:1698274/1‑100 (MQ=255)
                                            aaCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                         <  2:672689/100‑1 (MQ=255)
                                                gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                     >  1:3454465/1‑100 (MQ=255)
                                                    aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                 <  1:1001256/100‑1 (MQ=255)
                                                         aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                            >  2:3140789/1‑100 (MQ=255)
                                                         aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                            <  2:1197515/100‑1 (MQ=255)
                                                          aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGCCAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                           <  1:2768313/100‑1 (MQ=255)
                                                          aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                           <  1:3789589/100‑1 (MQ=255)
                                                           cAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                          >  1:979025/1‑100 (MQ=255)
                                                             gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAGCTCTTTAAACGGCTTTTACTGCATCGCCGATGCAgt                                        <  2:1443275/100‑1 (MQ=255)
                                                              cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                       <  2:1100881/100‑1 (MQ=255)
                                                                 aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtgt                                      <  2:214553/98‑1 (MQ=255)
                                                                 aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACCGCATCGCCGATGcagtaaaa                                    <  1:1957262/100‑8 (MQ=255)
                                                                  aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTACACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                                   <  2:327096/100‑1 (MQ=255)
                                                                  aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                                   <  2:2541768/100‑1 (MQ=255)
                                                                  aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                                   <  1:2945803/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:583083/100‑1 (MQ=255)
                                                                         tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                            <  2:2828777/100‑1 (MQ=255)
                                                                             cGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAc                        <  1:1237408/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)
                                                                                 gagaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc                    <  1:1372422/100‑1 (MQ=255)
                                                                                  gctaaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCCATGCGGTGTGCAAATTCAGTCACCTgc                    >  2:27612/4‑99 (MQ=37)
                                                                                  agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCACATTCAGTCACCTgcg                   >  1:3891568/1‑100 (MQ=255)
                                                                                  agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                   >  1:1276833/1‑100 (MQ=255)
                                                                                  agaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGGTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                   >  1:3325252/1‑100 (MQ=255)
                                                                                    aaaaaaGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCt                 >  1:3788822/1‑100 (MQ=255)
                                                                                     aaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                >  1:3084662/1‑100 (MQ=38)
                                                                                     aaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                >  1:996756/1‑100 (MQ=38)
                                                                                       aaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg              >  1:1000643/1‑100 (MQ=255)
                                                                                          ggCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt           >  1:3652066/1‑100 (MQ=255)
                                                                                           gCGACAATCGCCGGC‑TTTTTAGCCAGTTATCTAACGCTGTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCtt          >  2:2437134/1‑100 (MQ=255)
                                                                                            cGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGTATCGGCGATGCGGTGTGCAGATTCAGTAACCTGCGCTTCGtgtg          >  2:3577300/1‑96 (MQ=255)
                                                                                             gACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg        <  1:616868/100‑1 (MQ=255)
                                                                                             gACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg        >  1:3039946/1‑100 (MQ=255)
                                                                                               cAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc      >  2:1628179/1‑100 (MQ=255)
                                                                                               cAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc      >  1:986007/1‑100 (MQ=255)
                                                                                               cAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc      <  1:3286820/100‑1 (MQ=255)
                                                                                                 aTCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCt     >  2:1802704/1‑99 (MQ=255)
                                                                                                  tCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTc   <  2:260510/100‑1 (MQ=255)
                                                                                                  tCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTc   >  1:2265287/1‑100 (MQ=255)
                                                                                                   cgccgcC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcc  <  2:2996781/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
CAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCC  >  NC_000913/3322604‑3322802

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: ‑