Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,790 T→A 5.7% D258V (GAC→GTC)  thiI ← tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,750,7900TA 5.7% 201.0 / 2.6 70D258V (GAC→GTC) thiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base A (2/2);  total (36/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 = 8.41e-01

GTCGATGATCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGTTGAAG  >  NC_002947/5750642‑5750933
                                                                                                                                                    |                                                                                                                                               
gTCGATGATCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCtgttgt                                                                                                                                                 >  1:192320/1‑149 (MQ=255)
   gatcatCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGa                                                                                                                                              >  3:2756/1‑149 (MQ=255)
    atgatCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCGTGAGGGTCACGCCCATGTAGCTGTTGTCGAc                                                                                                                                             >  8:162589/1‑149 (MQ=255)
       atCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCtt                                                                                                                                          <  1:325091/149‑1 (MQ=255)
          gACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgtt                                                                                                                                       <  1:273231/149‑1 (MQ=255)
            cAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgtt                                                                                                                                       >  6:271677/1‑147 (MQ=255)
              ggTTCGGCAGGGTCTGTCTGGACACCTGGGAACTCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAgg                                                                                                                                   <  2:319601/149‑1 (MQ=255)
              ggTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAgg                                                                                                                                   <  2:192320/149‑1 (MQ=255)
                       gggTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGcc                                                                                                                           >  7:204921/1‑148 (MQ=255)
                          tCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGAGCTCGCCAAcc                                                                                                                       >  4:250920/1‑149 (MQ=255)
                              gctggACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGGCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCActt                                                                                                                   <  1:295653/149‑1 (MQ=255)
                                tggACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAAGCGGTCGGCCATATGGGCGGCGCCGCGCAGCAACATGCGCATGAGGGTCAAGCCCATGTAGCTGTTGTCCACCATGTTGAGGATCACGCCAACCAcataa                                                                                                                 >  1:275618/1‑145 (MQ=255)
                                  gACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAAGCActtctt                                                                                                                <  4:2756/148‑1 (MQ=255)
                                  gACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACttcttc                                                                                                               <  3:282183/149‑1 (MQ=255)
                                       cTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgt                                                                                                                                        <  1:351816/119‑1 (MQ=255)
                                       cTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgt                                                                                                                                        >  2:351816/1‑119 (MQ=255)
                                       cTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCACCCACTTCTTCGAACg                                                                                                          >  3:301692/1‑149 (MQ=255)
                                       cTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACg                                                                                                          >  7:304944/1‑149 (MQ=255)
                                           gAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGTAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCAc                                                                                                      <  6:173699/149‑1 (MQ=255)
                                               tCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCt                                                                                                   >  6:316028/1‑148 (MQ=255)
                                               tCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCt                                                                                                   >  6:371941/1‑148 (MQ=255)
                                                  cTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCtgatg                                                                                               <  8:304944/149‑1 (MQ=255)
                                                      gCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACa                                                                                           <  3:82664/149‑1 (MQ=255)
                                                      gCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACa                                                                                           >  3:60616/1‑149 (MQ=255)
                                                       ccGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGACGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAg                                                                                          >  3:125644/1‑149 (MQ=255)
                                                       ccGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGACGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAg                                                                                          >  3:131040/1‑149 (MQ=255)
                                                        cGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAg                                                                                          >  2:356426/1‑148 (MQ=255)
                                                        cGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGt                                                                                         <  7:81243/149‑1 (MQ=255)
                                                         ggTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTa                                                                                        >  2:288305/1‑149 (MQ=255)
                                                         ggTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTa                                                                                        >  2:110164/1‑149 (MQ=255)
                                                              ccAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgtt                                                                                                                                       <  6:21536/97‑1 (MQ=255)
                                                              ccAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCttgtt                                                                                                                                       >  5:21536/1‑97 (MQ=255)
                                                              ccAGCGCGTCAAACTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGAGCTTGAGGGTCACGCCCAAGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCAAACCACTTCTTCGAACGGCACGCTGATGAAAAGTACgcgc                                                                                   >  7:135967/1‑149 (MQ=255)
                                                               cAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCgct                                                                                  >  6:382688/1‑149 (MQ=255)
                                                                  cgcgTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCgctggc                                                                               <  5:51800/149‑1 (MQ=255)
                                                                     gTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGgctgc                                                                            >  1:202053/1‑149 (MQ=255)
                                                                          tCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGGTTGAGGGTCACGCCCATGTTGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCg                                                                       >  7:277959/1‑149 (MQ=255)
                                                                          tCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCAGGCTGCTGCCg                                                                       >  4:319231/1‑149 (MQ=255)
                                                                                   ggTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTc                                                              <  5:371941/149‑1 (MQ=255)
                                                                                   ggTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTc                                                              <  8:40948/149‑1 (MQ=255)
                                                                                   ggTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTc                                                              <  1:356426/149‑1 (MQ=255)
                                                                                       ggCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCca                                                            >  8:150113/1‑147 (MQ=255)
                                                                                           atatGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGt                                                       >  8:138288/1‑148 (MQ=255)
                                                                                            tatGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTaa                                                     >  2:370994/1‑149 (MQ=255)
                                                                                                ggcggcGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATgg                                                  <  4:60616/148‑1 (MQ=255)
                                                                                                ggcggcGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATgg                                                  <  5:316028/148‑1 (MQ=255)
                                                                                                ggcggcGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATgg                                                  <  2:275618/148‑1 (MQ=255)
                                                                                                 gcggcgCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCCGGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGc                                                >  5:124524/1‑149 (MQ=255)
                                                                                                 gcggcgCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATggca                                                >  8:308889/1‑147 (MQ=255)
                                                                                                    gcgcCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCAc                                             <  5:135867/149‑1 (MQ=255)
                                                                                                    gcgcCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCAc                                             >  2:169241/1‑149 (MQ=255)
                                                                                                      gccgcGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACtt                                           <  1:288305/149‑1 (MQ=255)
                                                                                                        cgcgCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTcc                                         <  5:271677/149‑1 (MQ=255)
                                                                                                         gcgcAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCa                                        <  4:301692/149‑1 (MQ=255)
                                                                                                              gcaTCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACg                                   <  2:37303/149‑1 (MQ=255)
                                                                                                                   atGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAg                              >  4:120924/1‑149 (MQ=255)
                                                                                                                   atGCGCTTGAGGGTCACGCACATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAg                              >  2:230456/1‑149 (MQ=255)
                                                                                                                           gAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTgg                       <  7:162589/148‑1 (MQ=255)
                                                                                                                           gAGGGTCACGCCCATGTAGCTGTTGACGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTgg                       <  4:131040/148‑1 (MQ=255)
                                                                                                                           gAGGGTCACGCCCATGTAGCTGTTGACGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTgg                       <  4:125644/148‑1 (MQ=255)
                                                                                                                             gggTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCa                    <  2:183289/149‑1 (MQ=255)
                                                                                                                              ggTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCAc                   >  6:175448/1‑149 (MQ=255)
                                                                                                                              ggTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCAc                   <  3:319231/149‑1 (MQ=255)
                                                                                                                                  aCGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAGCCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTGCACAGGTAATGGGCCACTGCGATTACGCCCAGCTCGTGGGCACGgcg               >  4:129109/1‑148 (MQ=255)
                                                                                                                                    gCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGgccgc             >  4:159698/1‑149 (MQ=255)
                                                                                                                                      ccATGTAGCTGTTGTTGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCAGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACACCCAGCTCGTGGGCACGgccgcc            >  1:119154/1‑148 (MQ=255)
                                                                                                                                      ccATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGgccgcc            >  1:320979/1‑148 (MQ=255)
                                                                                                                                        aTGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAg         <  3:120924/149‑1 (MQ=255)
                                                                                                                                         tGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAgg        <  5:28837/149‑1 (MQ=255)
                                                                                                                                         tGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAgg        <  1:370994/149‑1 (MQ=255)
                                                                                                                                           tAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGtt      <  5:382688/149‑1 (MQ=255)
                                                                                                                                               tgttgtCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGTTg     >  6:89895/1‑146 (MQ=255)
                                                                                                                                                 ttgtCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGTTgaag  >  8:163449/1‑147 (MQ=255)
                                                                                                                                                 ttgtCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGTTgaag  <  1:230456/147‑1 (MQ=255)
                                                                                                                                                   gtCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCa                    >  7:333417/1‑127 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTCGATGATCGACAGGTTCGGCAGGGTCTGGCTGGACACCTGGGAAATCGCTTCGCCGGTCACCAGCGCGTCAATCTGCAGGCGGTCGGCCATATGGGCGGCGCCGCGCAGCATCATGCGCTTGAGGGTCACGCCCATGTAGCTGTTGTCGACCTTGTTGAGGATCTCGCCAACCACTTCTTCGAACGGCACGCTGATGAACAGTACGCGCTGGCTGCTGCCGTATTTTTTCCACAGGTAATGGGCCACTTCCATTACGCCCAGCTCGTGGGCACGGCCGCCGAGGTTGAAG  >  NC_002947/5750642‑5750933

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