Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,927 G→T 6.5% A879D (GCC→GAC)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,147,9270GT 6.5% 294.1 / 5.8 93A879D (GCC→GAC) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base T (2/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CTTCTTGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTTGCTG  >  NC_002947/2147778‑2148076
                                                                                                                                                     |                                                                                                                                                     
cttcttGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTgg                                                                                                                                                       <  2:199219/150‑1 (MQ=255)
 ttcttGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGc                                                                                                                                                      >  1:882924/1‑150 (MQ=255)
 ttcttGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGc                                                                                                                                                      >  2:2007407/1‑150 (MQ=255)
 ttcttGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGc                                                                                                                                                      <  2:2147009/150‑1 (MQ=255)
 ttcttGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGc                                                                                                                                                      >  2:2166038/1‑150 (MQ=255)
     tGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAca                                                                                                                                                  <  1:968261/150‑1 (MQ=255)
      gAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAcag                                                                                                                                                 >  2:1433801/1‑150 (MQ=255)
       aGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGTCGCcagc                                                                                                                                                <  1:2239373/150‑1 (MQ=255)
          tttGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCGTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGTCGCCagcagc                                                                                                                                             <  2:193926/150‑1 (MQ=255)
             ggtggCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGTCGCCagcagcaac                                                                                                                                          <  1:2139281/150‑4 (MQ=255)
                 gCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGcc                                                                                                                                            >  1:1489446/1‑144 (MQ=255)
                 gCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGcc                                                                                                                                            <  2:1489446/144‑1 (MQ=255)
                  ccGCTTCGACCTGCTCGACCGGTTGGGCGGCTTGGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACgg                                                                                                                                     >  2:414736/1‑150 (MQ=255)
                   cGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGc                                                                                                                                    >  1:888377/1‑150 (MQ=255)
                   cGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGc                                                                                                                                    >  2:1071913/1‑150 (MQ=255)
                    gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAca                                                                                                                                                  <  1:478259/135‑1 (MQ=255)
                    gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAca                                                                                                                                                  >  2:478259/1‑135 (MQ=255)
                    gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCa                                                                                                                                   >  2:1842241/1‑150 (MQ=255)
                    gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCa                                                                                                                                   >  2:1511690/1‑150 (MQ=255)
                     cTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCaa                                                                                                                                  <  1:1842241/150‑1 (MQ=255)
                      ttCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAg                                                                                                                                 <  1:888007/150‑1 (MQ=255)
                      ttCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAg                                                                                                                                 >  2:145416/1‑150 (MQ=255)
                         gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGcc                                                                                                                              <  2:927780/150‑1 (MQ=255)
                         gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGc                                                                                                                               >  1:1920008/1‑149 (MQ=255)
                         gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGc                                                                                                                               <  2:1920008/149‑1 (MQ=255)
                            cTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGAACTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCcagc                                                                                                                           >  2:881186/1‑150 (MQ=255)
                             tGCTCTACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTACCACGGTCGGCGTGCTGGTTGGCCCGGGCTTCTGCCTCGGAGCTGACGGCGCTGCTCGCGACAGCAGCCGTCACGGCAAGGCcagca                                                                                                                          <  2:133785/150‑1 (MQ=255)
                                 cGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCa                                                                                                                      <  1:74563/150‑1 (MQ=255)
                                     cGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTc                                                                                                                  >  2:470492/1‑150 (MQ=255)
                                       gTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCa                                                                                                                                        >  2:1247058/1‑126 (MQ=255)
                                       gTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCa                                                                                                                                        <  1:1247058/126‑1 (MQ=255)
                                       gTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCgg                                                                                                                >  1:1970095/1‑150 (MQ=255)
                                            gcggcTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCAGGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGGGCTGATGTTGCTGCTGGCGACAGCAGACGTCACGGCAAGGCCAGCAGCCAGTTCGGCAccc                                                                                                           >  1:1083616/1‑150 (MQ=255)
                                                 ttCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCTGCGGTAGCACGCGCGGCCTGCTGGTTGGCCTGGGCTTCTTCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTAGGCACCCAGCTc                                                                                                      >  1:906641/1‑150 (MQ=255)
                                                    gCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGt                                                                                                   <  1:1071913/150‑1 (MQ=255)
                                                    gCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGt                                                                                                   <  1:1549884/150‑1 (MQ=255)
                                                    gCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTAGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGt                                                                                                   >  2:1736585/1‑150 (MQ=255)
                                                     cTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTg                                                                                                  <  2:1970095/150‑1 (MQ=255)
                                                     cTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTg                                                                                                  <  1:1625185/150‑1 (MQ=255)
                                                        ttgtCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCg                                                                                                     <  2:885586/144‑1 (MQ=255)
                                                        ttgtCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCg                                                                                                     >  1:885586/1‑144 (MQ=255)
                                                          gtCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTgg                                                                                                 >  1:2271967/1‑146 (MQ=255)
                                                          gtCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTgg                                                                                                 <  2:2271967/146‑1 (MQ=255)
                                                           tCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTgg                                                                                            >  1:1786640/1‑150 (MQ=255)
                                                           tCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTgg                                                                                            >  1:1049629/1‑150 (MQ=255)
                                                           tCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTgg                                                                                            >  2:1436028/1‑150 (MQ=255)
                                                            cGGTTTCTGCAACGGCAGCGACCTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGt                                                                                           <  1:881186/150‑1 (MQ=255)
                                                              gTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGc                                                                                         >  1:87740/1‑150 (MQ=255)
                                                              gTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGc                                                                                         >  1:616626/1‑150 (MQ=255)
                                                                    gCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCtgttgt                                                                                   <  1:470492/150‑1 (MQ=255)
                                                                    gCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTAACGGCAAGGCCAGCAGCCAGTTCGGCACCCGGCTCTGTGGCCTGGTGCtgttgt                                                                                   <  2:1836168/150‑1 (MQ=255)
                                                                     cAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTtgtg                                                                                  <  2:438461/150‑1 (MQ=255)
                                                                      aaCGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTgg                                                                                 <  1:1433801/150‑1 (MQ=255)
                                                                       aCGGCAGCGACTTCTGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGCCGACAGCACCCGTCACGGCAAGGCCAGCACCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGgc                                                                                <  1:52516/150‑1 (MQ=255)
                                                                              cGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCtt                                                                         <  1:145416/150‑1 (MQ=255)
                                                                                  ttCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCt                                                                     >  2:1676821/1‑150 (MQ=255)
                                                                                     ggcggTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGcc                                                                  <  1:1511690/150‑1 (MQ=255)
                                                                                     ggcggTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGcc                                                                  <  2:1956437/150‑1 (MQ=255)
                                                                                          gcgCACGCTCGCCCTGTTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCcttctt                                                             <  2:1083616/148‑1 (MQ=255)
                                                                                                gCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGt                                                                                           <  2:2359058/114‑1 (MQ=255)
                                                                                                gCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGt                                                                                           >  1:2359058/1‑114 (MQ=255)
                                                                                                    ggCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGcc                                                   >  1:1973014/1‑150 (MQ=255)
                                                                                                     gCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCg                                                  >  2:383216/1‑150 (MQ=255)
                                                                                                      ccTGCTGGGTGGGCTGGGCGTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGt                                                 <  2:2315038/150‑1 (MQ=255)
                                                                                                       ctgctgGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTc                                                <  2:1723130/150‑1 (MQ=255)
                                                                                                          ctgGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGAt                                             <  1:343516/150‑1 (MQ=255)
                                                                                                           tgGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATc                                            <  1:1676821/150‑1 (MQ=255)
                                                                                                              ttGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGc                                         <  1:353534/150‑1 (MQ=255)
                                                                                                                 gCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGa                                              <  1:1393176/142‑1 (MQ=255)
                                                                                                                 gCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGa                                              >  2:1393176/1‑142 (MQ=255)
                                                                                                                   cTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCa                                           <  2:2285802/143‑1 (MQ=255)
                                                                                                                   cTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCa                                           >  1:2285802/1‑143 (MQ=255)
                                                                                                                        cTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTgg                               <  1:2166038/150‑1 (MQ=255)
                                                                                                                        cTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTgg                               <  1:264154/150‑1 (MQ=255)
                                                                                                                        cTTCTGCGTCCGAGCTGATGTTGCTCCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCGGGTCCTGTTGTTGCTGCTCTTCGCTGCCTTCTTCTTCGCTGCCGTCGATCAGCTCGCCGTTgg                               >  2:709546/1‑150 (MQ=255)
                                                                                                                            tGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCcttc                                                               <  1:480519/114‑1 (MQ=255)
                                                                                                                            tGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCcttc                                                               >  2:480519/1‑114 (MQ=255)
                                                                                                                              cgtcgtAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGtt                                                                                                                   <  1:901589/60‑1 (MQ=255)
                                                                                                                              cGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGtt                                                                                                                   >  2:901589/1‑60 (MQ=255)
                                                                                                                                tCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTgcgc                       >  2:1658934/1‑150 (MQ=255)
                                                                                                                                    aGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCATCTCGCCGTTGGCATTGCGCTGAc                   <  2:2084791/150‑1 (MQ=255)
                                                                                                                                     gCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACg                  >  1:1920212/1‑150 (MQ=255)
                                                                                                                                        gtgacggctgctGTCGCCAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTc               >  2:449153/7‑150 (MQ=255)
                                                                                                                                        gATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCcagca                                                                                                                          <  2:513416/43‑1 (MQ=255)
                                                                                                                                        gATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCcagca                                                                                                                          >  1:513416/1‑43 (MQ=255)
                                                                                                                                        gATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTc               >  2:1802659/1‑150 (MQ=255)
                                                                                                                                        gATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTc               >  2:76355/1‑150 (MQ=255)
                                                                                                                                         aTGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCa              <  2:40131/150‑1 (MQ=255)
                                                                                                                                             ggctgctGTCGCCAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATc                                            >  2:752343/2‑116 (MQ=255)
                                                                                                                                             ggctgctGTCGCCAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATc                                            <  1:752343/115‑1 (MQ=255)
                                                                                                                                                tgctgGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCAc                                                       <  1:1095189/102‑1 (MQ=255)
                                                                                                                                                tgctgGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCAc                                                       >  2:1095189/1‑102 (MQ=255)
                                                                                                                                                tgctgGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGt       >  2:181590/1‑150 (MQ=255)
                                                                                                                                                   tgGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGTCGCTCACGACGGTtgc    >  2:557967/1‑150 (MQ=255)
                                                                                                                                                    ggCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTtgct   <  1:1736585/150‑1 (MQ=255)
                                                                                                                                                     gCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTtgctg  >  2:1474011/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCTTGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTTGCTG  >  NC_002947/2147778‑2148076

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