Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,540,539 A→T 5.4% intergenic (+59/+62) PP_4871 → / ← PP_4872 lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,540,5390AT 5.4% 258.4 / 3.3 86intergenic (+59/+62)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (51/27);  new base T (2/3);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTA  >  NC_002947/5540390‑5540688
                                                                                                                                                     |                                                                                                                                                     
gctgGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACaaa                                                                                                                                                       <  2:3821622/150‑1 (MQ=255)
    gggCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCAAAGGcc                                                                                                                                                   >  1:1303695/1‑149 (MQ=255)
        agcagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgg                                                                                                                                               <  2:1223932/150‑1 (MQ=255)
          cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                             <  1:4994755/150‑1 (MQ=255)
          cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                             <  1:4297696/150‑1 (MQ=255)
          cagcCGGACGAGATGCTGGAAGGGAGGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTAGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                             <  1:4431753/150‑1 (MQ=255)
            gcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGa                                                                                                                                           <  1:1852753/150‑1 (MQ=255)
                 aCGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGc                                                                                                                                      >  1:2918586/1‑150 (MQ=255)
                 aCGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCAGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCACTGGGGGAGCGGc                                                                                                                                      <  2:3019761/150‑1 (MQ=255)
                  cGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGcc                                                                                                                                     <  2:4191850/150‑1 (MQ=38)
                    aGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                                                                   <  1:4257435/150‑1 (MQ=255)
                     gctgctGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtg                                                                                                                                  >  2:2006642/1‑150 (MQ=255)
                      ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTtgt                                                                                                                                 <  1:1413616/150‑1 (MQ=38)
                          tgtaaGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcg                                                                                                                             >  1:383725/4‑150 (MQ=255)
                          tgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcg                                                                                                                             <  2:3263179/150‑1 (MQ=255)
                           ggACGGTATGGACAGCTTTGTGGCGCCGGTCGCGCCTACGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCGGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                            <  1:1262082/150‑1 (MQ=255)
                                ggATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGAGTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTTTGTCGCGATgg                                                                                                                       >  1:2484817/1‑150 (MQ=255)
                                    ggACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTg                                                                                                                   >  1:2489595/1‑150 (MQ=255)
                                     gACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGc                                                                                                                  >  1:3458208/1‑150 (MQ=255)
                                      aCAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCa                                                                                                                 >  2:2936713/1‑150 (MQ=255)
                                        cgCTTTGTGGCGCCGGTCGCGCCTACGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCATTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCaaa                                                                                                               <  2:879784/149‑1 (MQ=255)
                                        cgCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAACCACCATCACTGAAGCAGCGACTTTTGGGACCGCTACGCGGTCCTTTCGCGCCACCCGGCCGCGCCCACAAAGGCAATGGGGGAGGGGCGTTGGGTCGCGAGGGGCTGcttt                                                                                                               <  2:2753616/149‑4 (MQ=255)
                                           tttGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAagca                                                                                                            >  1:1830718/1‑150 (MQ=255)
                                                   gcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAAt                                                                                                    >  2:2954460/1‑150 (MQ=255)
                                                      ggTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGGCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGGGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGt                                                                                                 <  2:3858080/150‑1 (MQ=255)
                                                               tAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGTCAATGGGGGTGCGGCCTTGTGTCGCGATGGGCTGAAAAGCAGCCCCAATAGTCGATCAGtt                                                                                        >  1:246167/1‑150 (MQ=255)
                                                               tAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGtt                                                                                        <  2:2118585/150‑1 (MQ=255)
                                                                      gTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCgg                                                                                 <  2:1008797/150‑1 (MQ=255)
                                                                       tCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACACGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGa                                                                                <  2:3693629/150‑1 (MQ=255)
                                                                        cGACACGCCGCCTGACTGAAGCTGTTACTTTTGGGACCGCTGCGGGGCCCTTTCGCCACACACGGCCGCTCCCACAACGGCATTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGaa                                                                               <  1:3445609/150‑1 (MQ=37)
                                                                            acgcagcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGAc                                                                           <  1:4648292/148‑1 (MQ=255)
                                                                             agcagcCTCACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACa                                                                          <  2:3357660/150‑1 (MQ=255)
                                                                                    tgactgaAGCTGAGACCTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGAc                                                                   >  2:2604481/1‑150 (MQ=255)
                                                                                      actgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGAcgc                                                                 <  2:4564819/150‑1 (MQ=255)
                                                                                       ctgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGa                                                                                                                                           >  1:936444/1‑75 (MQ=255)
                                                                                       ctgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGa                                                                                                                                           <  2:936444/75‑1 (MQ=255)
                                                                                        ttaaGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGCCACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCgg                                                               <  1:1846397/148‑1 (MQ=255)
                                                                                        tgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGTCACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCACCCCCAATAGTCGGTCAGTTAACGCCGAATGACAGGCTGACGCgg                                                               >  1:2255962/1‑150 (MQ=255)
                                                                                        tgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCCATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCgg                                                               >  2:38774/1‑150 (MQ=255)
                                                                                        tgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTTCAAAGCAGCCCCAAAATTAGATCAGTTAAACCGGAATGACAGGCTGAAGCgg                                                               >  1:1332339/1‑150 (MQ=38)
                                                                                            gCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAAGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGgttc                                                           >  1:3961980/1‑150 (MQ=255)
                                                                                             cTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACCCAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttct                                                          <  1:4136258/150‑1 (MQ=255)
                                                                                              tGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGAAATGCAGACCCGGcttcat                                                         >  1:2693607/1‑148 (MQ=39)
                                                                                               gagaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAAAAGTCGATCAGTTAACGCGGAAAGACAGGCTGACGCGGCTTCttg                                                        >  1:2007833/1‑150 (MQ=255)
                                                                                               gagaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCCATAGTCGATCAGGTAACGCGGAATGACAGGCTGACGCGGCGTCttg                                                        >  1:610305/1‑150 (MQ=255)
                                                                                                    ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCGTGTTGcc                                                   >  2:2984917/1‑150 (MQ=255)
                                                                                                    ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCCATAGTCGATCCGTTAACGCGGAATGACAGGCTGACGCGGCGTCTTGTTGcc                                                   >  1:1953934/1‑150 (MQ=255)
                                                                                                    ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATCGTCGATCAGGTAACGCGGCATGACAGGCTGACGCGGCGTCTTGTTGcc                                                   >  1:4665476/1‑150 (MQ=39)
                                                                                                    ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAGGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGcc                                                   >  2:1872580/1‑150 (MQ=255)
                                                                                                    ttttGGGACCGCTTCGCGGACCTTTCGCCACACACGGCCGCTCCCACAATGGCAATGGGGGAGCGCCCTTGTGACGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGcc                                                   <  1:1421476/150‑1 (MQ=255)
                                                                                                     tttGGGACCGCTTCGCGGTCCTTCCGCGACACAAGGCCGCTCCCACAACGGCAATGGGGGAGGGGCCTTGTGGCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGccc                                                  <  1:3016358/150‑1 (MQ=255)
                                                                                                      aaaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATTGGCTGCAAAGCAGCCCCAAAAGTCGATAAGTTAAAGCGGAAAGACAGGCTGAAGCGGCGTCTTTGTGCCCt                                                 >  1:117061/4‑150 (MQ=255)
                                                                                                       tGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCtt                                                <  2:2898357/150‑1 (MQ=255)
                                                                                                         ggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTACCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCCATAGGCGCTCAGTTAACGCGGAATGACAGGCTGACGCGGGTTCGTGTTGCCCGTGc                                              >  1:4696275/1‑150 (MQ=38)
                                                                                                          gACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGcc                                             >  2:3064753/1‑150 (MQ=255)
                                                                                                                   gcgGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCt                                    >  2:626433/1‑150 (MQ=255)
                                                                                                                    cgGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTc                                   <  2:4649699/150‑1 (MQ=255)
                                                                                                                    cgGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTc                                   <  1:805489/150‑1 (MQ=255)
                                                                                                                      gTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAAGGGGAGCGGCCTTGTGTCGCGATGGGCAGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGAAAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCATGCTCGcc                                >  2:1388893/1‑150 (MQ=39)
                                                                                                                          tttCGCGACACAAGGCCGCTCCCACAAAGGCAAGGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAACGCAGCCCCCCTCGGCGCTCAGGTAACGCGGAATGACAGGCTGACGCGGCTTCTGGGTGCCCTTGCCACCCTTGCGCGCctct                             >  2:4314257/1‑150 (MQ=37)
                                                                                                                           ttCGCGACACAAGGCCGCTACCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCGTCTTGTTGCCCTTGCCAACCTTGCTCGCctcgt                            >  2:448482/1‑148 (MQ=255)
                                                                                                                           atCGTGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAAGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCtt                            >  1:11713/2‑150 (MQ=255)
                                                                                                                            tCGTGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCttt                           <  2:11713/150‑1 (MQ=255)
                                                                                                                              gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                         >  2:439737/1‑150 (MQ=255)
                                                                                                                              gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGAGCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                         >  2:3999679/1‑150 (MQ=255)
                                                                                                                                     aaGGCCGCTCCCACAACGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                  <  1:4843410/150‑1 (MQ=255)
                                                                                                                                     aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGACTGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                  >  1:2846408/1‑150 (MQ=255)
                                                                                                                                     aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                  >  1:411358/1‑150 (MQ=255)
                                                                                                                                     aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                  >  1:2846497/1‑150 (MQ=255)
                                                                                                                                      aGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                 >  1:400179/1‑150 (MQ=255)
                                                                                                                                      aGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                 >  2:4355000/1‑150 (MQ=255)
                                                                                                                                       ggCCGCTCCCACAAAGGCAATGGGGGTGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg                <  1:626433/150‑1 (MQ=255)
                                                                                                                                        gCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGc               >  2:1732257/1‑150 (MQ=38)
                                                                                                                                        gCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTCTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGc               >  2:3155387/1‑150 (MQ=38)
                                                                                                                                         ccGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCa              <  2:283228/150‑1 (MQ=255)
                                                                                                                                          cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAg             >  2:3107059/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGc            >  1:1018820/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGc            >  1:2097763/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGc            >  1:3314481/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGc            >  1:1018772/1‑150 (MQ=255)
                                                                                                                                            cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc           >  1:1512459/1‑150 (MQ=255)
                                                                                                                                            cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc           >  2:2747240/1‑150 (MQ=255)
                                                                                                                                            cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc           >  2:119747/1‑150 (MQ=255)
                                                                                                                                             tCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCct          >  2:2432952/1‑150 (MQ=255)
                                                                                                                                               ccACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgc        >  2:4382699/1‑150 (MQ=255)
                                                                                                                                                 acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg      >  2:2248198/1‑150 (MQ=255)
                                                                                                                                                 acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAAAAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCGTATTGTTGCCCGTGCCACCCTTGCTCGCCTCTTTGCGCTCCATGGCAGCatgctg      >  1:1771812/1‑150 (MQ=38)
                                                                                                                                                   aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt    >  2:3668227/1‑150 (MQ=255)
                                                                                                                                                   aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt    >  1:3592122/1‑150 (MQ=255)
                                                                                                                                                   aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt    >  1:4768326/1‑150 (MQ=255)
                                                                                                                                                   aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt    >  1:4768397/1‑150 (MQ=255)
                                                                                                                                                     aGGCAATGGGGGAGCGCCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTa  <  1:461909/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTA  >  NC_002947/5540390‑5540688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: