Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,540,544 A→T 6.1% intergenic (+64/+57) 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,5440AT 6.1% 263.2 / 5.9 92intergenic (+64/+57)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (55/30);  new base T (2/4);  total (57/35)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

AGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTACGCG  >  NC_002947/5540398‑5540692
                                                                                                                                                  |                                                                                                                                                    
agcagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgg                                                                                                                                                   <  2:1223932/150‑1 (MQ=255)
  cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                                 <  1:4297696/150‑1 (MQ=255)
  cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                                 <  1:4994755/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                                                                                                                                               <  2:936444/75‑1 (MQ=255)
                                                                               ctgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGa                                                                                                                                               >  1:936444/1‑75 (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                                       <  1:805489/150‑1 (MQ=255)
                                                                                                            cgGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTc                                       <  2:4649699/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                     >  2:4355000/1‑150 (MQ=255)
                                                                                                                              aGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                     >  1:400179/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:1018772/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)
                                                                                                                                    cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc               >  2:119747/1‑150 (MQ=255)
                                                                                                                                    cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc               >  1:1512459/1‑150 (MQ=255)
                                                                                                                                    cTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGcc               >  2:2747240/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:4768397/1‑150 (MQ=255)
                                                                                                                                           aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt        >  1:4768326/1‑150 (MQ=255)
                                                                                                                                           aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt        >  1:3592122/1‑150 (MQ=255)
                                                                                                                                             aGGCAATGGGGGAGCGCCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTa      <  1:461909/150‑1 (MQ=255)
                                                                                                                                               gCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAcg    >  1:2580550/1‑150 (MQ=255)
                                                                                                                                                cAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAcgc   <  2:1303695/150‑1 (MQ=255)
                                                                                                                                                 aaTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATTACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAcgcg  >  2:1319785/1‑150 (MQ=255)
                                                                                                                                                 aaTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAcgcg  >  2:2502695/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTACGCG  >  NC_002947/5540398‑5540692

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: