Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,601,876 C→A 5.6% D553Y (GAC→TAC)  PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,601,8760CA 5.6% 227.9 / 2.4 72D553Y (GAC→TAC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (2/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTC  >  NC_002947/5601728‑5602023
                                                                                                                                                     |                                                                                                                                                   
gCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTcc                                                                                                                                                    >  2:134988/1‑150 (MQ=255)
 cACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCa                                                                                                                                                   <  1:113558/150‑1 (MQ=255)
  aCTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAg                                                                                                                                                  >  1:1695615/1‑150 (MQ=255)
   cTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGc                                                                                                                                                 <  1:1166870/150‑1 (MQ=255)
    tGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCg                                                                                                                                                >  1:1613029/1‑150 (MQ=255)
               aCATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGTTGCCTTGa                                                                                                                                     >  1:1568359/1‑150 (MQ=255)
                cATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAg                                                                                                                                    <  1:1920359/150‑1 (MQ=255)
                  tcgtcgCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCa                                                                                                                                  >  2:1864295/1‑150 (MQ=255)
                       gcccgcccCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTg                                                                                                                             >  1:1191698/1‑150 (MQ=255)
                            ccccgccGTGATAGGTGTTGCTGCCCCCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGa                                                                                                                        <  2:606001/150‑1 (MQ=255)
                                gccgTGATAGTTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                                    <  2:109747/150‑1 (MQ=255)
                                gccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                                    <  1:646265/150‑1 (MQ=255)
                                gccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                                    <  2:492585/150‑1 (MQ=255)
                                 ccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTTGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                                   <  1:383141/150‑1 (MQ=255)
                                  cgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGTTGCCTTGAGCACTCTGGCTGATGTAgg                                                                                                                  <  2:1568359/150‑1 (MQ=255)
                                      aTAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATTATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCt                                                                                                              >  1:847347/1‑150 (MQ=255)
                                       tAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCtt                                                                                                             <  2:656981/150‑1 (MQ=255)
                                        aGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTcc                                                                                                                                                    >  1:1906318/1‑110 (MQ=255)
                                        aGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTcc                                                                                                                                                    <  2:1906318/110‑1 (MQ=255)
                                         ggTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCg                                                                                                           <  1:1864295/150‑1 (MQ=255)
                                          gtgtTGCTGCCACCGCCACCAACCAGCCACCTGCCTTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                                          <  2:15284/150‑1 (MQ=255)
                                          gtgtTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                                          <  1:1301461/150‑1 (MQ=255)
                                          gtgtTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                                          >  1:899151/1‑150 (MQ=255)
                                               gctgcCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGt                                                                                                     >  1:1580285/1‑150 (MQ=255)
                                                   ccaccgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGc                                                                                                 >  1:257157/1‑150 (MQ=255)
                                                     accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGccc                                                                                               >  2:867216/1‑150 (MQ=255)
                                                     accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGccc                                                                                               >  2:1671126/1‑150 (MQ=255)
                                                         ccaccaACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCtt                                                                                           <  2:1049378/150‑1 (MQ=255)
                                                              aaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                                   >  2:772718/1‑121 (MQ=255)
                                                              aaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                                   <  1:772718/121‑1 (MQ=255)
                                                              aaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGa                                                                                      >  2:942918/1‑150 (MQ=255)
                                                               accaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                                   >  1:418925/1‑120 (MQ=255)
                                                               accaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAg                                                                                                                   <  2:418925/120‑1 (MQ=255)
                                                                        cTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTaga                                                                            <  2:937046/150‑1 (MQ=255)
                                                                           cccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTg                                                                         >  1:668395/1‑150 (MQ=255)
                                                                           cccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTCCCCACTTGGGTACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTg                                                                         >  2:59385/1‑150 (MQ=255)
                                                                            ccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGAc                                                                                        <  1:6188/134‑1 (MQ=255)
                                                                            ccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGAc                                                                                        >  2:6188/1‑134 (MQ=255)
                                                                              tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTc                                                                                                                                                     >  1:1596919/1‑71 (MQ=255)
                                                                              tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTc                                                                                                                                                     <  2:1596919/71‑1 (MQ=255)
                                                                              tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                                          >  2:696702/1‑114 (MQ=255)
                                                                              tctGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCgg                                                                                                          <  1:696702/114‑1 (MQ=255)
                                                                                     ccGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGATGATGTAAGAGCTTCGGATAGAGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcctcct                                                               >  1:918464/1‑150 (MQ=255)
                                                                                      cgcgGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTg                                                              <  2:1205683/150‑1 (MQ=255)
                                                                                      cgcgCTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTg                                                              <  2:1763499/150‑1 (MQ=255)
                                                                                       gcgGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTAAGACGAGCCCCACGGGGTTCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCCCCCTCCTgg                                                             <  2:50896/150‑1 (MQ=255)
                                                                                         ggTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGc                                                           <  1:134988/150‑1 (MQ=255)
                                                                                           tCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCt                                                         <  1:942918/150‑1 (MQ=255)
                                                                                             aaGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGaa                                                                                     <  2:299289/120‑1 (MQ=255)
                                                                                             aaGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGaa                                                                                     >  1:299289/1‑120 (MQ=255)
                                                                                                  atcatcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCg                                                  >  2:986240/1‑150 (MQ=255)
                                                                                                  atcatcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCg                                                  >  1:896640/1‑150 (MQ=255)
                                                                                                  atcatcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCg                                                  >  1:896657/1‑150 (MQ=255)
                                                                                                   tcatcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTTAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGa                                                 <  1:1866004/150‑1 (MQ=255)
                                                                                                    catcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGAt                                                <  1:1671126/150‑1 (MQ=255)
                                                                                                     atcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGAtt                                               >  1:1456406/1‑150 (MQ=255)
                                                                                                      tcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTa                                              >  2:1766623/1‑150 (MQ=255)
                                                                                                           cgcgCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGAc                                         <  1:867216/150‑1 (MQ=255)
                                                                                                           cgcgCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGAc                                         <  1:986240/150‑1 (MQ=255)
                                                                                                             cgcTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACgg                                       <  2:257157/150‑1 (MQ=255)
                                                                                                              gcTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGt                                      <  2:1613029/150‑1 (MQ=255)
                                                                                                                tGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCATCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAAGCCCGTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGCTTACAGACGGTAt                                    >  2:501620/1‑150 (MQ=255)
                                                                                                                     tAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAg                                                                                                                                    <  2:491393/49‑1 (MQ=255)
                                                                                                                     tAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAg                                                                                                                                    >  1:491393/1‑49 (MQ=255)
                                                                                                                         gCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGaa                           <  2:1580285/150‑1 (MQ=255)
                                                                                                                         gCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATACCGACGGTATCGGCGAGaaa                          >  1:478577/1‑150 (MQ=255)
                                                                                                                              tggaCGCCGCAGCGCTCGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                                         <  2:247449/80‑1 (MQ=37)
                                                                                                                              tggaCGCCGCAGCGCTCGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                                         >  1:247449/4‑83 (MQ=37)
                                                                                                                              tggaCGCCGCAGCGCTCGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTc                      >  2:861488/4‑150 (MQ=255)
                                                                                                                                   gCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAg                                                                                >  1:1366688/1‑87 (MQ=255)
                                                                                                                                   gCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAg                                                                                <  2:1366688/87‑1 (MQ=255)
                                                                                                                                   gCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGc                 <  2:1601915/150‑1 (MQ=255)
                                                                                                                                            cTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGc                                >  2:649598/1‑127 (MQ=255)
                                                                                                                                            cTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGc                                <  1:649598/127‑1 (MQ=255)
                                                                                                                                               cggcgTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTccc      <  1:861488/150‑1 (MQ=255)
                                                                                                                                                  cgTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTc  <  2:478577/150‑1 (MQ=255)
                                                                                                                                                   gTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTaga                                                                            <  2:74717/76‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTC  >  NC_002947/5601728‑5602023

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