Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,445,251 C→G 5.2% G117R (GGG→CGG) ‡ nicE ← maleate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,445,2510CG 5.2% 232.4 / 3.6 77G117R (GGG→CGG) ‡nicEmaleate isomerase
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base G (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGC  >  NC_002947/4445105‑4445399
                                                                                                                                                  |                                                                                                                                                     
ggggCTTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGGAGTTCACCACCCTCTCGGTCGGCGGCTTCATGTACGGGGCGACGATGGCTATGCGCTTGGCCTGCCTCACCTGGAAGGCCTCGGCCAGCGCCCggg                                                                                                                                                     >  8:412520/1‑149 (MQ=255)
ggggATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCgg                                                                                                                                                     >  6:307891/1‑149 (MQ=255)
ggggATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCgg                                                                                                                                                     >  8:385959/1‑149 (MQ=255)
 gtgATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGgc                                                                                                                                                    <  2:103175/147‑1 (MQ=255)
 gggATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGgc                                                                                                                                                    <  2:97220/149‑1 (MQ=255)
   gATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGgcgc                                                                                                                                                  <  5:87572/149‑1 (MQ=255)
       tCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGt                                                                                                                                              >  2:4527/1‑149 (MQ=255)
       tCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACACGCAAGGCCTCGACCAGCGCCGGGGCGCTgg                                                                                                                                               <  5:395249/148‑1 (MQ=255)
        ccAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGt                                                                                                                                              <  2:127594/148‑1 (MQ=255)
        ccAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGt                                                                                                                                              <  2:127598/148‑1 (MQ=255)
           gcgcgcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTGGtcgag                                                                                                                                          <  5:358707/149‑5 (MQ=255)
                 gcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGt                                                                                                                                    >  1:29774/1‑149 (MQ=255)
                 gcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGt                                                                                                                                    >  2:260363/1‑149 (MQ=255)
                   cAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑Gc                                                                                                                                 >  8:103531/1‑149 (MQ=255)
                    aGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑Gc                                                                                                                                 >  6:35316/1‑148 (MQ=255)
                        cTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTg                                                                                                                            >  8:402063/1‑149 (MQ=255)
                           cACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCATCTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAAGGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGcat                                                                                                                         <  5:300194/149‑1 (MQ=255)
                           cACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCATGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGc                                                                                                                           <  4:310706/147‑1 (MQ=255)
                              cTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCATCGGCTTCATGTACGGCGCGACCATGGCAATCCGCTTGGCCTTCATCACATGCAAGGCCTCGACCACCGCCCCGCCGCTGGTGATCACCGCT‑GCCAGTGcatcat                                                                                                                      >  4:97219/1‑149 (MQ=255)
                              cTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGcatcat                                                                                                                      >  6:13274/1‑149 (MQ=255)
                                  taGCCTTGCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGAGGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACGAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTc                                                                                                                  <  3:97219/148‑1 (MQ=255)
                                  aaGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTc                                                                                                                  <  1:250850/149‑1 (MQ=255)
                                   aGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTcg                                                                                                                 <  7:18423/149‑1 (MQ=255)
                                      cTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTcggcg                                                                                                              <  7:402063/149‑1 (MQ=255)
                                       ttCCTCGCGGATGTAGTTCACCACCAACTCGGACAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGAGCTGGTGATCACAGGT‑GCCAGTGCATCATTGTcggcgg                                                                                                             >  8:194508/1‑149 (MQ=255)
                                         ccTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGt                                                                                                            <  6:120049/148‑1 (MQ=255)
                                          cTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTg                                                                                                           <  5:35316/148‑1 (MQ=255)
                                           tCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCAGT‑GCCAGTGCATCATTGTCGGCGGTGGc                                                                                                         >  7:205530/1‑149 (MQ=255)
                                             gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTCGCCAGTGCATCATTGTCGGCGGTGGc                                                                                                         <  1:4527/148‑1 (MQ=255)
                                                gATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCctgct                                                                                                    <  3:367889/149‑1 (MQ=255)
                                                  tGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCctgctg                                                                                                   >  6:292116/1‑148 (MQ=255)
                                                    tAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAg                                                                                                >  5:322721/1‑149 (MQ=255)
                                                      gTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAgg                                                                                               <  2:45643/148‑1 (MQ=255)
                                                          accaccAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtt                                                                                          >  8:189404/1‑149 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                        >  5:141090/1‑149 (MQ=255)
                                                              ccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCa                                                                                      <  1:418335/149‑1 (MQ=255)
                                                              ccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCACtttt                                                                                        >  4:338632/1‑147 (MQ=255)
                                                               caACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                        >  7:71707/1‑146 (MQ=255)
                                                                  cTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGAtt                                                                                  <  1:260363/149‑1 (MQ=255)
                                                                   tCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTg                                                                                 >  2:272215/1‑149 (MQ=255)
                                                                     ggTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGc                                                                               <  7:385959/149‑1 (MQ=255)
                                                                      gTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCg                                                                              <  7:103531/149‑1 (MQ=255)
                                                                       tCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                             <  6:141090/149‑1 (MQ=255)
                                                                          gCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGAAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATAACCGGT‑GCAAGTGCATCATTGTCGGAGGTGGCCTGCTGCAGGCGCTTATCACATTGGCTGTGa                                                                          >  7:80204/1‑149 (MQ=255)
                                                                                      aCGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGcc                                                                                                        >  5:223007/1‑107 (MQ=255)
                                                                                      aCGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGcc                                                                                                        <  6:223007/107‑1 (MQ=255)
                                                                                         gcgcGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATg                                                           >  3:393058/1‑149 (MQ=255)
                                                                                          cgcgACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATgg                                                          >  8:138397/1‑149 (MQ=255)
                                                                                             gacgaTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       <  1:417493/149‑1 (MQ=255)
                                                                                                 aTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGAt                                                   >  7:250556/1‑149 (MQ=255)
                                                                                                  tGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       >  1:422821/1‑144 (MQ=255)
                                                                                                  tGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       <  2:422821/144‑1 (MQ=255)
                                                                                                   ggaaaTGGGCTTGGCCTTCATCACATGCAATGCCTCGACCAGCGCGCCGGCGCTGGTGATCACCGGT‑GCTAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATc                                                  <  8:80204/145‑1 (MQ=255)
                                                                                                    gCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGTCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGc                                                <  5:307891/149‑1 (MQ=255)
                                                                                                      aaTGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCca                                              <  2:29774/149‑1 (MQ=255)
                                                                                                            cTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAagg                                         >  5:398468/1‑148 (MQ=255)
                                                                                                             ttGGCCTTCCTCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAaggca                                       <  1:170742/149‑1 (MQ=255)
                                                                                                               ggCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCAGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAaggcag                                      <  8:205530/148‑1 (MQ=255)
                                                                                                                  cTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCAt                                  <  6:322721/149‑1 (MQ=255)
                                                                                                                   ttCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAggcaggca                                   >  5:401955/1‑147 (MQ=255)
                                                                                                                    tcatcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCGATGGCCa                                                       <  2:127734/126‑1 (MQ=255)
                                                                                                                    tcatcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       >  1:127734/1‑126 (MQ=255)
                                                                                                                      atcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGcc                              >  7:206562/1‑149 (MQ=255)
                                                                                                                      atcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGAATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGcc                              >  7:272073/1‑149 (MQ=255)
                                                                                                                         acaTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCa                            <  8:206562/148‑1 (MQ=255)
                                                                                                                         acaTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCa                            <  7:138397/148‑1 (MQ=255)
                                                                                                                         acaTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAg                           >  3:64302/1‑149 (MQ=255)
                                                                                                                          caTGCAAGGCCTCGACCAGCGCCGCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGACAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGc                          <  2:125407/149‑1 (MQ=255)
                                                                                                                           aTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGca                         >  1:322574/1‑149 (MQ=255)
                                                                                                                            tGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGAGGTGATAGCCCAGCCCCAAGGCCATGATCGCCACAAGGCAGGCATAGCCCAGca                         >  5:2559/1‑148 (MQ=255)
                                                                                                                               aaGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATc                     >  5:438846/1‑149 (MQ=255)
                                                                                                                                aGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCg                    <  6:401955/149‑1 (MQ=255)
                                                                                                                                  gCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGAc                  <  8:343247/149‑1 (MQ=255)
                                                                                                                                       gACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTg              <  8:272073/148‑1 (MQ=255)
                                                                                                                                       ctccAGCGCCGGGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGTTTGGCGGTGATAGCGCAGCCCCATGGGAATGATCGGCACAAGGCAGGCATAGCCCATCCAATCGACGTTgg             >  6:363661/3‑149 (MQ=255)
                                                                                                                                        aCCAGCGCCGGGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGc            >  3:251117/1‑149 (MQ=255)
                                                                                                                                        aCCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGc            <  5:13274/149‑1 (MQ=255)
                                                                                                                                             cgcCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCg        >  7:404922/1‑148 (MQ=255)
                                                                                                                                                cccGGCGCTCGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCGTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATACCCCAGCACAACGACCTTCGCGTCGGAc     >  8:123891/1‑148 (MQ=255)
                                                                                                                                                 ccGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAg   >  3:426765/1‑149 (MQ=255)
                                                                                                                                                  cGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGTCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGc  <  8:250556/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGT‑GCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGC  >  NC_002947/4445105‑4445399

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