Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,350,908 T→A 6.4% intergenic (+1/‑157) PP_2945 → / → arfB two‑component system sensor histidine kinase/response regulator/peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,350,9080TA 6.4% 182.4 / 2.7 63intergenic (+1/‑157)PP_2945/arfBtwo‑component system sensor histidine kinase/response regulator/peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (17/40);  new base A (2/2);  total (20/43)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGA  >  NC_002947/3350761‑3351051
                                                                                                                                                   |                                                                                                                                               
cccGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATg                                                                                                                                                >  8:30402/1‑149 (MQ=255)
         cAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCAcgcg                                                                                                                                       <  6:378168/149‑1 (MQ=255)
            ggTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTc                                                                                                                                    >  1:208067/1‑149 (MQ=255)
            ggTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTc                                                                                                                                    <  5:96089/149‑1 (MQ=255)
                 cGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGACGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCg                                                                                                                               >  5:386684/1‑149 (MQ=255)
                   ggCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCgct                                                                                                                             <  5:249291/149‑1 (MQ=255)
                   ggCTAGGCCTGTCACTTTCATATGGCATTGTAAAGAATCACGTTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCTTTCTACCgct                                                                                                                             <  6:386684/149‑1 (MQ=255)
                        ggCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCgctggct                                                                                                                         <  5:392971/148‑1 (MQ=255)
                         gCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGctggctgg                                                                                                                       >  4:67860/1‑149 (MQ=255)
                         gCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGctggctgg                                                                                                                       >  6:397876/1‑149 (MQ=255)
                          ccTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCTCGCGTTCTACCGCtggctggc                                                                                                                      <  2:189863/149‑1 (MQ=255)
                              tCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCAt                                                                                                                    <  7:340915/147‑1 (MQ=255)
                                ccTTTCCTATGGCATTGTAAAGAAACACGGTGCGGAGATTACCGTATCAAGCCAGGTCTGAGTGGGAACAAGGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCGGTGATGGCGCACGCGTTCTACCGGTGGCTGGCATGaa                                                                                                                 <  6:363164/148‑1 (MQ=255)
                                 cTTTCCTATGGCATTGTAAAGAAACACGGTGGGTAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGt                                                                                                               <  3:13466/149‑1 (MQ=255)
                                      cTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGAGAGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACg                                                                                                          >  4:188122/1‑149 (MQ=255)
                                         tGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACgg                                                                                                         <  8:137631/147‑1 (MQ=255)
                                         tGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGcc                                                                                                       <  1:99507/149‑1 (MQ=255)
                                              ttGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGa                                                                                                   <  3:41182/148‑1 (MQ=255)
                                              ttGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGAc                                                                                                  <  2:208067/149‑1 (MQ=255)
                                                 tAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGa                                                                                                <  7:406923/148‑1 (MQ=255)
                                                   aaGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGACGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTg                                                                                             >  4:215867/1‑149 (MQ=255)
                                                    aGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTg                                                                                             >  7:240477/1‑148 (MQ=255)
                                                     gAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGa                                                                                           <  5:330451/149‑1 (MQ=255)
                                                        acacGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCAAGGCCCTGACGAGTGGAggg                                                                                        >  4:383495/1‑149 (MQ=255)
                                                         cacGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGATCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGc                                                                                       <  7:314764/149‑1 (MQ=255)
                                                                gggAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCtt                                                                                 <  2:132756/148‑1 (MQ=255)
                                                                gggAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTg                                                                                >  2:302145/1‑149 (MQ=255)
                                                                   agaTAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGt                                                                             >  3:418852/1‑149 (MQ=255)
                                                                   agaTAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGt                                                                             <  4:32812/149‑1 (MQ=255)
                                                                       aaCCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGc                                                                         >  8:386003/1‑149 (MQ=255)
                                                                        aCCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCt                                                                        >  8:360207/1‑149 (MQ=255)
                                                                        aCCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCt                                                                        >  4:88366/1‑149 (MQ=255)
                                                                         ccGTATCAAGCCAGGGCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCtt                                                                       <  8:292652/149‑1 (MQ=255)
                                                                          cGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTg                                                                      <  7:386003/149‑1 (MQ=255)
                                                                          cGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTg                                                                      <  8:240477/149‑1 (MQ=255)
                                                                            tATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTgcg                                                                    <  1:352737/149‑1 (MQ=255)
                                                                               cAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                 >  2:188756/1‑149 (MQ=255)
                                                                                   ccAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGAt                                                             <  7:360207/149‑1 (MQ=255)
                                                                                    cAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGAtt                                                            <  3:88366/149‑1 (MQ=255)
                                                                                     aGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTc                                                           <  3:383495/149‑1 (MQ=255)
                                                                                      ggTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCa                                                          <  2:143656/149‑1 (MQ=255)
                                                                                       gTCGGAGTGGGAACAACGTGTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaa                                                         <  5:317121/149‑1 (MQ=255)
                                                                                        tCGGAGTGGGAACTACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTAGCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaaa                                                        <  3:288700/149‑1 (MQ=255)
                                                                                        tCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaaa                                                        <  3:67860/149‑1 (MQ=255)
                                                                                            aGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAAc                                                    <  4:418852/149‑1 (MQ=255)
                                                                                                ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACTCTGTGAAGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                 >  3:336981/1‑132 (MQ=255)
                                                                                                ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGAAGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                 <  4:331142/132‑1 (MQ=255)
                                                                                                ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGAAGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                 >  3:331142/1‑132 (MQ=255)
                                                                                                ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGAAGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                 <  4:336981/132‑1 (MQ=255)
                                                                                                       cGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATtgtg                                         >  3:119188/1‑149 (MQ=255)
                                                                                                       cGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATtgtg                                         <  7:30402/149‑1 (MQ=255)
                                                                                                            cGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTAt                                     <  1:188756/148‑1 (MQ=255)
                                                                                                             gtgtACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATg                                    <  8:32784/148‑1 (MQ=255)
                                                                                                                tACAACTTCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGGTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                 <  5:128394/148‑1 (MQ=255)
                                                                                                                tACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                 <  5:122222/148‑1 (MQ=255)
                                                                                                                   aaCTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGc                             >  2:293190/1‑149 (MQ=255)
                                                                                                                        ccGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGca                         <  4:413055/148‑1 (MQ=255)
                                                                                                                        ccGGTAAAGCAGAGCGCAACGCTGTGACGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGc                          <  3:215867/147‑1 (MQ=255)
                                                                                                                              aaGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCtg                  <  3:202052/149‑1 (MQ=255)
                                                                                                                                gCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCtgt                 >  4:425179/1‑148 (MQ=255)
                                                                                                                                 cAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGa               <  2:295910/149‑1 (MQ=255)
                                                                                                                                 cAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGa               <  5:397876/149‑1 (MQ=255)
                                                                                                                                                tgATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGa  <  1:302145/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGA  >  NC_002947/3350761‑3351051

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