Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,930 A→C 6.6% V878G (GTC→GGC)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,147,9300AC 6.6% 249.4 / 9.6 91V878G (GTC→GGC) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base C (2/4);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

TGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTTGCTG  >  NC_002947/2147783‑2148076
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                            <  2:1489446/144‑1 (MQ=255)
            gCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGcc                                                                                                                                            >  1:1489446/1‑144 (MQ=255)
             ccGCTTCGACCTGCTCGACCGGTTGGGCGGCTTGGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACgg                                                                                                                                     >  2:414736/1‑150 (MQ=255)
              cGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGc                                                                                                                                    >  2:1071913/1‑150 (MQ=255)
              cGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGc                                                                                                                                    >  1:888377/1‑150 (MQ=255)
               gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAca                                                                                                                                                  <  1:478259/135‑1 (MQ=255)
               gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGAca                                                                                                                                                  >  2:478259/1‑135 (MQ=255)
               gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCa                                                                                                                                   >  2:1511690/1‑150 (MQ=255)
               gCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCa                                                                                                                                   >  2:1842241/1‑150 (MQ=255)
                cTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCaa                                                                                                                                  <  1:1842241/150‑1 (MQ=255)
                 ttCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAg                                                                                                                                 >  2:145416/1‑150 (MQ=255)
                 ttCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAg                                                                                                                                 <  1:888007/150‑1 (MQ=255)
                    gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGcc                                                                                                                              <  2:927780/150‑1 (MQ=255)
                    gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGc                                                                                                                               <  2:1920008/149‑1 (MQ=255)
                    gACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGc                                                                                                                               >  1:1920008/1‑149 (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                                                                                                                                        <  1:1247058/126‑1 (MQ=255)
                                  gTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCa                                                                                                                                        >  2:1247058/1‑126 (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                                                                                                  <  1:1625185/150‑1 (MQ=255)
                                                cTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTg                                                                                                  <  2:1970095/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                                                                                            >  2:1436028/1‑150 (MQ=255)
                                                      tCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTgg                                                                                            >  1:1786640/1‑150 (MQ=255)
                                                      tCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTgg                                                                                            >  1:1049629/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                                                                  <  2:1956437/150‑1 (MQ=255)
                                                                                ggcggTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGcc                                                                  <  1:1511690/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                                           >  1:2285802/1‑143 (MQ=255)
                                                                                                              cTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCa                                           <  2:2285802/143‑1 (MQ=255)
                                                                                                                   cTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTgg                               <  1:2166038/150‑1 (MQ=255)
                                                                                                                   cTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTgg                               <  1:264154/150‑1 (MQ=255)
                                                                                                                   cTTCTGCGTCCGAGCTGATGTTGCTCCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCGGGTCCTGTTGTTGCTGCTCTTCGCTGCCTTCTTCTTCGCTGCCGTCGATCAGCTCGCCGTTgg                               >  2:709546/1‑150 (MQ=255)
                                                                                                                       tGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCcttc                                                               >  2:480519/1‑114 (MQ=255)
                                                                                                                       tGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCcttc                                                               <  1:480519/114‑1 (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                                                                                                                          >  1:513416/1‑43 (MQ=255)
                                                                                                                                   gATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCcagca                                                                                                                          <  2:513416/43‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
TGAGCTTTGGTGGCCGCTTCGACCTGCTCGACCGGTTGGGCGGCTTCGCTGTTGTCGGTTTCTGCAACGGCAGCGACTTCGGCGGTAGCACGCTCGGCCTGCTGGTTGGCCTGGGCTTCTGCGTCGGAGCTGATGTTGCTGCTGGCGACAGCAGCCGTCACGGCAAGGCCAGCAGCCAGTTCGGCACCCAGCTCGGTGGCCTGGTGCTGTTGTGGCTGCTCTTCGCTGCCTTCTTCTTCACTGCCGTCGATCAGCTCGCCGTTGGCATTGCGCTGACGCTCACGACGGTTGCTG  >  NC_002947/2147783‑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: ‑