Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,439 T→G 6.8% intergenic (+118/+156) PP_0635 → / ← PP_0636 group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947743,4390TG 6.8% 216.6 / 6.5 74intergenic (+118/+156)PP_0635/PP_0636group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base G (2/3);  total (36/38)
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.97e-01

ATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGG  >  NC_002947/743295‑743582
                                                                                                                                                          |                                                                                                                                               
aTTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTc                                                                                                                                               <  5:99610/147‑1 (MQ=255)
aTTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTc                                                                                                                                               >  6:99610/1‑147 (MQ=255)
aTTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCgg                                                                                                                                             <  7:340294/149‑1 (MQ=255)
   aaCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGtt                                                                                                                                           <  8:251229/148‑1 (MQ=255)
     ccGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCg                                                                                                                                        >  6:47063/1‑149 (MQ=255)
      cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGt                                                                                                                                       >  8:27639/1‑149 (MQ=255)
        gCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTa                                                                                                                                     >  4:80594/1‑149 (MQ=255)
            aGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACa                                                                                                                                   <  5:282773/147‑1 (MQ=255)
             gCCGCGTAACCTGAACCGCCGTGTACGGAACCGCACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCGAAGGTTGCTTCACATa                                                                                                                                 <  1:47865/148‑1 (MQ=18)
                   tAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGAGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGgca                                                                                                                          <  6:293340/149‑1 (MQ=18)
                       cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCTTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                        >  7:286472/1‑147 (MQ=255)
                       cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                        >  7:314962/1‑147 (MQ=255)
                         gAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGc                                                                                                                    >  8:8175/1‑149 (MQ=255)
                          aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGTGGACGGCGGGTGTGAATCCCACCCCCTCCCCCCCCTCTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGa                                                                                                                            <  1:296334/149‑1 (MQ=18)
                          aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcg                                                                                                                   >  1:307493/1‑149 (MQ=255)
                              gccgTGTACGTTAACGTACGGACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTAGATTTCGCGAAATCAGAGGTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               <  6:218801/149‑1 (MQ=37)
                              gccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               <  7:8175/149‑1 (MQ=255)
                              gccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               <  7:27639/149‑1 (MQ=255)
                                tctttACGGAACCTTACGCACTGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTTCGCGAAATCAGACCTGCCTTGTGTGCCCGGTATTTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               <  1:165140/143‑1 (MQ=37)
                                    tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTg                                                                                                         >  7:18598/1‑149 (MQ=255)
                                    tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACAATCGGGTGCGTTACATAGCGGACCCACAGCGTCGTGAcatc                                                                                                         >  2:165140/1‑146 (MQ=37)
                                      cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGaca                                                                                                                          <  2:76739/130‑1 (MQ=255)
                                      cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGaca                                                                                                                          >  1:76739/1‑130 (MQ=255)
                                       ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               >  7:322465/1‑140 (MQ=255)
                                       ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                                               <  8:322465/140‑1 (MQ=255)
                                       ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg                                                                                                      <  4:276004/149‑1 (MQ=255)
                                         aaCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCggg                                                                                                    <  1:40462/149‑1 (MQ=255)
                                            cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGc                                                                                                                                         <  1:331606/109‑1 (MQ=255)
                                            cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGc                                                                                                                                         >  2:331606/1‑109 (MQ=255)
                                                       tggtgTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTc                                                                                      <  8:213620/149‑1 (MQ=255)
                                                       tggtgTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTc                                                                                      >  1:141848/1‑149 (MQ=255)
                                                         gtgtGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgc                                                                                    >  7:320859/1‑149 (MQ=255)
                                                          tgtgAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgcg                                                                                   <  8:234260/149‑1 (MQ=255)
                                                           gtgAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCgcgg                                                                                  >  4:247485/1‑149 (MQ=255)
                                                                     ggcggGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg                                                                         >  3:270018/1‑148 (MQ=255)
                                                                      gcggGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCa                                                                       <  3:59962/149‑1 (MQ=255)
                                                                           tgtgAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCTTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGAc                                                                  <  8:286472/149‑1 (MQ=255)
                                                                               aaT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTg                                                                >  4:233358/1‑147 (MQ=255)
                                                                               aaT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTg                                                                >  2:116965/1‑147 (MQ=255)
                                                                                aT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGAAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg                                                             <  8:314962/149‑1 (MQ=255)
                                                                                       cGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAgggg                                                          <  8:18598/149‑1 (MQ=255)
                                                                                       cGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                           <  1:116965/148‑1 (MQ=255)
                                                                                       cGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                           <  3:233358/148‑1 (MQ=255)
                                                                                        gCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTGTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                           >  1:235473/1‑148 (MQ=255)
                                                                                                     cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                                  <  5:47063/148‑1 (MQ=255)
                                                                                                      tCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                                  >  7:265411/1‑147 (MQ=255)
                                                                                                        gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGt                                              >  3:318417/1‑149 (MQ=255)
                                                                                                         aTTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTg                                                                >  5:82703/1‑130 (MQ=255)
                                                                                                         aTTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTg                                                                <  6:82703/130‑1 (MQ=255)
                                                                                                            gcgcGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGc                                           <  1:55198/148‑1 (MQ=255)
                                                                                                              gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAgg                                                            <  8:46635/129‑1 (MQ=255)
                                                                                                              gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAgg                                                            >  7:46635/1‑129 (MQ=255)
                                                                                                              gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGt                                        >  6:100675/1‑149 (MQ=255)
                                                                                                               cgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGtt                                       >  2:294938/1‑149 (MQ=255)
                                                                                                                gAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAgaga                                                                   >  5:268568/1‑120 (MQ=255)
                                                                                                                gAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAgaga                                                                   <  6:268568/120‑1 (MQ=255)
                                                                                                                   aTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTc                                     >  6:340409/1‑147 (MQ=255)
                                                                                                                   aTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                   <  7:114179/149‑1 (MQ=255)
                                                                                                                      aGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagc                                 >  5:149085/1‑148 (MQ=255)
                                                                                                                              cTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                        <  5:340409/149‑1 (MQ=255)
                                                                                                                                    gtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCa                   >  6:181649/1‑149 (MQ=255)
                                                                                                                                     tgtccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCAc                  >  4:291634/4‑149 (MQ=37)
                                                                                                                                     tgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCa                   <  2:235473/148‑1 (MQ=255)
                                                                                                                                        ccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                                          >  5:309775/1‑74 (MQ=255)
                                                                                                                                        ccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                                          <  6:309775/74‑1 (MQ=255)
                                                                                                                                        ccGCTATGTAACGCAACCGAAGGTTGCTTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                        <  1:271566/140‑1 (MQ=38)
                                                                                                                                        ccGCTATGTAACGCAACCGAAGGTTGCTTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGAAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                        >  2:271566/1‑140 (MQ=38)
                                                                                                                                               gtgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCg         >  2:52097/1‑148 (MQ=255)
                                                                                                                                               gtgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCg         <  1:52097/148‑1 (MQ=255)
                                                                                                                                                tgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa        <  6:149085/148‑1 (MQ=255)
                                                                                                                                                       aCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa        <  3:138101/141‑1 (MQ=255)
                                                                                                                                                       aCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa        >  4:138101/1‑141 (MQ=255)
                                                                                                                                                       aCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAgg  <  3:276834/147‑1 (MQ=255)
                                                                                                                                                       aCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAgg  >  4:276834/1‑147 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
ATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAAT‑‑‑‑CCGCCT‑‑‑‑‑CCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGG  >  NC_002947/743295‑743582

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