Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,898,660 T→C 10.0% intergenic (+66/‑105) PP_t16 → / → PP_t17 tRNA‑Asp/tRNA‑Val

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,898,6600TC10.0% 158.8 / 6.5 60intergenic (+66/‑105)PP_t16/PP_t17tRNA‑Asp/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base C (3/3);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCGGGAGAGCATCTGCC  >  NC_002947/1898513‑1898796
                                                                                                                                                   |                                                                                                                                         
ataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCgtttt                                                                                                                                          >  6:49823/1‑149 (MQ=21)
    cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCATCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAACCGACCTTAGAGTCGCTTTTTTtgt                                                                                                                                      <  3:191334/149‑1 (MQ=25)
    cgcAGCGGTAGTTCAGGCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAAAAAGCGACTCTAAGGTCGCTTTCTTtgt                                                                                                                                      <  7:226230/149‑1 (MQ=11)
      cAGCTGTAGTTCAGTCGGTTAGAATACCGGCCTGTCATGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGc                                                                                                                                    <  6:280458/149‑1 (MQ=255)
       aggggTAGTTCAGTCTGTTATAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTtgtg                                                                                                                                     <  8:258465/144‑1 (MQ=255)
       aGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTtgtgaa                                                                                                                                   >  7:258465/1‑147 (MQ=37)
                 cAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTgg                                                                                                                         >  1:124994/1‑149 (MQ=255)
                    tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTACAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGgtgt                                                                                                                      >  8:7948/1‑149 (MQ=255)
                      ggTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGttt                                                                                                                    <  2:299452/149‑1 (MQ=255)
                      ggTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGtt                                                                                                                     <  7:180790/148‑1 (MQ=255)
                       gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGttt                                                                                                                    >  5:165731/1‑148 (MQ=255)
                       gTTAGAATACCGGCCTGTCACGCCGCGGCTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGttt                                                                                                                    >  8:271764/1‑148 (MQ=255)
                        ttAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGATTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGa                                                                                                                  <  7:276111/149‑1 (MQ=255)
                            aaTACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCg                                                                                                               >  6:90548/1‑148 (MQ=255)
                              tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTtgt                                                                                                                                      >  3:222309/1‑123 (MQ=37)
                              tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTtgt                                                                                                                                      <  4:222309/123‑1 (MQ=37)
                              tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAg                                                                                                            >  8:183544/1‑149 (MQ=255)
                               cccGGCGTGACAGGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAgt                                                                                                           >  6:267699/2‑149 (MQ=255)
                                   gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGtgtg                                                                                                        >  2:226499/1‑148 (MQ=255)
                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGaaaa‑a                                                                                                  >  3:136485/1‑148 (MQ=255)
                                          aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACt                                                                                               >  2:1576/1‑149 (MQ=255)
                                              cGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGa                                                                                           <  5:49823/149‑1 (MQ=255)
                                              cGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGa                                                                                           <  4:168792/149‑1 (MQ=255)
                                                      gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTGTTTTTTTGGGTGATTTTTGGTGTTTGTCCGAAGTGTGAAAA‑AACTGGGATAAGTgga                                                                                   >  6:240769/1‑147 (MQ=18)
                                                       cgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCtt                                                                                  <  7:183544/149‑1 (MQ=255)
                                                          gTTCGAGTCCCGTCCGCTGCGGCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGcc                                                                               <  4:280758/149‑1 (MQ=255)
                                                                          cTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGcc                                                               >  8:87308/1‑149 (MQ=255)
                                                                              gcCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTa                                                           <  4:136485/149‑1 (MQ=255)
                                                                              gcCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTa                                                           <  5:90548/149‑1 (MQ=255)
                                                                                cATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGa                                                         <  1:1576/149‑1 (MQ=255)
                                                                                  tatTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGc                                                                                <  3:289463/124‑1 (MQ=255)
                                                                                   atTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCGTTTTTTGTGCCGGATTTTTGGGGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATg                                                      >  4:104155/1‑149 (MQ=255)
                                                                                      ccACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATgcgc                                                   >  1:12623/1‑149 (MQ=255)
                                                                                       cACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAACAACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATgcgc                                                   <  1:226499/149‑1 (MQ=255)
                                                                                            ccAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCcaca                                             >  5:90047/1‑149 (MQ=255)
                                                                                             cAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCacac                                            <  2:124994/149‑1 (MQ=255)
                                                                                              aGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACAcc                                           >  3:209548/1‑149 (MQ=255)
                                                                                               gggCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCg                                          <  2:12623/149‑1 (MQ=255)
                                                                                               gggCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCg                                          >  4:282729/1‑149 (MQ=255)
                                                                                                 gCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCga                                         <  7:271764/148‑1 (MQ=255)
                                                                                                  cccACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCgaga                                       <  4:209548/149‑1 (MQ=255)
                                                                                                         aaCACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACAt                                  <  7:87308/147‑1 (MQ=255)
                                                                                                         aaCACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATgg                                >  3:179879/1‑149 (MQ=255)
                                                                                                         aaCACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATgg                                >  7:252055/1‑149 (MQ=255)
                                                                                                         aaCACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATgg                                <  4:179879/149‑1 (MQ=255)
                                                                                                          acacCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATggg                               <  6:165731/149‑1 (MQ=255)
                                                                                                            acCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGGCGGATTTTTTGTGGTTGACCGAACTGTGAAAA‑AACTTCGATAAGTGGTTGCCAGAGCGCCAAGTAGCCCCTTGAATTCGCGCCACACCCAGAGCCATTGGTg                             >  5:17463/1‑149 (MQ=255)
                                                                                                             cctggAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGa                            <  1:260636/146‑1 (MQ=255)
                                                                                                             cccGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGa                            <  5:267699/149‑1 (MQ=255)
                                                                                                                    gTAACAAAGAAAGCGACCTTAGAGTCTCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTc                     <  7:7948/149‑1 (MQ=255)
                                                                                                                      aaCAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAg                   >  5:308860/1‑149 (MQ=255)
                                                                                                                       aCAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTAAGc                  >  7:96469/1‑149 (MQ=255)
                                                                                                                       aCAAAAAAAGCGACTCTAAGGTCGCTTTCTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGAt                           <  8:113703/140‑1 (MQ=39)
                                                                                                                       aCAAAAAAAGCGACTCTAAGGTCGCTTTCTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGAt                           >  7:113703/1‑140 (MQ=39)
                                                                                                                       aCAAAAAAAGCGACTCTAAGGTCGCTTTCTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAg                        <  4:55761/143‑1 (MQ=38)
                                                                                                                       aCAAAAAAAGCGACTCTAAGGTCGCTTTCTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAg                        >  3:55761/1‑143 (MQ=38)
                                                                                                                         aaagaaagCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCg                <  8:252055/149‑1 (MQ=255)
                                                                                                                         aaagaaagCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCg                >  1:284493/1‑149 (MQ=255)
                                                                                                                           aaaaagCGACTCTAAGGTCGCTTTCTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCgg               >  2:75128/3‑148 (MQ=38)
                                                                                                                                    ccTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCGGGAGAGCATCt     <  5:240769/149‑1 (MQ=255)
                                                                                                                                        agagTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCGGGAGAGCATCTGcc  <  3:282729/148‑1 (MQ=255)
                                                                                                                                                tttttttGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTa                         <  7:263427/117‑1 (MQ=255)
                                                                                                                                                tttttttGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTa                         >  8:263427/1‑117 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
ATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAACAAAGAAAGCGACCTTAGAGTCGCTTTTTTTGTGCCTGATTTTTGGTGTTTGACCGAAGTGTGAAAA‑AACTTCGATAAGTGCTTGCCAGAGCGCCAAGTCGCCCGTAGAATGCGCGCCACACCGAGAGACATGGGTGATTAGCTCAGCCGGGAGAGCATCTGCC  >  NC_002947/1898513‑1898796

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