Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,127,094 T→G 12.2% intergenic (+98/‑84) PP_4511 → / → PP_4512 AraC family transcriptional regulator/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,127,0940TG12.2% 479.8 / 3.4 180intergenic (+98/‑84)PP_4511/PP_4512AraC family transcriptional regulator/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (60/82);  new base G (21/2);  total (96/84)
Fisher's exact test for biased strand distribution p-value = 7.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

GGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTC  >  NC_002947/5126945‑5127243
                                                                                                                                                     |                                                                                                                                                     
ggCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCt                                                                                                                                                       <  1:3423403/150‑1 (MQ=255)
    ttCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAt                                                                                                                                                   <  2:4808460/150‑1 (MQ=255)
       gggTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                                                                <  2:3871198/150‑1 (MQ=255)
       gggTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                                                                <  1:1839584/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  1:1815258/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  2:3553645/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  2:5178762/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  2:6531742/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  2:2595522/150‑1 (MQ=255)
         gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCAGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                                              <  2:4374682/150‑1 (MQ=255)
          tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTtgtg                                                                                                                                             <  1:1654313/150‑1 (MQ=255)
                ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAc                                                                                                                                       <  1:394551/150‑1 (MQ=255)
                  tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCt                                                                                                                                     <  2:7357582/150‑1 (MQ=255)
                   ttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTg                                                                                                                                    <  2:2154770/150‑1 (MQ=255)
                   ttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTg                                                                                                                                    <  2:3351105/150‑1 (MQ=255)
                   gaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTg                                                                                                                                    >  1:113832‑M1/84‑150 (MQ=255)
                   gaagaacagcgcgccaatgtgcagggcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTg                                                                                                                                    >  1:8282114‑M1/84‑150 (MQ=255)
                    tGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                                                   <  1:9135027/150‑1 (MQ=255)
                     ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                                                  <  1:7170362/150‑1 (MQ=255)
                      gCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAg                                                                                                                                 <  2:5006563/150‑1 (MQ=255)
                      gCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAg                                                                                                                                 <  1:9115002/150‑1 (MQ=255)
                        tttGCGGCGAGAGTCTGGCGTTAACTGACAGCAGTCTTGCGTCAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGa                                                                                                                               <  1:6028106/150‑1 (MQ=255)
                        tttGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCGCCTACAGGATGCCCCCGCCGGAATGGTGGGCGAAAACAGTAggt                                                                                                                               >  1:8176820/1‑149 (MQ=255)
                        tttGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCACCAACACGATGCCCCCCCCGGAATGTTGTGCAAAAACGGTAggc                                                                                                                               >  2:1719682/1‑149 (MQ=255)
                          tGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                                                             <  1:3160504/150‑1 (MQ=255)
                          agcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                                                             >  2:4186924‑M1/77‑150 (MQ=255)
                            cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCgg                                                                                                                           <  2:3710202/150‑1 (MQ=255)
                            cggcgAGAGTCTGGCGATAACCGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCgg                                                                                                                           <  1:7773706/150‑1 (MQ=255)
                             ggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCggg                                                                                                                          <  2:4754747/150‑1 (MQ=255)
                             ggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCggg                                                                                                                          <  2:1497865/150‑1 (MQ=255)
                             ggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCTACAGGATTCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCggg                                                                                                                          >  2:7094474/1‑150 (MQ=255)
                              gcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGGAAACCCCGCCTACACGCGGCCCCCCCCGGGATGTTGTGCAAAAACAGTGGGAGCgggg                                                                                                                         >  1:4740837/1‑149 (MQ=255)
                                 agagTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTa                                                                                                                      <  2:5140005/150‑1 (MQ=255)
                                   atgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTAcc                                                                                                                    >  2:827418‑M1/68‑150 (MQ=255)
                                     tCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCcg                                                                                                                  <  1:4454702/150‑1 (MQ=255)
                                     tCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCTACACGATGCCCCCGCCTGAAAGTTGTGCCAAACCGGGAGGAGCCGGTTTAcccc                                                                                                                  >  2:4912558/1‑149 (MQ=255)
                                      gtGGCGATATCTGACAGCAGTCTTGCGACAGATTGTGTTTGCCTGAGGCTTTTCACTTTGCCTTTTCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATTCCCCCGCCTGTATGTTGTGCTAAACCTGTAGGAGCGGGTTTACCcgc                                                                                                                 <  2:5315653/149‑1 (MQ=38)
                                      cTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCcgc                                                                                                                 <  2:407188/150‑1 (MQ=255)
                                          cGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                             <  2:5906708/150‑1 (MQ=255)
                                           gATAACTGCCAGCAGTCTTGCGACTGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTCCCGGATGCCCCCGCCTGTATGTTGTGCAAAACCTGTAGGAGGGGGTTTACCCGCGaagg                                                                                                            <  2:7486456/150‑3 (MQ=255)
                                               ccTGACACCAGTCTTGCGACTGACTGTGTTCGCCTGAGCCTTTACACTTTCCCTTCGCCGCCCCCTGCGCGGGTAAACCCCCTCCTACAGTATGCCCCCCCCTTAATGTTGTGCAAAACCTGTAGTAGCGGGTTTACCCGCGAAAGGGcc                                                                                                        <  2:8863587/149‑1 (MQ=25)
                                               aCTGACAGCCGTCTTGCGACTGATTGTGTTTGCCTGAGGCTTTACACTTTGCCTTGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTACAGGATGCCCCCGCCTGAATTTTGTGCAAAACCTGAAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                                                        <  2:673693/150‑1 (MQ=38)
                                                cTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCGGAATGGTGGGCAAAACCTGGAGGCGCGGGGTTGCCCGCGAAAGGGCCg                                                                                                       >  2:5761357/1‑150 (MQ=39)
                                                cTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTACTACAGGATTCCCCCGCCGGAATTTTGTGCACAACCTGTTGGAGCGGGTTTTCCCGCGCAAGGGccc                                                                                                       >  2:1672406/1‑149 (MQ=255)
                                                    ccgccggcTTGCGACCGACGGTGTGGGCCGGCGGCTTTACACTTTGCCTTGGCCGGCCCCTGCGCGGGTACACCCGCTCCCACCGGAGGCCCCCGCCTGCCTGTTGGGCCAACCCTGTCGGCGGGGGTTTACCCGCGAAAGGGCCGGCAc                                                                                                   <  2:2072893/143‑1 (MQ=255)
                                                          tCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAGCCGCTGGTACAGGATGCGCCCGCCTGAATGTTGTGCTAAAGCTGTTGGAGCGGGTTTAAGCGCGTATGGGGCGGGACTGCCTg                                                                                             >  2:9143400/1‑150 (MQ=255)
                                                          tCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCGACAGGCTGCCCCCGCCGGGATGTTGTGCCAAACCGGTAGGGGCGGGGTTACCCGCGCAAGGGCCGGCACGGCCTg                                                                                             >  2:6536093/1‑150 (MQ=255)
                                                            ttGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAACACGATGCCCCCCCCAGAAAGTTTTGCAAAAACAGGAGGAGCGGGGTTAACCCCCAAAGGGCCGGCACTGCCGGTa                                                                                           >  1:8215912/1‑150 (MQ=255)
                                                                gACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAAGATGCCCCCGCCAGAATGTTTTGCAAAAACAGTAGGAGCGGGGTTACCCCCCCAAGGGCCGGCACAGCCGGgagaga                                                                                       >  1:8355724/1‑150 (MQ=255)
                                                                   aGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGATGCCCCCCCCGGAATGTTGTGCAAAACCGGTAGGGGCGGGTTTACCCGCGCAAGGGCCGGCCCTGCCCGTAGAGAtct                                                                                    >  2:2613761/1‑150 (MQ=255)
                                                                       tgtgTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGATGCCCCCGCCAGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCGGTAGAGATCTCCAt                                                                                >  2:7317355/1‑150 (MQ=255)
                                                                         tgtTTGCCTGTGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGTATGCCCCCGCCTGATTGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCa                                                                              <  2:2032847/150‑1 (MQ=255)
                                                                           tttGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACACCCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCGAAACCGGTGGGAGCGGGTTTAACCCCGAAAGGGCCGGCACTGCCGGTAGGGATCGACCTCAAg                                                                            >  2:6996818/1‑150 (MQ=255)
                                                                              gCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCGCCAACACGATTCCCCCCCCAGAAAGTTGTGCAAAAACAGTAGGAGCGGGTTTACCCCCCAAAGGGCCGGCCCAGCCGGGAGAGAAAACCATCAAGGAc                                                                         >  1:9017595/1‑150 (MQ=255)
                                                                                 tGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGACAGGCTGCCCCCCCCGGAATGTTGTGCCAAACCGGGAGGAGCGGGGTTACCCGCCAAAGGGCCCGCACAGCCGGGAGAGATCTCCAGCCACGAcccc                                                                      >  2:2182829/1‑149 (MQ=255)
                                                                                       tttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTACTAAAGGATGCCCCCCCCGGAATTTTGTTCAAAAACCGTAGGAGAGGGTTTAACCCCGAAAGGGCCCGCACAGACGGTAGATAAAGCCAAAAAAGACCCACCAAGa                                                                >  2:2514993/1‑150 (MQ=255)
                                                                                        ttACACTCTGCCTTTGCCGGGCCCTGCGCGGGTACACCCGCGCCTCCAGTATGGCCCCTCCTTATTGTTGTGAAAAACGTTTAGGTGGGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGaa                                                               <  2:6102042/150‑1 (MQ=255)
                                                                                          acacTTTGCCTTTGCCGGCCCCTTCGCGGGGAAAACCCCGCCTACACGAATCCCCCCCCCGCAATTTGTGCAAAAACAGGAGGGGCGGGGTTACCCCCGCAAGGGGCCGCCCAGCCGGTAGAGGGATACCAGAAGGACCCACCAAGaaac                                                             >  2:312588/1‑148 (MQ=255)
                                                                                           tctacaggcagTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                                                   >  1:7575495‑M1/12‑79 (MQ=255)
                                                                                           tctacaggcagTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                                                   <  2:7575495‑M1/68‑1 (MQ=255)
                                                                                                           GCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCAAAACCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCAGTAGAGATCTCCCTCAAGGACCCACCATGAACCACCGCACCCCCCTCGgc                                            <  2:5284764‑M1/150‑1 (MQ=255)
                                                                                                             gtaggagcgGGTTTACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGgcgc                                          >  1:2416042/7‑150 (MQ=255)
                                                                                                               cTTCGCGGGTAAACCCGCTCCTCCAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc                                        <  2:6018263/150‑1 (MQ=255)
                                                                                                               cTTCGCGGGTAAACCCGCTCCTACCGGATGCCCCCGCCTGACTGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc                                        <  2:3493541/150‑1 (MQ=255)
                                                                                                                   gcgGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGca                                    >  1:7847258/1‑150 (MQ=255)
                                                                                                                      ggTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGcacct                                 >  2:7995852/1‑148 (MQ=255)
                                                                                                                       gTAAACCCGCTCCTACAGGATGCCCCCGCCTTAATGTTGTGCAAAACCTGTAGGTGGGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAtt                                <  2:4587084/150‑1 (MQ=255)
                                                                                                                       ccaaaCCCGCGCCTCCCGGAGCCCCCCCCCTCAATGTGGCGCAAACCCGTCAGGCGGGGTTTTCCCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAtt                                <  1:2723797/148‑1 (MQ=255)
                                                                                                                         aaaCCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCACACCGCCCTCGGCGCACTGCACATTgg                              >  1:4147680/1‑150 (MQ=255)
                                                                                                                          aaCCCGCTCCTCATGGAGGCCCCCGCTTTTTTTTTGGGAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgc                             <  1:3645632/150‑1 (MQ=37)
                                                                                                                            cccGCTCCTCCAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgc                           <  2:6248756/150‑1 (MQ=255)
                                                                                                                             ccGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcg                          >  1:2832431/1‑150 (MQ=255)
                                                                                                                             ccGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcg                          >  2:3806256/1‑150 (MQ=255)
                                                                                                                             ccGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcg                          >  2:2933745/1‑150 (MQ=255)
                                                                                                                              cGCTACTTCAGTATGGCGCCCGCTGTATGTTGCGCTATGCTTGTATGTTCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcgc                         <  2:6841887/150‑1 (MQ=255)
                                                                                                                               gCTCCTACAGGAGGCCCCCGCCGGAATGTTGTGCAAACCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCt                        <  2:1278144/150‑1 (MQ=255)
                                                                                                                                 tCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGt                      >  1:6209835/1‑150 (MQ=255)
                                                                                                                                 tCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGt                      >  1:4796254/1‑150 (MQ=255)
                                                                                                                                    tACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                   >  1:2583753/1‑150 (MQ=255)
                                                                                                                                    tACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                   >  1:931715/1‑150 (MQ=255)
                                                                                                                                    tACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                   >  1:1441198/1‑150 (MQ=255)
                                                                                                                                      ccgtaTGCCCCCTCCTGATTGTTGGGCAAACCCTTACGGAGGGGGTTTCCCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                 <  1:2954127/146‑1 (MQ=255)
                                                                                                                                      ccggAGGCCCCCGCATTATTGTTGTCCAAACCGTGTAGGAGGGGTTTAACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                 <  1:212526/148‑1 (MQ=255)
                                                                                                                                      cAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                 >  1:1151437/1‑150 (MQ=255)
                                                                                                                                      cAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                 >  1:5878964/1‑150 (MQ=255)
                                                                                                                                        ggATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCgg               >  1:5061320/1‑150 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGTTTGTTGTGCAAAACCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              <  2:2854118/150‑1 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              >  2:6522633/1‑150 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              >  2:4297620/1‑150 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              >  2:3512117/1‑150 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              >  2:8124809/1‑150 (MQ=255)
                                                                                                                                         gATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc              >  2:8423340/1‑150 (MQ=255)
                                                                                                                                             cccccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGCCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTa          <  2:3370074/150‑1 (MQ=255)
                                                                                                                                              ccccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAc         >  1:7363107/1‑150 (MQ=255)
                                                                                                                                               cccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAcc        >  2:2213381/1‑150 (MQ=255)
                                                                                                                                                    cTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGt   >  1:3912055/1‑150 (MQ=255)
                                                                                                                                                    cTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGt   >  1:3912280/1‑150 (MQ=255)
                                                                                                                                                    cTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGt   >  1:6910219/1‑150 (MQ=255)
                                                                                                                                                    cTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGt   >  2:3170463/1‑150 (MQ=255)
                                                                                                                                                     tGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTc  >  2:5019008/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTC  >  NC_002947/5126945‑5127243

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: