Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,698 A→G 15.2% intergenic (‑193/+35) 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,6980AG15.2% 394.2 / 58.7 172intergenic (‑193/+35)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (71/73);  new base G (10/16);  total (82/90)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 174 total aligned reads displayed.

CTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCG  >  NC_000913/3322599‑3322796
                                                                                                   |                                                                                                   
 tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATAAAAAAAAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGAc                                                                                                    <  2:2010529/100‑1 (MQ=255)
   aCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACaa                                                                                                  >  1:3270679/1‑100 (MQ=255)
     cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATc                                                                                                >  1:2044454/1‑100 (MQ=255)
       aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              <  1:935728/100‑1 (MQ=255)
         gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                            <  1:2261331/100‑1 (MQ=255)
           gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                          >  1:1704834/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                                                                                    <  1:3814657/100‑1 (MQ=38)
                aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑ttttt                                                                                    >  1:503645/1‑100 (MQ=255)
                aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑ttttt                                                                                    <  1:3792834/100‑1 (MQ=255)
                 ccTGCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑tttttt                                                                                   <  2:2359274/100‑2 (MQ=38)
                  cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttc                                                                                  <  1:724525/100‑3 (MQ=38)
                  attCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAGAACAGCAAAACCGCTGTGTTCGCTGAGAAAAACGGCGACAATCGCCGCC‑TTTTTAg                                                                                  >  2:2562780/2‑100 (MQ=255)
                   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:2414792/100‑5 (MQ=37)
                    tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctc                                                                                <  2:3034287/100‑5 (MQ=37)
                    tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGcc                                                                                >  2:1987543/1‑100 (MQ=255)
                     cGCAAGGGCAATTTTAGCGCATTGTATAAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctca                                                                               <  2:1141756/100‑6 (MQ=25)
                      gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                              <  2:3216198/100‑7 (MQ=37)
                      gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                              <  1:2145477/100‑7 (MQ=37)
                      gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCC‑ttttttctcag                                                                              <  2:3628856/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)
                           ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATc                                                                         >  2:1170245/1‑100 (MQ=255)
                           ggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATc                                                                         >  1:1332217/1‑100 (MQ=255)
                            gCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTAGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCt                                                                        <  1:2273553/100‑1 (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:2541785/1‑100 (MQ=255)
                                    aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCttt                                                                <  1:2162227/100‑1 (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)
                                              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)
                                                       ggAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAACAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATc                                             <  2:1015258/100‑1 (MQ=255)
                                                         aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                           <  1:1001256/100‑1 (MQ=255)
                                                         aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                           >  1:306240/1‑100 (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)
                                                             aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                       <  1:987548/100‑1 (MQ=255)
                                                              aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                      >  2:3140789/1‑100 (MQ=255)
                                                               aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGAATCGGCGATGc                                     >  2:1423565/1‑100 (MQ=255)
                                                                cAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                                    >  1:979025/1‑100 (MQ=255)
                                                                 aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                   >  1:3322041/1‑100 (MQ=255)
                                                                  gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                  >  2:1813235/1‑100 (MQ=255)
                                                                   cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                 <  2:2169534/100‑1 (MQ=255)
                                                                      aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGcagtaaaa                              <  1:2650053/100‑8 (MQ=255)
                                                                       aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTACACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                             <  2:327096/100‑1 (MQ=255)
                                                                       aCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                             <  1:2945803/100‑1 (MQ=255)
                                                                        ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCaa                            >  1:2370281/1‑100 (MQ=255)
                                                                         cGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACCGCATCGCCGATGcagtaaaagcc                           <  1:1144351/100‑11 (MQ=38)
                                                                          gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAt                          <  2:1332208/100‑1 (MQ=255)
                                                                          gATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAt                          <  2:1331233/100‑1 (MQ=255)
                                                                            tgtgTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCAAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                        <  2:2518295/100‑1 (MQ=255)
                                                                             gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTc                        <  2:620510/99‑1 (MQ=255)
                                                                             gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                       <  2:784515/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                      <  2:1198634/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                      <  2:2828777/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                      <  1:583083/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                      >  1:3506007/1‑100 (MQ=255)
                                                                                 tCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCa                   <  2:1467221/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)
                                                                                    ctggctaaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCt                >  1:2368138/7‑100 (MQ=37)
                                                                                     tGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTg               >  2:3850129/1‑100 (MQ=255)
                                                                                      ggctaaaaaGGCGGCGATTGTCGCCTTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTg               >  1:3544193/5‑100 (MQ=25)
                                                                                      gagaAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc              <  1:1372422/100‑1 (MQ=255)
                                                                                       gctaaaaaGGCGGCGATTGTCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg             >  2:1297318/4‑100 (MQ=38)
                                                                                       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‑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:3084662/1‑100 (MQ=38)
                                                                                            acaGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAACCGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg        >  2:1750423/3‑100 (MQ=255)
                                                                                            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:616868/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCC‑TTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCG  >  NC_000913/3322599‑3322796

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