Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,641 G→C 7.5% noncoding (35/76 nt) PP_t59 ← tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,6410GC 7.5% 208.0 / 8.3 67noncoding (35/76 nt)PP_t59tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base C (3/2);  total (35/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 = 5.82e-01

GGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAG  >  NC_002947/4632499‑4632786
                                                                                                                                              |                                                                                                                                                 
gggCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                                                                                                                            <  1:1834556/150‑1 (MQ=255)
  gttgagTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGt                                                                                                                                          <  2:1077103/145‑1 (MQ=255)
  gCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCgg                                                                                                                                          <  2:92945/150‑2 (MQ=255)
     gggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGct                                                                                                                                       <  2:1145181/150‑1 (MQ=255)
         aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAc                                                                                                                                   <  2:780228/150‑1 (MQ=255)
         aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAc                                                                                                                                   <  2:749116/150‑1 (MQ=255)
          ttGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAcc                                                                                                                                  >  1:825496/1‑150 (MQ=255)
              ggCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACt                                                                                                                              <  1:1277295/150‑1 (MQ=255)
               gCAGAATGCATTCGAACGAAAAAAGGCCCCGCAGGGCGTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCGTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTg                                                                                                                             >  1:155445/1‑150 (MQ=255)
               gCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTg                                                                                                                             >  2:1786174/1‑150 (MQ=255)
                cAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  2:825496/150‑1 (MQ=255)
                  gAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGATCGGGAAACGGGGCTCGAATTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGc                                                                                                                          <  1:1008792/150‑1 (MQ=255)
                        aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAg                                                                                                                           >  1:1035041/1‑143 (MQ=255)
                        aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAg                                                                                                                           <  2:1035041/143‑1 (MQ=255)
                         ttCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTccc                                                                                                                   >  1:1338739/1‑150 (MQ=255)
                              acgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTGTTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggggtc                                                                                                                                      <  2:641564/126‑6 (MQ=255)
                              acgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTGTTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggggtc                                                                                                                                      >  1:641564/1‑121 (MQ=255)
                              acgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtc                                                                                                              <  2:1338739/150‑1 (MQ=255)
                                gaaAAAAGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtcgt                                                                                                            >  2:863568/1‑150 (MQ=255)
                                gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                                                                                                                <  2:604159/130‑1 (MQ=255)
                                gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                                                                                                                >  1:604159/1‑130 (MQ=255)
                                   aaaaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGAt                                                                                                         >  1:554352/1‑150 (MQ=255)
                                     aaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg                                                                                                       <  2:554352/150‑1 (MQ=255)
                                     aaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg                                                                                                       >  1:646140/1‑150 (MQ=255)
                                     aaGGCCCCGAAGGGCCTTTTTTCAACGAATTCCAGCATGCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg                                                                                                       <  2:646140/150‑1 (MQ=21)
                                                gggCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTa                                                                                            <  2:555463/150‑1 (MQ=255)
                                                   ccTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCg                                                                                         >  2:1047179/1‑150 (MQ=255)
                                                    cTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGa                                                                                        <  1:45460/150‑1 (MQ=255)
                                                        tttCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGATCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTc                                                                                    >  1:1414574/1‑150 (MQ=255)
                                                                  ttACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCg                                                                          >  1:1938467/1‑150 (MQ=255)
                                                                  ttACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCg                                                                          >  1:698662/1‑150 (MQ=255)
                                                                       gACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAg                                                                                                                           <  2:68246/96‑1 (MQ=255)
                                                                       gACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAg                                                                                                                           >  1:68246/1‑96 (MQ=255)
                                                                        aCTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAAcc                                                                    >  2:1771104/1‑150 (MQ=255)
                                                                          ttCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTAGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAAcccc                                                                  >  2:758709/1‑150 (MQ=255)
                                                                               cTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGa                                                              >  1:1862203/1‑149 (MQ=255)
                                                                               cTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGa                                                              <  2:1862203/149‑1 (MQ=255)
                                                                                   aaCCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCaa                                                         >  1:1526332/1‑150 (MQ=255)
                                                                                      cTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCaaaaa                                                      >  2:88183/1‑150 (MQ=255)
                                                                                       tGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTaaa                                                                                >  1:1294010/1‑123 (MQ=255)
                                                                                       tGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTaaa                                                                                <  2:1294010/123‑1 (MQ=255)
                                                                                         tAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg                                                                                                       <  1:623546/98‑1 (MQ=255)
                                                                                         tAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg                                                                                                       >  2:623546/1‑98 (MQ=255)
                                                                                            aTCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACtttt                                                >  1:291282/1‑150 (MQ=255)
                                                                                             tCTTTATTTGTAGCGGGAAACGAGTCTCGAACTCCGGGCCCGGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTa                                               <  2:155445/150‑1 (MQ=37)
                                                                                              cTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAAc                                                    >  2:129916/1‑144 (MQ=255)
                                                                                              cTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAAc                                                    <  1:129916/144‑1 (MQ=255)
                                                                                               tttATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCggg                                                                                                                                         >  2:221087/1‑56 (MQ=25)
                                                                                               tttATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCggg                                                                                                                                         <  1:221087/58‑3 (MQ=25)
                                                                                               tttATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTAtt                                             >  1:532007/1‑150 (MQ=255)
                                                                                                         gCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGa                                   <  2:532007/150‑1 (MQ=255)
                                                                                                           gggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGc                                 <  2:1414574/150‑1 (MQ=255)
                                                                                                                 cGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTcgc                           <  1:863568/150‑1 (MQ=255)
                                                                                                                   agaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCGAAAAACTTTTATTTTTTTTCCGAGCCTTCgcgc                         <  2:1916098/150‑1 (MQ=255)
                                                                                                                   agaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCgcgc                         <  2:1598798/150‑1 (MQ=255)
                                                                                                                       tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGa                                   >  2:5962/1‑136 (MQ=255)
                                                                                                                       tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGa                                   <  1:5962/136‑1 (MQ=255)
                                                                                                                       tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGt                     >  1:1170628/1‑150 (MQ=255)
                                                                                                                             tCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggcc               >  1:1004180/1‑150 (MQ=255)
                                                                                                                             tCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggcc               >  1:880206/1‑150 (MQ=255)
                                                                                                                                  agcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcg          >  2:153544/5‑150 (MQ=255)
                                                                                                                                   gcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcgg         >  1:1404226/4‑150 (MQ=255)
                                                                                                                                    cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc        <  2:698662/150‑1 (MQ=255)
                                                                                                                                     ccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggcc       >  1:1799382/1‑150 (MQ=255)
                                                                                                                                      cGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggcca      >  1:539821/1‑150 (MQ=255)
                                                                                                                                          cTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCcagcag  >  2:1837341/1‑150 (MQ=255)
                                                                                                                                          cTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCcagcag  <  1:1047179/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAG  >  NC_002947/4632499‑4632786

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