Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,052,430 C→A 6.1% T23T (ACG→ACT)  PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,052,4300CA 6.1% 212.5 / 3.3 66T23T (ACG→ACT) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ACACCACCAACATCCCTGGGGCCGACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAGGCCATCG  >  NC_002947/4052284‑4052578
                                                                                                                                                  |                                                                                                                                                    
acaccacCAACATCCCTGGAGCCGACAGTGGCGAACAGGATGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCAGCTCACCGTg                                                                                                                                                   >  2:20126/1‑150 (MQ=255)
     accaACATCCCTGGGGCCGACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCaaa                                                                                                                                              <  1:287555/150‑1 (MQ=255)
                       gACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAACACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCcag                                                                                                                            >  2:868334/1‑150 (MQ=255)
                       gACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTg                                                                                                                                                   <  2:1582361/127‑1 (MQ=255)
                       gACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTg                                                                                                                                                   >  1:1582361/1‑127 (MQ=255)
                        aCAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCcagc                                                                                                                           >  2:146137/1‑150 (MQ=255)
                         cAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCcagca                                                                                                                          <  2:1001954/150‑1 (MQ=255)
                              gCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAgcgcg                                                                                                                     >  1:319414/1‑150 (MQ=255)
                                gAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGAt                                                                                                                                        >  1:518143/1‑129 (MQ=255)
                                gAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGAt                                                                                                                                        <  2:518143/129‑1 (MQ=255)
                                 aaCAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAGAGCCGATCAGCCGGGACAGCAGCGCGCGa                                                                                                                  >  1:743445/1‑150 (MQ=255)
                                      gTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGc                                                                                                             >  1:1058868/1‑150 (MQ=255)
                                      gTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACAGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGc                                                                                                             >  2:1292070/1‑150 (MQ=255)
                                      gTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACAGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGc                                                                                                             >  2:1292158/1‑150 (MQ=255)
                                       ttGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCa                                                                                                            >  2:1262821/1‑150 (MQ=255)
                                       ttGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACAGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGcc                                                                                                            <  1:1292070/150‑2 (MQ=255)
                                       ttGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACAGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCa                                                                                                            <  1:1292158/150‑1 (MQ=255)
                                           ccTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAAtt                                                                                                        >  2:1165733/1‑150 (MQ=255)
                                              tGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCAc                                                                                                     <  2:326149/150‑1 (MQ=255)
                                               ggTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACt                                                                                                    >  2:585232/1‑150 (MQ=255)
                                                          ggTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGc                                                                                         <  2:743445/150‑1 (MQ=255)
                                                          ggTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCAAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGc                                                                                         <  1:516542/150‑1 (MQ=255)
                                                           gTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCa                                                                                        >  2:1400360/1‑150 (MQ=255)
                                                              aTGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCActc                                                                                     >  1:23952/1‑148 (MQ=255)
                                                                ggAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCa                                                                                   >  2:283629/1‑150 (MQ=255)
                                                                 gAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCAt                                                                                  >  2:1009723/1‑150 (MQ=255)
                                                                 gAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCAt                                                                                  >  2:972793/1‑150 (MQ=255)
                                                                    cAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGc                                                                               <  1:1009723/150‑1 (MQ=255)
                                                                      gCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAg                                                                                >  1:375538/1‑147 (MQ=255)
                                                                      gCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAg                                                                                <  2:375538/147‑1 (MQ=255)
                                                                      gCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCtt                                                                             <  2:1058868/150‑1 (MQ=255)
                                                                        tCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTgg                                                                           >  2:1358714/1‑150 (MQ=255)
                                                                             cTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCgcgg                                                                      >  1:1541277/1‑150 (MQ=255)
                                                                              tCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCgcggc                                                                     <  1:531249/150‑1 (MQ=255)
                                                                              tCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCgcggc                                                                     <  1:146137/150‑1 (MQ=255)
                                                                                cACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGcggcgg                                                                   >  2:297297/1‑150 (MQ=255)
                                                                                     ccAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCggg                                                              >  1:737513/1‑150 (MQ=255)
                                                                                        gCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTg                                                           >  1:1025153/1‑150 (MQ=255)
                                                                                         cTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGa                                                          <  2:791704/150‑1 (MQ=255)
                                                                                               taGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGc                                                                                                                                            <  1:1504238/61‑1 (MQ=255)
                                                                                               cAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGc                                                                                                                                            >  2:1504238/1‑62 (MQ=255)
                                                                                               cAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTg                                                                                                               >  1:199780/1‑91 (MQ=255)
                                                                                               cAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTg                                                                                                               <  2:199780/91‑1 (MQ=255)
                                                                                                 ggTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCCGCGCGCGACTGGCAAATTTTGACATCGCCCGACACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCgg                                                  <  1:283629/150‑1 (MQ=255)
                                                                                                 ggTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCgg                                                  <  1:1400360/150‑1 (MQ=255)
                                                                                                 ggTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCgg                                                  <  1:585232/150‑1 (MQ=255)
                                                                                                 ggTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCgg                                                  <  1:1165733/150‑1 (MQ=255)
                                                                                                  gTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATc                                                                                                <  2:1213452/103‑1 (MQ=255)
                                                                                                  gTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATc                                                                                                >  1:1213452/1‑103 (MQ=255)
                                                                                                  gTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCggg                                                 <  1:1358714/150‑1 (MQ=255)
                                                                                                          ggATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGtt                                                            <  2:1197534/131‑1 (MQ=255)
                                                                                                          ggATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGtt                                                            >  1:1197534/1‑131 (MQ=255)
                                                                                                                aaCCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCAc                                                                                                     <  2:727971/84‑1 (MQ=255)
                                                                                                                aaCCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCAc                                                                                                     >  1:727971/1‑84 (MQ=255)
                                                                                                                aaCCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGc                                   >  1:1125935/1‑150 (MQ=255)
                                                                                                                            ccAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGa                       <  1:1046995/150‑1 (MQ=255)
                                                                                                                            ccAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGa                       <  1:297297/150‑1 (MQ=255)
                                                                                                                               ggTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAg                    <  2:1541277/150‑1 (MQ=255)
                                                                                                                                   ggAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCATCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCtt                <  1:20126/150‑1 (MQ=255)
                                                                                                                                     aGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGc              >  1:1208153/1‑150 (MQ=255)
                                                                                                                                      gTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGcc             <  1:1480210/150‑1 (MQ=255)
                                                                                                                                          cgGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAg         <  1:972793/150‑1 (MQ=255)
                                                                                                                                              tCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGc              <  1:934207/141‑1 (MQ=255)
                                                                                                                                              tCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGc              >  2:934207/1‑141 (MQ=255)
                                                                                                                                                aCCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAGGCCATc   <  2:319414/150‑1 (MQ=255)
                                                                                                                                                aCCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAGGCCATc   <  2:737513/150‑1 (MQ=255)
                                                                                                                                                 ccGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAGGCCATCg  >  1:930235/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACACCACCAACATCCCTGGGGCCGACAGTGGCGAACAGGTTGGCCTTGGTGCCCTCCGGGTTATGGAACAGCTCGCACTCCACCCCCAGCTCGGCCAGGTAATCACGGATAAAACCGATCAGCGCCAGGTTGGAGTCGCGGCTCACCGTGGCAAAGCCGATCAGCCGGGCCAGCAGCGCGCGACTGGCAAATTCACTCATCGCCCGGCACCCCATAGCTTGGCGCGGCGGTCGGGTTGAGGGCGCGGGTGACGTAATCCTGCATCTGCGGGCGATAGGCTTGCCACAGGCCATCG  >  NC_002947/4052284‑4052578

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: ‑