Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,101,957 A→C 6.7% intergenic (+73/+80) PP_4488 → / ← phhR hypothetical protein/sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,101,9570AC 6.7% 495.8 / 4.6 174intergenic (+73/+80)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (65/91);  new base C (12/3);  total (77/97)
Fisher's exact test for biased strand distribution p-value = 5.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGT  >  NC_002947/5101809‑5102106
                                                                                                                                                    |                                                                                                                                                     
gACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCag                                                                                                                                                      <  1:2550342/150‑1 (MQ=255)
 aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                     <  2:5813652/150‑1 (MQ=255)
 aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                     <  1:2896718/150‑1 (MQ=255)
  ccTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCagag                                                                                                                                                    <  1:663543/150‑1 (MQ=255)
       gAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAAGCGCGCAGCCAGAGCACAg                                                                                                                                               >  2:1121090/1‑150 (MQ=255)
        aGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGc                                                                                                                                              <  1:8032122/150‑1 (MQ=255)
          gACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCGATCGCGACACAAGGCCGCTACCACACATAATTCGCTCAGCCCGAGCACAGCGt                                                                                                                                            >  2:1130739/1‑150 (MQ=255)
               ccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTAc                                                                                                                                       <  1:5013117/150‑1 (MQ=255)
               ccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCCGGCCCCAAAACCCAAAAAAAGCCCCACCCACACAAAATACCCACACCCACAGCACAACGTaaaac                                                                                                                                       >  1:4469164/1‑146 (MQ=255)
                   ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                                                                                   <  1:4104535/150‑1 (MQ=255)
                   ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                                                                                   <  1:6152939/150‑1 (MQ=255)
                    gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                                                                                  <  2:3712365/150‑1 (MQ=255)
                    gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                                                                                  <  1:6558101/150‑1 (MQ=255)
                    gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                                                                                  <  2:3598138/150‑1 (MQ=255)
                    gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                                                                                  <  2:408284/150‑1 (MQ=255)
                     ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCACCCACAGATAATGCGCTCAGCCCCCGCACAGCGTATTACCAGTgg                                                                                                                                 >  2:5032951/1‑150 (MQ=255)
                     ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCACCCACAGATAAAGCGCTAAGCCAGAGCACACCGTATTACCGGTgg                                                                                                                                 >  2:1854688/1‑150 (MQ=255)
                        gAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAg                                                                                                                              <  2:4631605/150‑1 (MQ=255)
                         taGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGc                                                                                                                             <  1:9139455/149‑1 (MQ=255)
                         aaGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGc                                                                                                                             <  2:1430806/150‑1 (MQ=255)
                          aGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCg                                                                                                                            <  1:1359422/150‑1 (MQ=255)
                           gCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCAGCGGGCCCCGATCGCGACCCAACGCCCCGGCCACCGCTAATGCGCCCACCCCCAGCACCGCCGATTAGCGGGGGGGGCgg                                                                                                                           >  1:4248873/1‑150 (MQ=255)
                              caTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGCATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCCGAGCACAGCGTATTACCTGTGGGCGCGGccg                                                                                                                        >  2:3527658/1‑149 (MQ=255)
                               aTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCtt                                                                                                                       <  1:9140698/150‑1 (MQ=255)
                                ttCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCCCTCCCACAGAAAAAGAGCTCAGCCAGAGCACACCGTATTAACTGTGGGAGCGGCCAtg                                                                                                                      >  2:7191111/1‑150 (MQ=255)
                                 tctcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCAGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTtgt                                                                                                                     >  2:7051089/1‑150 (MQ=255)
                                   tcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTgtgt                                                                                                                   <  2:7864127/150‑1 (MQ=255)
                                          ccGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGaaa                                                                                                            <  2:6545354/150‑1 (MQ=255)
                                           cgcgTGTTGGGTATCGGCTGCTCCAACGGCTGATACGTTGGATTCGGGGGCCCTGCGGGCCCCTATCGCGACACACTGCCGATCCCACAGATATTGCGCTCAGCCTGAGCACAGAGTTTTACCTGTGGGATGGGCCTTTTGTCGCGAAAg                                                                                                           <  2:5589643/150‑1 (MQ=255)
                                           cgcgTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAg                                                                                                           <  2:8233803/150‑1 (MQ=255)
                                           cgcgTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAg                                                                                                           <  2:2103824/150‑1 (MQ=255)
                                             cgTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGa                                                                                                         <  2:756717/150‑1 (MQ=255)
                                                       aaCGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCa                                                                                               <  2:3751109/150‑1 (MQ=255)
                                                         ccccgattCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGGCACAAGGCCGCGCCCCCAGAAAATGCGCGCAGCCAGCGCACACCGGATTAACGGTGGGGGCGGGCGTGTGGGGCGGAAGGACGGCCACGCCGt                                                                                             >  2:1745261/8‑150 (MQ=255)
                                                            ctgctCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATACCGACACAAGGCCGCGCCCACAGATAATTCGCGCACCCCGAGCCCCCCCTATTACCGGTGGGAGCGGCCGGGTGTCGCGGAACGACTGCAACGCACTccc                                                                                          >  2:6145881/1‑150 (MQ=255)
                                                               ctCCACCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcctttc                                                                                       <  2:2230014/150‑5 (MQ=38)
                                                               ctCCAACGGCTGATACCTGGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACCCCCGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTCCCTGGGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcctttc                                                                                       <  2:1959044/150‑5 (MQ=255)
                                                                tCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGcaca                                                                                                                                                >  2:7448171/1‑92 (MQ=255)
                                                                tCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCGGTTGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGt                                                                                      >  2:2874696/1‑150 (MQ=255)
                                                                 ccAACGGCTGATACCTTGGAATCGGGGGGCCGGCGGGCCCCGAACGCGACCCCAGGCCGCGCCCCCACAAAAATCGCGCACCCCCCGCACAGCGGAGTAACGGGGGGGGCGGCCGTGGGTGGCGGAAGGGCGGCAAAGCAGACCCACGtt                                                                                     >  2:4254003/1‑150 (MQ=255)
                                                                 ccAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCGATCGCGCCACAAGGCCGCTCCCACACATAATTCGCTCAGCCCGAGCACAGCCTATTAACGGTGGGGGGGGGCGTGTGTGGCGCAAGGACGGCCACCCCGTCCCACGtt                                                                                     >  2:4254277/1‑150 (MQ=255)
                                                                  cAACGGCTGATACCTGGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAAGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCCGTCCCAGGtgc                                                                                    >  2:7456354/1‑148 (MQ=255)
                                                                   aaCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCt                                                                                   <  2:7517913/150‑1 (MQ=255)
                                                                   aaCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATAGCGACACAACGCCGCTCCCACAGATAAAGCGCGCAGCCAGAGCACAGCGGAATACCAGTGGGAGCGGCCTTGTGTCGCGAAAGGACGGCAACCCCGACCCAGGTTcg                                                                                   >  2:8572802/1‑149 (MQ=255)
                                                                           aTACCTTGGAATCGGGGGCCCGGCGGGCCCCAATCGCGACACAAGGCCCCGCCCACACAAAAATCGCGCACCCAGAGCACACCGGATTAACAGTGGGAGCGGCCGTGTGTCGCGAAAAGACAGCAACCCACTCCCAGGGTCTTTGAAcaa                                                                           >  2:1671055/1‑150 (MQ=255)
                                                                            tACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaac                                                                          <  2:2754379/150‑1 (MQ=255)
                                                                                ttGGAATCGGGGGCCCGGCGGGCCCCTAACGCGACACAAGGCCGCTCCCACAGATAATGCGCGCAGCCAGAGCCCAGCGTATTACCGGTGGGAGCGGCCGGGTGTCGCGAAAGGACTGCAACGCAGTCCCACGTTCTTTGACCaacaaca                                                                      >  2:4557179/1‑149 (MQ=255)
                                                                                    aaTCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCAAAGAAAAATCGCTCAGCCAGAGCACAGCGTATTAACTGTTGGAGCGGCCTTGTGTAGCGAAAGGACAGCAACGCAGTCCCAGGTTCTTTGACCAAAAACaaaac                                                                  >  2:5489043/1‑148 (MQ=255)
                                                                                        gggggCCCGGCGGGCCCCTAGCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCCGAGCACCGCTTATTACGTGTGGGCGGGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgcc                                                              <  1:1738133/150‑1 (MQ=255)
                                                                                           ggCCCTGCGGGCCCCTATCGCGACACAAGGCCCCGCCCCCACATAATGCGCGCAGCCAGAGCACAGCCTATTACCGGTGGGAGCGGCCGTGTGTCGCGGAAGGACGGCAACGCAGTCCCAGGGTCTTTGTCCAAAAACTAATccgccgcc                                                           >  1:2785971/1‑148 (MQ=255)
                                                                                              ccTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCGTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATACGCCGACTTg                                                        >  2:441187/1‑150 (MQ=255)
                                                                                               cTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGt                                                       <  2:9095126/150‑1 (MQ=255)
                                                                                                ggCGGGCCCCTTTCGCGACACAAGGCCGCTCCCACAGATACTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGtt                                                      <  2:6055068/149‑1 (MQ=255)
                                                                                                    ggCCCCTTTCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                  <  2:8606043/150‑1 (MQ=255)
                                                                                                    ggCCCCTATCGCGACACAAGGGCGCTCCCACAGATACGGCGCTCAGCCAGTGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGCAAGGACTGCAACGCAGTCCCTGGTTCTTTGACCATCAACGAATCCGCCGACATGTTAAGc                                                  <  2:746458/150‑1 (MQ=38)
                                                                                                    ggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                  <  1:5264939/150‑1 (MQ=255)
                                                                                                    ggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                  <  2:5055937/150‑1 (MQ=255)
                                                                                                    ggCCCCTATCGCGACACAAGGCCGCACCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCAGTGGGAGCGGCCGTGTGTCGCGAAATGAATGCAACGCAGTCCCAGGTTCTTTTACCAACAACTAAACCGCCGACATGTTAAGc                                                  >  1:2088447/1‑150 (MQ=255)
                                                                                                     gCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAATCCCTCAGCCAGAGCACAGCGTATTAACTGTGGGAGCGGCCATGTGTCGCGAAAGGACAGCAACCCAGTCCCAGGTTCTTTGTCCAACAACTAATACGCCGACTTGATAAGCg                                                 >  1:6705441/1‑150 (MQ=38)
                                                                                                      ccccTATCGCGCCACAAGGCCGCTCCCACCGATAATGCGCTCCGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTGGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACACCTAATCCGCCGACTTGTTAAGCGa                                                <  2:7744778/150‑1 (MQ=255)
                                                                                                      ccccTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  2:5641328/150‑1 (MQ=255)
                                                                                                      ccccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCATAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  1:8253399/150‑1 (MQ=255)
                                                                                                      ccccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  2:7313201/150‑1 (MQ=255)
                                                                                                      ccccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  2:5072049/150‑1 (MQ=255)
                                                                                                      ccccTATCGCGAACCAAGGCGGCTACCACTGATACTGCGGTCAGCCAGAGCAGAGCGTATTACCTGTGGGAGCGGCGTTGTGTAGCGACAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  2:521850/150‑1 (MQ=255)
                                                                                                      ccccTATCGCCCCACACGGCCGCTCCCACAGTTACTGCGCTCCGCCAGAGCACCGCTTATTACGTGTGGGCGCGGCCTTGTGCCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  1:1671055/150‑1 (MQ=255)
                                                                                                       cccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAg                                               >  2:8273374/1‑150 (MQ=255)
                                                                                                            tCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                <  1:8320950/144‑1 (MQ=255)
                                                                                                            tCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                                >  2:8320950/1‑144 (MQ=255)
                                                                                                             cgcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATcc                                         <  2:8123542/150‑1 (MQ=255)
                                                                                                               ccacacAACCCCGCTCCCACACATAATGCCCACACCCACACCACACCGAATTACCTGTGGGAGCGGCCTTGTGACGCGAAAGGACACCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCt                                       <  2:3824618/148‑1 (MQ=255)
                                                                                                                  cacaAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCaa                                    >  2:7952199/1‑150 (MQ=255)
                                                                                                                     aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCt                                                                                                                                     <  2:1105438/50‑1 (MQ=255)
                                                                                                                     aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCt                                                                                                                                     >  1:1105438/1‑50 (MQ=255)
                                                                                                                       ggCCCCCCCCACATCTAAGGCGCTCAGCCAGCGCAAAGCGTATTACCTGTGGGCGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgt                               <  1:1812488/150‑1 (MQ=255)
                                                                                                                        gCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTAGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgtt                              <  2:5559519/150‑1 (MQ=255)
                                                                                                                         ccGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgttg                             <  1:8885750/150‑1 (MQ=255)
                                                                                                                          cGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCTTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTgg                            <  1:5477983/150‑1 (MQ=255)
                                                                                                                                cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                      >  2:8030035/1‑150 (MQ=255)
                                                                                                                                 acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                     >  2:3216591/1‑150 (MQ=255)
                                                                                                                                 acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                     >  2:899482/1‑150 (MQ=255)
                                                                                                                                 acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                     >  1:8345159/1‑150 (MQ=255)
                                                                                                                                    gacaaTGCGCTCAGCCAGAGCACAGCGTTATACCTCGGGGAGCTGCCATGTGACGCGAAAGGACTGCAACGCTGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTTAGCGAGTAATCCCTCAACTTGTGGGCAATTGTCgt                  <  2:2780453/147‑1 (MQ=255)
                                                                                                                                    gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                  >  1:1017409/1‑150 (MQ=255)
                                                                                                                                    gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                  <  2:6351914/150‑1 (MQ=255)
                                                                                                                                    gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                  <  1:7535038/150‑1 (MQ=255)
                                                                                                                                         tGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCg             >  1:5851044/1‑150 (MQ=255)
                                                                                                                                            gcTCAGCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGaaa          <  1:3502399/150‑1 (MQ=255)
                                                                                                                                            gcTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCGACTTGTTGGCAATTGTCGTGTGCGaaa          >  1:6968083/1‑150 (MQ=255)
                                                                                                                                            gcTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGaaa          <  2:8126286/150‑1 (MQ=255)
                                                                                                                                              tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAaca        >  2:241882/1‑150 (MQ=255)
                                                                                                                                              tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAaca        >  2:2112403/1‑150 (MQ=255)
                                                                                                                                                cgCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc      <  1:7294916/149‑1 (MQ=255)
                                                                                                                                                aGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc      >  1:1487672/1‑150 (MQ=255)
                                                                                                                                                aGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc      >  1:1868567/1‑150 (MQ=255)
                                                                                                                                                 gCCAGAGCACAGCGTATTACCTGTGGGTGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAccc     <  1:7045561/150‑1 (MQ=255)
                                                                                                                                                  ccAGAGCACAGCTTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa    <  1:4057895/150‑1 (MQ=255)
                                                                                                                                                  ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa    <  2:3700995/150‑1 (MQ=255)
                                                                                                                                                    agagCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGt  >  1:3671314/1‑150 (MQ=255)
                                                                                                                                                    agagCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGt  >  1:2605039/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGT  >  NC_002947/5101809‑5102106

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