Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 450,627 G→C 7.2% intergenic (+19/‑177) PP_0370 → / → PP_mr05 acyl‑CoA dehydrogenase family protein/RsmY ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947450,6270GC 7.2% 194.0 / 3.9 66intergenic (+19/‑177)PP_0370/PP_mr05acyl‑CoA dehydrogenase family protein/RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base G (37/24);  new base C (3/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

CTATTCGGGCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTCACCACTGTTC  >  NC_002947/450482‑450773
                                                                                                                                                 |                                                                                                                                                  
cTATTCGGGCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTgcg                                                                                                                                                  <  4:297380/148‑1 (MQ=255)
  aTTCGGGCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAg                                                                                                                                               >  6:191975/1‑149 (MQ=255)
     cGGGCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGACTTTTGGGGCTGCTGCGCAGcc                                                                                                                                             <  4:84505/148‑1 (MQ=255)
       ggCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCaa                                                                                                                                          >  8:64800/1‑149 (MQ=255)
       ggCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCGCAGCAGcccca                                                                                                                                          <  2:309161/149‑3 (MQ=255)
                ggTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACg                                                                                                                                  >  7:231732/1‑148 (MQ=255)
                        gcgGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGc                                                                                                                         >  1:28925/1‑149 (MQ=255)
                              cTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGATGCGCAGCCCAATCGCGACGCAAGGCCGCTCCca                                                                                                                    >  5:212720/1‑148 (MQ=255)
                                 cTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCcac                                                                                                                   <  8:231732/146‑1 (MQ=255)
                                      cgcgGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCg                                                                                                           <  7:292271/149‑1 (MQ=255)
                                       gcgGTTGCTGCCACGGGTGGATTCGCTGCTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGc                                                                                                          >  5:313860/1‑149 (MQ=255)
                                           ttGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTaa                                                                                                      <  5:64454/149‑1 (MQ=255)
                                               tgcCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGttt                                                                                                  >  8:287377/1‑149 (MQ=255)
                                               tgcCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGttt                                                                                                  >  2:56393/1‑149 (MQ=255)
                                                 ccACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTg                                                                                                 >  5:61091/1‑148 (MQ=255)
                                                       gTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGt                                                                                          >  1:341026/1‑149 (MQ=255)
                                                       gTGGATTCGCGGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGt                                                                                          >  7:3508/1‑149 (MQ=255)
                                                            ttCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAg                                                                                     >  2:60275/1‑149 (MQ=255)
                                                            ttCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCAACGGTATCGCGTAAGTTTGGATAgtgggagg                                                                                     >  5:251943/1‑142 (MQ=39)
                                                               gCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCgg                                                                                  >  1:260323/1‑149 (MQ=255)
                                                                cTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGc                                                                                 <  2:341026/149‑1 (MQ=255)
                                                                    tggCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTg                                                                                         >  5:357577/1‑137 (MQ=255)
                                                                    tggCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTg                                                                                         <  6:357577/137‑1 (MQ=255)
                                                                    tggCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCa                                                                             <  2:28925/149‑1 (MQ=255)
                                                                      gCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAg                                                                            <  6:313860/148‑1 (MQ=255)
                                                                      gCTTCGGTGAAGGCGGGGAGGGAGGCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGc                                                                           <  4:318467/149‑1 (MQ=255)
                                                                        ttCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGGTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGcc                                                                          >  6:21082/1‑148 (MQ=255)
                                                                        ttCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCa                                                                         >  1:13994/1‑149 (MQ=255)
                                                                            gTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCg                                                                      >  2:258977/1‑148 (MQ=255)
                                                                             tGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCa                                                                    <  6:61091/149‑1 (MQ=255)
                                                                             tGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTAAGCCATCGCa                                                                    >  7:359733/1‑149 (MQ=255)
                                                                                    ggggAGGGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCtgt                                                             <  8:267967/149‑1 (MQ=255)
                                                                                      ggAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCtgtg                                                            >  1:250544/1‑148 (MQ=255)
                                                                                       gagtgagtCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGtgaa                                                          >  6:203225/1‑149 (MQ=255)
                                                                                            agtCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGtgaatgaa                                                      >  7:151123/1‑148 (MQ=255)
                                                                                             gtCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGAGAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAAtc                                                    >  2:340749/1‑149 (MQ=255)
                                                                                                 ccttACCTGCCGGCTGCCGATCATTTCCGCTCGTTTTTTGGGGCCGCCGCGCATCCCAATCCCGACGCAATGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCt                                                <  6:251943/147‑1 (MQ=255)
                                                                                                  tttACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCt                                                >  6:242387/1‑148 (MQ=255)
                                                                                                    tACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCt                                                <  5:242387/146‑1 (MQ=255)
                                                                                                             gATGCCGAGCAGTTCTGGTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCtttt                                    <  5:21082/149‑1 (MQ=255)
                                                                                                              aTGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCttttt                                   <  2:13994/149‑1 (MQ=255)
                                                                                                                  cGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCt                                                <  1:5170/132‑1 (MQ=255)
                                                                                                                  cGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCt                                                >  2:5170/1‑132 (MQ=255)
                                                                                                                    aGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGACTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTaca                             <  8:151123/149‑1 (MQ=255)
                                                                                                                         ttCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACacgacg                        >  3:25845/1‑149 (MQ=255)
                                                                                                                            tGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtgg                     >  3:213982/1‑149 (MQ=255)
                                                                                                                             gATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtggt                    >  8:194788/1‑149 (MQ=255)
                                                                                                                             gATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtggt                    >  6:217803/1‑149 (MQ=255)
                                                                                                                              aTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtggt                    >  1:382096/1‑148 (MQ=255)
                                                                                                                              aTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtggt                    >  3:283527/1‑148 (MQ=255)
                                                                                                                              aTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtggt                    <  7:287377/148‑1 (MQ=255)
                                                                                                                              aTCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCAAATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACgtgg                     <  8:3508/147‑1 (MQ=255)
                                                                                                                               tCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGa                  >  5:332037/1‑149 (MQ=255)
                                                                                                                               tCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGa                  >  3:183971/1‑149 (MQ=255)
                                                                                                                               tCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGa                  >  7:194571/1‑149 (MQ=255)
                                                                                                                                 gATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGAc                 <  7:194788/148‑1 (MQ=255)
                                                                                                                                  aTTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGc                                                                                                          >  1:72758/1‑58 (MQ=255)
                                                                                                                                  aTTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGc                                                                                                          <  2:72758/58‑1 (MQ=255)
                                                                                                                                  aTTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACtt               <  5:203225/149‑1 (MQ=255)
                                                                                                                                     gattgggCTGCGCAGCAGCCCCATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGtgaa                                                          >  2:359718/5‑103 (MQ=38)
                                                                                                                                     gattgggCTGCGCAGCAGCCCCATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGtgaa                                                          <  1:359718/99‑1 (MQ=38)
                                                                                                                                     gattgggCTGCGCAGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACtttt             >  2:12727/5‑148 (MQ=255)
                                                                                                                                     gattgggCTGCGCAGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGACCGCGCGAATATGTGAAAGAAACTCATGTGTAAGCTTTTTAATACACGACATGGTGACtttt             >  1:323821/5‑148 (MQ=255)
                                                                                                                                       tGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTca           <  4:183971/148‑1 (MQ=255)
                                                                                                                                           gctgctGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTa                               <  6:62805/124‑1 (MQ=255)
                                                                                                                                           gctgctGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTa                               >  5:62805/1‑124 (MQ=255)
                                                                                                                                           gctgctGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTCACCACt      >  1:73506/1‑149 (MQ=255)
                                                                                                                                             tgctgcGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTcac          >  2:182684/1‑143 (MQ=255)
                                                                                                                                             tgctgcGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTcac          <  1:182684/143‑1 (MQ=255)
                                                                                                                                              gctgcGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTCACCACTGtt   >  7:356747/1‑149 (MQ=255)
                                                                                                                                               ctgcGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTCACAACTGTTc  >  4:289825/1‑149 (MQ=255)
                                                                                                                                                 gcgcAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGaa                                                                                      >  4:335544/1‑63 (MQ=255)
                                                                                                                                                 gcgcAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGaa                                                                                      <  3:335544/63‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTATTCGGGCAAGCTGGGTACGGCGCGGTTCTACTTTGCGCGGTTGCTGCCACGGGTGGATTCGCTGGTGGCTTCGGTGAAGGCGGGGAGTGAGTCGCTTTACCTGCTGGATGCCGAGCAGTTCTGATCGATTTTTGGGGCTGCTGCGCAGCCCAATCGCGACGCAAGGCCGCTCCCACAGGTATCGCCTAAGTTTGGATAGGTGGAAGCGGCCTCAGCCATCGCATATCTGTGAATGAATCTCCTGTGTAAGCTTTTTCCTACACGACGTGGTGACTTTTCACCACTGTTC  >  NC_002947/450482‑450773

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: