Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,531,156 G→T 6.6% G118C (GGC→TGC)  PP_2219 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,531,1560GT 6.6% 239.0 / 4.9 76G118C (GGC→TGC) PP_2219hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base T (3/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

CTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGA  >  NC_002947/2531007‑2531296
                                                                                                                                                      |                                                                                                                                            
ctgagctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTATACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAg                                                                                                                                              <  1:2215187/150‑1 (MQ=255)
 tgagctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAgg                                                                                                                                             >  2:363779/1‑150 (MQ=255)
           tGCCGACCGGCGGATTGTGAGCAGGGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGc                                                                                                                                   <  1:931558/150‑1 (MQ=255)
               gACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACAt                                                                                                                               >  1:1834951/1‑150 (MQ=255)
                  cggcggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            >  1:873753/1‑150 (MQ=255)
                    gcggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATc                                                                                                                          <  2:1472755/150‑1 (MQ=255)
                     cggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         <  2:1682540/150‑1 (MQ=255)
                      ggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         >  2:1441551/1‑150 (MQ=255)
                      ggATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCg                                                                                                                         >  2:1443509/1‑150 (MQ=255)
                              aGCAGTTGAAGCAAATTCGCACCCATCTGCAGTTGCTCCACGGGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATCCCGCTGGCTATGCATAGGCATCGACGGCACATGAATCCCCCAAATGc                                                                                                                <  2:995399/150‑1 (MQ=255)
                              aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAATAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGc                                                                                                                <  1:2237773/150‑1 (MQ=255)
                              aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCCGTCGATGCCTATGCATCgccagc                                                                                                                                   >  1:919236/1‑126 (MQ=255)
                              aGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCCGTCGATGCCTATGCATCgccagc                                                                                                                                   <  2:919236/131‑6 (MQ=255)
                               gCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGAGCAGCTTGTTCTTAGCAGAAAGCCGCCTGCGGTGCATAGGCCTCGACGGCACATGAATCGCCCCAATGcc                                                                                                               >  1:877193/1‑149 (MQ=255)
                                      aaGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACAttttt                                                                                                        >  2:165959/1‑150 (MQ=255)
                                      aaGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACAttttt                                                                                                        >  1:2012567/1‑150 (MQ=255)
                                          aaaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCAt                                                                                                    >  2:3029/1‑150 (MQ=255)
                                           aaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATa                                                                                                   <  2:1834951/150‑1 (MQ=255)
                                           aaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATa                                                                                                   <  1:770453/150‑1 (MQ=255)
                                             ttCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATaga                                                                                                 >  1:776048/1‑150 (MQ=255)
                                                  aCCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAt                                                                                            >  1:1582592/1‑150 (MQ=255)
                                                  aCCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAt                                                                                            >  2:2418040/1‑150 (MQ=255)
                                                    ccATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCg                                                                                          >  2:486874/1‑150 (MQ=255)
                                                    ccATTTGCAGTCGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCg                                                                                          >  2:292652/1‑150 (MQ=255)
                                                      aTTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGAt                                                                                        <  1:1716570/150‑1 (MQ=255)
                                                        ttGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  2:776048/150‑1 (MQ=255)
                                                                  tCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACa                                                                            >  2:858042/1‑150 (MQ=255)
                                                                      cgcgGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTg                                                                        >  1:985523/1‑150 (MQ=255)
                                                                      cgcgGCGTGGCGCCGTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCCCAGCAGTCGACGGCACCTGAATCGCCCAAATGCACCTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTg                                                                        <  1:1993741/150‑1 (MQ=255)
                                                                          gcgTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTgg                                                                    <  1:2418040/150‑1 (MQ=255)
                                                                          gcgTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTgg                                                                    <  2:2012567/150‑1 (MQ=255)
                                                                              ggCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCAt                                                                >  2:1016088/1‑150 (MQ=255)
                                                                               gcgcCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATg                                                               <  2:985523/150‑1 (MQ=255)
                                                                                 gcCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGGCAATGAATCGCCCACATGCACATTTTTTCATAGAGACATCGATGACCAAGGGACAGTTGCTGGCCATGTc                                                             >  2:790245/1‑150 (MQ=255)
                                                                                  ccTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACCAATTTTCATAGTGAGAGCTATTATCAAGTAACAGTTTCTGGCCATGTCg                                                            >  2:1592029/1‑150 (MQ=255)
                                                                                   cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  1:1215132/123‑1 (MQ=255)
                                                                                   cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      >  2:1215132/1‑123 (MQ=255)
                                                                                   cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      <  1:1732272/123‑1 (MQ=255)
                                                                                   cTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGa                                                                                      >  2:1732272/1‑123 (MQ=255)
                                                                                      aGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGcc                                                         <  1:1441551/150‑1 (MQ=255)
                                                                                      aGACTGCCGCCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGcc                                                         <  1:1443509/150‑1 (MQ=255)
                                                                                       gACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGa                                                       >  1:1045056/1‑150 (MQ=255)
                                                                                         cTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGc                                                                       >  1:1034962/1‑132 (MQ=255)
                                                                                         cTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGc                                                                       <  2:1034962/132‑1 (MQ=255)
                                                                                          ggccg‑ccCCCTGTAAAAGAGTGCCGGGTGGCGCTTGAACAATTCCGCTGGCGATGCATAGGCAGCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGa                                                    <  2:877193/149‑1 (MQ=255)
                                                                                              g‑ccCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTAc                                                >  1:540495/1‑150 (MQ=255)
                                                                                                  cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            <  2:1825164/71‑1 (MQ=255)
                                                                                                  cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGaa                                                                                                                            >  1:1825164/1‑71 (MQ=255)
                                                                                                  cccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATg                                             >  1:577694/1‑150 (MQ=255)
                                                                                                   ccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGa                                            >  2:302098/1‑150 (MQ=255)
                                                                                                     tGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAAc                                          <  1:858042/150‑1 (MQ=255)
                                                                                                     tGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAAGCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAAc                                          <  1:1022785/150‑1 (MQ=255)
                                                                                                       tAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGc                                        <  1:3029/150‑1 (MQ=255)
                                                                                                          aaGAGTGCCTGTTTGCGCTTGAAGAATTCCGCTTGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGTAACACTTGCTGACCATATCGGCCCATGACTACATTAACGCCGa                                     <  2:873753/150‑1 (MQ=255)
                                                                                                          aaGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAAt                                                                                                                           <  1:534048/64‑1 (MQ=255)
                                                                                                          aaGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAAt                                                                                                                           >  2:534048/1‑64 (MQ=255)
                                                                                                           agagTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCAGCTACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGAc                                    <  1:854694/150‑1 (MQ=255)
                                                                                                               tGCCGGTTGGCGCTTGAAGAATTTCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGc                                >  1:1400676/1‑150 (MQ=255)
                                                                                                               tGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGc                                <  2:540495/150‑1 (MQ=255)
                                                                                                                 ccGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACCTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCtt                              <  1:165959/150‑1 (MQ=255)
                                                                                                                        gcgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttct                       >  1:1894647/1‑150 (MQ=255)
                                                                                                                         cgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttctt                      <  2:1111615/150‑1 (MQ=255)
                                                                                                                         cgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCcttctt                      <  1:1735404/150‑1 (MQ=255)
                                                                                                                              gaagaaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCCAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTg                 >  2:802335/1‑150 (MQ=255)
                                                                                                                                  aaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTTAgct             >  1:638407/1‑150 (MQ=255)
                                                                                                                                     tCCGCTGGCGATGCATAGGCATCGACGGCATATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTg                                                                     >  2:1791236/1‑91 (MQ=255)
                                                                                                                                     tCCGCTGGCGATGCATAGGCATCGACGGCATATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTg                                                                     <  1:1791236/91‑1 (MQ=255)
                                                                                                                                      ccGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGctgctg         <  2:1136567/150‑1 (MQ=255)
                                                                                                                                        gCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCa       <  2:244242/150‑1 (MQ=255)
                                                                                                                                         ccgtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAg      >  2:2118331/5‑150 (MQ=255)
                                                                                                                                           gtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    >  1:1668607/3‑150 (MQ=255)
                                                                                                                                           gtcGATGCCTATGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    <  2:1668607/148‑1 (MQ=255)
                                                                                                                                           ggCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGc    <  2:1894647/150‑1 (MQ=255)
                                                                                                                                             cGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGa  <  1:486874/150‑1 (MQ=255)
                                                                                                                                              gATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGt                                                              >  2:1286163/1‑89 (MQ=255)
                                                                                                                                              gATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGt                                                              <  1:1286163/89‑1 (MQ=255)
                                                                                                                                                   aTAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgct             >  2:1634697/1‑133 (MQ=255)
                                                                                                                                                   aTAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgct             <  1:1634697/133‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
CTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCG‑CCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCAAATGCACATTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGA  >  NC_002947/2531007‑2531296

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