Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,104,932 T→G 9.3% intergenic (+21/+132) 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,104,9320TG 9.3% 448.7 / 7.9 169intergenic (+21/+132)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (32/114);  new base G (13/7);  total (48/121)
Fisher's exact test for biased strand distribution p-value = 1.69e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTG  >  NC_002947/5104783‑5105079
                                                                                                                                                     |                                                                                                                                                   
   aTTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCg                                                                                                                                                  <  1:4227719/150‑1 (MQ=255)
    tttCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCgg                                                                                                                                                 >  1:4686300/1‑150 (MQ=255)
    tttCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCgg                                                                                                                                                 <  1:1621961/150‑1 (MQ=255)
    tttCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCCCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCgg                                                                                                                                                 <  2:5076830/150‑1 (MQ=255)
     ttCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCggg                                                                                                                                                <  2:2831521/150‑1 (MQ=255)
      tCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGc                                                                                                                                               <  2:1012406/150‑1 (MQ=255)
       ccATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGcc                                                                                                                                              <  1:5511422/150‑1 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                             >  1:6407658/1‑150 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                             <  2:4156191/150‑1 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGcgc                                                                                                                                             >  1:6407878/1‑148 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGccc                                                                                                                                             >  2:5612847/1‑150 (MQ=255)
         aTGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGCAGGGcccc                                                                                                                                            <  1:4309042/150‑1 (MQ=255)
         aTGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGCAGGGcccc                                                                                                                                            <  2:647322/150‑1 (MQ=255)
          tGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCt                                                                                                                                           <  2:3326336/150‑1 (MQ=255)
              cGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATcg                                                                                                                                       <  2:2064757/150‑1 (MQ=255)
                aaGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGGCCGGCGGGCCCCAATcgcc                                                                                                                                     >  1:5003337/1‑149 (MQ=255)
                 aGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGa                                                                                                                                    <  1:5966637/150‑1 (MQ=255)
                 aGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGa                                                                                                                                    <  2:1342408/150‑1 (MQ=255)
                 aGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGa                                                                                                                                    <  2:1342487/150‑1 (MQ=255)
                        aaTGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAGCGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAgg                                                                                                                             <  1:354586/150‑1 (MQ=255)
                        aaTGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAgg                                                                                                                             <  1:259631/150‑1 (MQ=255)
                          tgtgtgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGcc                                                                                                                           <  1:1861629/150‑1 (MQ=255)
                             gtgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCt                                                                                                                        <  2:3094137/150‑1 (MQ=255)
                             gtgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCt                                                                                                                        <  2:6980741/150‑1 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGc                                                                                                                                               >  2:1721454/1‑126 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGc                                                                                                                                               >  2:5492159/1‑126 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGcgc                                                                                                                       >  2:7190716/1‑148 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCcaca                                                                                                                  <  2:1983468/150‑1 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCcaca                                                                                                                  <  2:5922115/150‑1 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCcaca                                                                                                                  <  2:4584828/150‑1 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCcaca                                                                                                                  <  1:5218998/150‑1 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCcaca                                                                                                                  <  2:106999/150‑1 (MQ=255)
                                   aGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAACGCCGCTCCcaca                                                                                                                  <  2:3188136/150‑1 (MQ=255)
                                     tCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGAAGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGa                                                                                                                <  2:6884685/150‑1 (MQ=255)
                                      cGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGAt                                                                                                               <  2:5374910/150‑1 (MQ=255)
                                       gTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCca                                                                                                                    >  2:770552/1‑144 (MQ=255)
                                        tGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATaa                                                                                                             <  1:5140532/150‑1 (MQ=255)
                                         gCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAt                                                                                                            <  1:6700914/150‑1 (MQ=255)
                                         gCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAt                                                                                                            <  1:5353104/150‑1 (MQ=255)
                                          ccGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATg                                                                                                           <  2:2194955/150‑1 (MQ=255)
                                          ccGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATg                                                                                                           >  2:6184046/1‑150 (MQ=255)
                                            ggCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGCCCTGCGGGCCCCTATCGCGGCACAATGCCGCGGACACAGATaaatagg                                                                                                        >  1:3246355/1‑145 (MQ=255)
                                            ggCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGCCCGGCGGGCCCCTATCGCGACAAAAGGCCCCGCCCACAGGTaaatggc                                                                                                        >  1:3245667/1‑145 (MQ=255)
                                             gCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATgcgc                                                                                                        >  2:578079/1‑150 (MQ=255)
                                              cccGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCt                                                                                                       <  1:2339908/150‑1 (MQ=255)
                                               ccGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCCAGCTaaatcgcgc                                                                                                      >  1:573783/1‑143 (MQ=255)
                                                cGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCa                                                                                                     <  1:2134199/150‑1 (MQ=255)
                                                cGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCa                                                                                                     <  2:6244811/150‑1 (MQ=255)
                                                 gCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGCCACAAGGCCGCGCCCCCACAAAAAGCGCGCAg                                                                                                    >  2:2485270/1‑150 (MQ=255)
                                                  cATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGc                                                                                                   <  1:5876556/150‑1 (MQ=255)
                                                   aTGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGcc                                                                                                  <  1:2440155/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                               <  1:1088293/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                               <  1:7176724/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                               <  1:338756/150‑1 (MQ=255)
                                                       ccTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCagag                                                                                              >  2:2842667/1‑150 (MQ=255)
                                                       ccTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCagag                                                                                              >  2:2842867/1‑150 (MQ=255)
                                                           cGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGcaca                                                                                          <  1:3190054/150‑1 (MQ=255)
                                                            gAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAg                                                                                         >  2:6709307/1‑150 (MQ=255)
                                                            gAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGAGACACAAGGCCGCTCCCACAGATAATGAGCTCAGCCAGAGCACAg                                                                                         >  2:4847165/1‑150 (MQ=255)
                                                             aGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGc                                                                                        <  1:2780281/150‑1 (MQ=255)
                                                              ggACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCAGCGGGCCCCTAACGCGACACAAGGCCGCACCCACAGATAAAGCGCTAAGCCAGAGCACAGCg                                                                                       >  2:3551918/1‑150 (MQ=255)
                                                                aCTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTa                                                                                     <  1:4842170/150‑1 (MQ=255)
                                                                aCTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTa                                                                                     <  1:1044444/150‑1 (MQ=255)
                                                                 cTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTAt                                                                                    <  2:7135790/150‑1 (MQ=255)
                                                                   cccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTa                                                                                  <  1:1030018/150‑1 (MQ=255)
                                                                   cccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTa                                                                                  <  2:2081226/150‑1 (MQ=255)
                                                                     cAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTAcc                                                                                <  2:5336190/150‑1 (MQ=255)
                                                                        ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                             <  1:3551918/150‑1 (MQ=255)
                                                                        ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                             <  2:3917685/150‑1 (MQ=255)
                                                                         gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                            <  2:7211314/150‑1 (MQ=255)
                                                                         gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                            <  2:6821487/150‑1 (MQ=255)
                                                                          ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                           <  2:5480833/150‑1 (MQ=255)
                                                                          ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCAATCGCGACACAAGGCCGCACCCACAGAAAATGCGCGCAGCCAGAGCACAGCCTATTAccggggg                                                                           >  1:7419115/1‑145 (MQ=255)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGGTAATACGCTGTGCTCTGGCTGAGCGCATTATCTGTGGGa                                                                         <  1:7235373/150‑1 (MQ=21)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGa                                                                         <  1:2842867/150‑1 (MQ=255)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGa                                                                         <  1:2842667/150‑1 (MQ=255)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGa                                                                         <  2:2456/150‑1 (MQ=255)
                                                                              aaGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGc                                                                       <  2:3522981/150‑1 (MQ=255)
                                                                              aaGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGGCACAAGCCCCCACCCACACAAAAAGCGCACAACCCGAGCACAGCGGAATACCAGTGGGAGc                                                                       >  1:1135948/1‑150 (MQ=255)
                                                                               aGCACATTCTCGCAGCCCGCGTGATGGGGCACGGCTGCTCCACCGGCTGAGACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCg                                                                      <  2:5988099/150‑1 (MQ=255)
                                                                                       ctctAATACGGGGTGTTGGGGAAGGGCTGCTCCAAGGGCTGTTACTTTGTAATGGGGGGCCCTGGGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCCGCCTGAGCACCGCGTTTTACCTGTGGGAGCGGCCTtgtg                                                              <  1:5612847/150‑1 (MQ=37)
                                                                                           aaTACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGACTCGGGGGCCCTGCGGGACCCTATCGCGACACACGGCCGCTCCAACTGATAATGCGCTCAGCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTcgc                                                          <  1:6747936/150‑1 (MQ=255)
                                                                                           aaTACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTcgc                                                          <  2:332768/150‑1 (MQ=255)
                                                                                               ccGCGTGCTGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGaaa                                                      <  1:1090066/150‑1 (MQ=255)
                                                                                               ccGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGaaa                                                      <  1:6184046/150‑1 (MQ=255)
                                                                                                          ggAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGGATTACCTGTGGGAGCGGCCATGTGTCGCGAAAGGACTGCAACg                                           >  2:6729220/1‑150 (MQ=255)
                                                                                                            aaCGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACCCACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCa                                         <  2:3668704/150‑1 (MQ=255)
                                                                                                                  tgctCCACCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCa                                   <  2:4223375/150‑1 (MQ=255)
                                                                                                                  tgctCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGGCACAAGGCCGCTCCCACAGATAATGCGCTCAGCCCGAGCACAGCGGATTACCTGTGGGAGCGGCCGTGTGTCGCGAAAGGACTGCAACGCAGTCCCa                                   >  2:3538871/1‑150 (MQ=255)
                                                                                                                        aaCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAAGCGCTCAGCCAGCGCACCGCGTATTACCTGGGGGAGCGGCCGTGTGGCGCGGAAGGACTGCAACGCCGGCCCAGGTTCt                             >  2:6375402/1‑150 (MQ=255)
                                                                                                                        aaCGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGCCACAAGGCCGCTCCCACAGATAAGGCGCTCAGCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCt                             <  1:3692496/150‑1 (MQ=255)
                                                                                                                          cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                           <  2:3872783/150‑1 (MQ=255)
                                                                                                                            gcggATACCTTGGACTGGGGGGCCCGGGGGGCCCCTATCGCGACACACGGCCGCGCCCACAGATACGGCGCACCGCCAGCGCACAGCTTATTACGTGGGGGTGCGGCCTTGTGCCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGa                         <  1:820353/147‑1 (MQ=255)
                                                                                                                            gCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTACCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAGCGCAGTCCCAGGTTCTTTGa                         <  2:5753931/150‑1 (MQ=255)
                                                                                                                             cTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCGCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGATCTTTGAc                        >  2:6150001/1‑150 (MQ=255)
                                                                                                                             cTGATACCTGGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACCCAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGAc                        <  2:4799137/150‑1 (MQ=255)
                                                                                                                              tGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGAcc                       <  2:4454315/150‑1 (MQ=255)
                                                                                                                              tGATACCTTGGAATCGGGGGCCCATCGGGCCCCAATCGCCACACAAGGCCGCAACCACAAATAAAGCGCAAACCCAGAGCACAGCGTATTACCTGTGGGATCGGCCATGTTTCGCCCAAGGACAGCAAAGCAGTACCAGGTTAATTGaac                       >  2:3586296/1‑148 (MQ=255)
                                                                                                                                       ggattggggggCCGGGGGGGCCCCTTTGGCGCCACACGGCCGCTCCCACATATCTGGCGCTCCGCCTGCGCACGGGGTTTTCGGTGTGGGTGGGGCGTTGTGCCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTaa              <  1:2056129/144‑1 (MQ=255)
                                                                                                                                              ggggCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgccg       >  2:7441032/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTG  >  NC_002947/5104783‑5105079

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