Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,609 A→T 5.4% intergenic (+117/‑91) pyk → / → PP_4302 pyruvate kinase/urea transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,6090AT 5.4% 376.9 / 7.9 112intergenic (+117/‑91)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base A (50/56);  new base T (3/3);  total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 116 total aligned reads displayed.

CTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4895461‑4895757
                                                                                                                                                     |                                                                                                                                                    
cTAACAGGCTGCGGATCGAGACTTTGATCTGAGATGGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCAGGCGGGCACAAGGGAATGCAGTCCCTGTGGGCGCGGCCTTGTGTCGCGATTGGAGGGCaaa                                                                                                                                                     <  2:2578735/150‑1 (MQ=255)
cTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCaaa                                                                                                                                                     <  2:4419253/150‑1 (MQ=255)
         tGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCCCTCCCACAGGGCCACGCTGGCACACGGGAATGCAGGCCCTGGGGGAGCGGCCTTGTGGCGCGATTGGAGGGCAAAGccctccct                                                                                                                                            <  2:4605904/150‑1 (MQ=255)
          gCGGATCGAGACTTGGATCTGAGTTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCCCCCACCGGGCCACGCTGGCACACGGGAATGCAGACCCTGTGGGAGCGGCCTTGTGTCGCGATGGGAGGGCTTTGccctccaat                                                                                                                                           <  1:4080154/150‑4 (MQ=255)
             gATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCATGTGTCGCGATTGGAGGGCAAAGCCCTACCAGAAt                                                                                                                                        >  2:1626298/1‑150 (MQ=255)
              aTCGAGACTTTGATCTGAGATGGCCGGGGGCGCGTTGCGCCCCTTTCGCCACACACGCCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGACGCGATTGGAGGGCAAAGCCCTCCCTGAATa                                                                                                                                       <  2:2080104/150‑1 (MQ=255)
               tCGAGACTTTTATCTGAGATTGCCGGGGGCGGGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACTGGGCCACGCTGGCACATGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAc                                                                                                                                      <  1:4078330/150‑1 (MQ=39)
               tCGAGACTTTGTTCTGAGTTTGCCGGGGGCGCGTTGCGCCCCTTGCGCGACACACGGCCGCGCCCACAGGGCCACGCGGGCACACGGGAAGGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAc                                                                                                                                      <  1:3518209/150‑1 (MQ=38)
                   gaCTTTTATCTTAGTTTGCGGGGGGCGCTTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACATGGTACTGCAGTCCCTGTGGGCGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGct                                                                                                                                  <  1:782367/150‑1 (MQ=255)
                   gaCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGct                                                                                                                                  <  2:4465847/150‑1 (MQ=255)
                      tttGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTTGAGGGCAAAGCCCTCCCTGAATACCGCTCTa                                                                                                                               >  1:4485724/1‑150 (MQ=255)
                       tttaTCTGAGATTGCCGGGGGCGCTTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACATGGGATTGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTAc                                                                                                                              <  1:1329358/147‑1 (MQ=39)
                           tCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCgc                                                                                                                          >  2:4060596/1‑149 (MQ=255)
                           tCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCGGTGGGAGCGGCCATGTGTCGCGATTTGAGGGCAAAGCCCACCCAGAAAAACGCGCTAcccca                                                                                                                          >  1:4458079/1‑148 (MQ=38)
                           tCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCGGGCACAAGGGAATGCAGTCCCAGTGGGAGCGGCCTTGTGTCGCGATTTGAGGGCAAAGCCCTCCCTGAATACCGCTCTAACCGa                                                                                                                          >  2:3327034/1‑150 (MQ=255)
                           tCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCACCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCACCCTGAAAACCCCGCTACCCGa                                                                                                                          >  1:1799619/1‑150 (MQ=38)
                                   ggcgggggggGCGTTGCGCCCCTTGCGCCACACCCGGCCGCCCCCCCTGGGCCACGCTGGAACAAGGGAATGCAGGCCGTGGGGGCGGGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGaa                                                                                                                  <  1:2958551/146‑1 (MQ=255)
                                     ccGGGGGCGCGTTGCGCCCCTTTCGCGGCACAAACCCCCACCCACAGGGCCACCCGCGCACAAGGGAATGCACTGCCAGAGGGGGCGGGCATGTTTAGCGTTTGGTGGGGAAAAACCACCCTGAATAACGCACAACCCCACCCCGGcaaa                                                                                                                >  1:1152671/1‑150 (MQ=255)
                                         gggCGCTTTGCGCCCCTTTCGCTACACAAGGCCGCTCCCACAGGGCCACGTGGGCCCACGGGAATGCAGTCCGTGTGGGCGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  1:4511314/150‑1 (MQ=255)
                                         gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  2:3316713/150‑1 (MQ=255)
                                          ggCGCGTTGCGCCCCTTTCGCGACCCCCGGCCTCACCCACGGGGCCCCGCGGCCACATGGTATGCACGCCAGTGTGGGAGCGGCTTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCTCGGAAAACCGAg                                                                                                           <  1:2862130/150‑1 (MQ=255)
                                          ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACCGGGCCACGCTGGCACAAGGGATTGCAGTCCCTGTGGGCGCGGCCTTGTGACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  1:2184759/150‑1 (MQ=255)
                                           gcgcGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                                                          <  2:5213842/150‑1 (MQ=255)
                                             gcgTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGcc                                                                                                        <  2:3023268/150‑1 (MQ=255)
                                                 tGCGCCCCTTTCGCTACACACGGCCGCCCCCCCGGGCCCAGGGGGGCACATGGGACTGCTGTCCGTGTGGGGGGGGCTTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                                                                    <  1:637718/150‑1 (MQ=255)
                                                 tGCGCCCCTTTCGCCACACACGGCCGCTCCCACGGGGCCACGCTGGCACACGGTAATTCCGTCCGTGTGGGCGCGGCCTTGTTTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                                                                    <  1:2056143/150‑1 (MQ=255)
                                                  gcgcCCCTTTCGCTACAAACGGCCGCACCCACTGGGCCACGCTGGAACATGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAag                                                                                                   <  1:1887977/150‑1 (MQ=38)
                                                       ccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGTGCGGGCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAGCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGcc                                                                                              >  2:3387294/1‑150 (MQ=255)
                                                         tttCGCGACACAAGGCCGCTCCCACCGGGCCACGCTGGCACAAGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                            <  2:1414615/150‑1 (MQ=255)
                                                         tttCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                            <  2:4162108/150‑1 (MQ=255)
                                                             GCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGCCCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGc                                                                                        <  2:1347374‑M1/150‑1 (MQ=255)
                                                                 CACACGGCCGCCCCCCCCGGGCCAGGCTGCCACATGGTATTGCTGTCCGTGTGGGCGGGGCGTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTccc                                                                                    >  1:1073024‑M1/1‑150 (MQ=255)
                                                                   acaaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                                                                   >  1:2924999‑M2/73‑150 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:4389398‑M1/150‑1 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  1:4016612‑M1/150‑1 (MQ=255)
                                                                     aaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcc                                                                                     <  1:5054035‑M2/76‑1 (MQ=255)
                                                                     AGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCAt                                                                                <  1:4194799‑M1/150‑1 (MQ=255)
                                                                      GGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                                                         >  1:984067‑M1/1‑124 (MQ=255)
                                                                      GGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               <  1:1090386‑M1/150‑1 (MQ=255)
                                                                         CGCTCCCACTGGCCCACGCTGCCAAAAGGTATTGCAGCCCGTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            >  1:1626298‑M1/1‑150 (MQ=255)
                                                                         CGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            <  2:418914‑M1/150‑1 (MQ=255)
                                                                            TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         <  2:4783201‑M1/150‑1 (MQ=255)
                                                                            TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         <  2:2220045‑M1/150‑1 (MQ=255)
                                                                             tcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcggTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         >  2:5052648‑M2/63‑150 (MQ=255)
                                                                             tcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         >  1:4464569‑M2/63‑150 (MQ=255)
                                                                             tcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         >  1:302735‑M2/63‑150 (MQ=255)
                                                                             CCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGTGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGTATCCCGCTCTACCCGAGCGCGGAAAACCGAGGCCCAAAGAGGACCACCGCTCCCCCATCAACCTcc                                                                        >  2:3071316‑M1/1‑150 (MQ=255)
                                                                             CCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTcc                                                                        <  1:185508‑M1/150‑1 (MQ=255)
                                                                              CCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                                                       <  2:4232705‑M1/150‑1 (MQ=255)
                                                                               ccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                                                       >  2:431796‑M2/61‑150 (MQ=255)
                                                                               ccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                                                       <  2:2924999‑M2/90‑1 (MQ=255)
                                                                               CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              >  1:4095023‑M1/1‑110 (MQ=255)
                                                                               CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              <  2:4095023‑M1/110‑1 (MQ=255)
                                                                               CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  2:5025293‑M1/150‑1 (MQ=255)
                                                                                cacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGACAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:4214171‑M2/60‑150 (MQ=255)
                                                                                cacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:1619928‑M2/60‑150 (MQ=255)
                                                                                cacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:4214285‑M2/60‑150 (MQ=255)
                                                                                cacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:3676768‑M2/60‑150 (MQ=255)
                                                                                actgggCCCCGCTGGCACTTGGTATTGCAGTCCGTGTGGTAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     <  1:572653/147‑1 (MQ=255)
                                                                                 acagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAc                                                                                           >  2:1095131‑M2/59‑128 (MQ=255)
                                                                                 acagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAc                                                                                           <  1:1095131‑M2/70‑1 (MQ=255)
                                                                                  AGGGCCAGGCTGGCAAATGGTATTGCAGTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactga                                                                   >  1:3040766‑M1/1‑150 (MQ=255)
                                                                                    GGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgact                                                                 >  2:3923522‑M1/1‑150 (MQ=255)
                                                                                       cACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACcgcg                                                                                                                     >  1:2956798/1‑95 (MQ=255)
                                                                                              gCACACGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:589059/150‑1 (MQ=255)
                                                                                              gCACACGGGAAGGCAGTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:588958/150‑1 (MQ=255)
                                                                                              gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       >  1:4730665/1‑150 (MQ=255)
                                                                                              gCACAAGGGAAGGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:589026/150‑1 (MQ=255)
                                                                                                acaAGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:5105305/150‑1 (MQ=255)
                                                                                                acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  2:1663877/1‑150 (MQ=255)
                                                                                                acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:4864934/150‑1 (MQ=255)
                                                                                                  aaGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                   >  2:418097/1‑150 (MQ=255)
                                                                                                     ggAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaa                                                >  2:2591761/1‑150 (MQ=255)
                                                                                                               tCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                       >  1:4416442/1‑150 (MQ=255)
                                                                                                               tCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                       >  1:1386990/1‑150 (MQ=255)
                                                                                                                   tgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAAcc                                   >  1:3014127‑M2/25‑150 (MQ=255)
                                                                                                                     tgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                                                         <  1:3482083‑M2/56‑1 (MQ=255)
                                                                                                                     tgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                                                         >  2:3482083‑M2/23‑78 (MQ=255)
                                                                                                                      gggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                >  2:3466702‑M2/22‑150 (MQ=255)
                                                                                                                        gagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGc                              <  1:418914‑M2/131‑1 (MQ=255)
                                                                                                                              ccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGAc                        <  1:4550068‑M2/137‑1 (MQ=255)
                                                                                                                                 tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATTCCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:4139802‑M2/140‑1 (MQ=255)
                                                                                                                                 tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  1:5202149‑M2/140‑1 (MQ=255)
                                                                                                                                 tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:4446219‑M2/140‑1 (MQ=255)
                                                                                                                                 tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  1:3825969‑M2/140‑1 (MQ=255)
                                                                                                                                       gcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACg               >  2:4936169‑M2/5‑150 (MQ=255)
                                                                                                                                       gcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACg               >  1:1026461‑M2/5‑150 (MQ=255)
                                                                                                                                          cagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  2:1812187/4‑109 (MQ=255)
                                                                                                                                          cagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  2:1812506/4‑109 (MQ=255)
                                                                                                                                          cagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:1812187/106‑1 (MQ=255)
                                                                                                                                          cagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:1812506/106‑1 (MQ=255)
                                                                                                                                             ggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTa         >  1:3422999/1‑150 (MQ=255)
                                                                                                                                               aGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt       >  2:573050/1‑150 (MQ=255)
                                                                                                                                               aGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt       >  2:2238748/1‑150 (MQ=255)
                                                                                                                                                gggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTc      <  2:4458079/150‑1 (MQ=255)
                                                                                                                                                  gCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCaa    >  2:4029436/1‑150 (MQ=255)
                                                                                                                                                  gCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCaa    >  1:823833/1‑150 (MQ=255)
                                                                                                                                                   cAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAAc   >  1:5103280/1‑150 (MQ=255)
                                                                                                                                                    aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  2:4857723/1‑150 (MQ=255)
                                                                                                                                                    aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATGTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  2:5383781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4895461‑4895757

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: