Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,635 A→T 5.3% intergenic (+116/‑92) 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,898,6350AT 5.3% 681.1 / 19.6 210intergenic (+116/‑92)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base A (102/92);  new base T (8/3);  total (114/96)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 218 total aligned reads displayed.

GACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898486‑4898784
                                                                                                                                                     |                                                                                                                                                     
gACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCa                                                                                                                                                       >  2:1379985/1‑150 (MQ=255)
  cTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCaaa                                                                                                                                                     >  2:4363716/1‑150 (MQ=255)
   tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAATGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTACCAGTGGGAGCGGCCGTGTGTCGCGGTTGGAGGGCAAAg                                                                                                                                                    >  1:5548632/1‑150 (MQ=255)
   tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAg                                                                                                                                                    >  2:2964453/1‑150 (MQ=255)
   tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAg                                                                                                                                                    >  2:510747/1‑150 (MQ=255)
   tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCCCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCAAAg                                                                                                                                                    >  2:2442266/1‑150 (MQ=255)
         gCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGGCACAAAGCCCCGCCCCCAGGGCCCCCCCGGGACAACGGGATGCAGTCCCAGTGGGGGCGGCCGTGTGTGGCGATTTGTGGGGAAAGcccacc                                                                                                                                              >  1:254239/1‑150 (MQ=255)
           tGCGGATCGAGACTTTGATGTGAGTTGGCCGGGGGCGCGTTGCGCCCCTTTCGCCACACACGGCCGCTCCCCCCGGGCCACGCTGGCAAACGGTAATGCAGTCCCTGTGGGAGGGGCCTTGTGCCGCGATTGGAGGGCTTTGccctccaa                                                                                                                                            <  1:5478770/150‑3 (MQ=255)
             cGGATCGAGACTTTGATCTGAGATGGCCGGGGGCGCGTTGCGCCCCTGTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGGCCCTGTGGGAGCGGCCTTGTGTCGCGATGGGAGGGCTTTGccctccaatc                                                                                                                                          <  2:2127867/150‑5 (MQ=38)
              gtaTCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCATTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGaa                                                                                                                                         <  2:1413763/148‑1 (MQ=255)
              ggATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCAGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCGCCCGGaa                                                                                                                                         >  1:996950/1‑150 (MQ=255)
                 tCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAc                                                                                                                                      <  2:3231428/150‑1 (MQ=255)
                   gagaCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAACGGCATGCAGTCCCAGTGGGAGCGGCCATGTGTCGCGATTGGAGGGCAAAGCCCACCCTGaaaaacg                                                                                                                                    >  2:2346229/1‑145 (MQ=38)
                      aCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctc                                                                                                                                 >  2:5499780/1‑150 (MQ=255)
                       cTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAACGCCCCGCCCACCGGGCCACCCGGGCCCCACGGCAAGCAGTCCCGGGGGGCGCGGCCGTGTGTCGCGCTTGGAGGGCAAACCCCACCCGGAAAAACGcaca                                                                                                                                >  1:2029786/1‑147 (MQ=255)
                                 agaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGAAACAATGCCGCTCCCACACGGCCACGCAGGCACAAGGGAAATCAGTCCCAGGTGGAGCGGCCATGTGTGGCGAATGGAGGGCAAAGCCCGCCCAGAAAACCGCACGACccgcccgc                                                                                                                      >  1:3850456/1‑150 (MQ=255)
                                       cggggggCGCGTTGCGCCCCTTTCGCCACACACGCCCGCTCCCACCGGGCCCCGGTGGCACACGGGAATGCAGTCCCTGTGGGTGGGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGaaaa                                                                                                                <  1:1282827/148‑1 (MQ=255)
                                         aggggCGCAACGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCGCCCGGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              >  2:410695/2‑150 (MQ=39)
                                         aggggCGCAACGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCAGGCACAAGGGAATGCACTCCCAGTGGGAGCGGCCGTGTGGCGCCATTGGAGGGCAAAGCCCTCCCAGAATACCGCTCTACCCGACCCCGGAAAAcc                                                                                                              >  1:6145668/2‑150 (MQ=255)
                                           gggCGCTTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACCGGGCCCCGCTGGAACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  1:2612962/150‑1 (MQ=255)
                                            ggCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGATGGCCGGCCCTGTGGGTGCGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  1:1550139/150‑1 (MQ=255)
                                            ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  2:5408204/150‑1 (MQ=255)
                                            ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  2:3913492/150‑1 (MQ=255)
                                             gcggGTTGCGCCCCTTTCGCGACACACGGCCGCGCCCACTGGGCCACGCTGGCACAAGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                                                          <  1:2814782/150‑1 (MQ=39)
                                             gcgcTTGGCCCCCCTTGCGCCACACCCGCCCGCCCCCACCGGGCCACGCGGGCACACGGGATTGCTGACCCTGGGGGAGCGCCCTTTGGACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                                                          <  1:3020727/150‑1 (MQ=255)
                                             gcgcGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGa                                                                                                                          <  2:5611806/134‑1 (MQ=255)
                                              cgcgTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCATTCCCTGTGGGAGCGGCATTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                                                         <  1:4531154/150‑1 (MQ=255)
                                                         ccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGcc                                                                                              >  2:2622941/1‑150 (MQ=255)
                                                           aatcgcgacacaaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                            >  1:5509602‑M2/82‑150 (MQ=255)
                                                              CGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                         >  1:1055577‑M1/1‑150 (MQ=255)
                                                              CGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                         >  2:6803033‑M1/1‑150 (MQ=255)
                                                               GCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCccactgc                                                                                        <  2:3935979‑M1/150‑4 (MQ=255)
                                                                    ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                                                                   <  1:3661883‑M1/150‑1 (MQ=255)
                                                                     CAAGGCCGCTCCCACTGGGCCACGCTGGCACATGGTAATGCTGTCCCTTTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTccccc                                                                                  >  1:6790720‑M1/1‑150 (MQ=255)
                                                                      AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:6340526‑M1/150‑1 (MQ=255)
                                                                      AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:754090‑M1/150‑1 (MQ=255)
                                                                      AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  1:442906‑M1/150‑1 (MQ=255)
                                                                      AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:4783050‑M1/150‑1 (MQ=255)
                                                                        GGCCGCTCCCACCGGGCCACGCTGGCACACGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               >  2:1246571‑M1/1‑150 (MQ=255)
                                                                        GGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               >  2:2107302‑M1/1‑150 (MQ=255)
                                                                           CGCTCCCCCAGGGCCACGCTGGCACAAGGGAATGCCGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            >  2:642844‑M1/1‑150 (MQ=255)
                                                                           CGCTCCCACCGGGCCACGCTGGCACATGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            >  2:4060079‑M1/1‑150 (MQ=255)
                                                                                 CACAGGGCCGCGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATGGGAGGGCtttgccctccaatc                                                                                                                                          <  2:4523461‑M1/82‑5 (MQ=255)
                                                                                 CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  1:5419756‑M1/150‑1 (MQ=255)
                                                                                 CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  2:5218009‑M1/150‑1 (MQ=255)
                                                                                  cccGGGCCACGCTGGCACACGGGATTGCAGTCCGTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     >  2:1767370‑M1/4‑150 (MQ=255)
                                                                                  ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTCCCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     <  2:1708856‑M1/150‑1 (MQ=255)
                                                                                    AGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactga                                                                   >  1:925412‑M1/1‑150 (MQ=255)
                                                                                      GGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccct                                                                                                                                                >  1:1079029‑M1/1‑71 (MQ=255)
                                                                                        ccACGCTGGCACAAGGGTATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGc                                                               <  2:1858363/150‑1 (MQ=255)
                                                                                            gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                                                           <  2:5873682/150‑1 (MQ=255)
                                                                                            gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                                                           >  1:2178612/1‑150 (MQ=255)
                                                                                             cTGGCACATGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCAt                                                          <  2:2387546/150‑1 (MQ=255)
                                                                                             cTGCCACACGGGAATGCAGTCCCTGGGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCAt                                                          <  1:1899919/150‑1 (MQ=255)
                                                                                                gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:5279986/150‑1 (MQ=255)
                                                                                                  acaAGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:6238417/150‑1 (MQ=255)
                                                                                                    atgggATTGCAGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                   <  1:3080127/148‑1 (MQ=255)
                                                                                                     aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  2:925153/1‑150 (MQ=255)
                                                                                                     aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  1:4797331/1‑150 (MQ=255)
                                                                                                        gAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaa                                               >  1:3977781/1‑150 (MQ=255)
                                                                                                           tGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAAtt                                            >  1:583685/1‑150 (MQ=255)
                                                                                                            gCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCACAAAttt                                           >  1:5676373/1‑150 (MQ=255)
                                                                                                            cgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAAttt                                           >  1:2104523‑M2/33‑150 (MQ=255)
                                                                                                             ctttCCCTGGGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  1:2477343/147‑1 (MQ=255)
                                                                                                             cAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  2:3977668/150‑1 (MQ=255)
                                                                                                             cAGTCCCTGTGGGAGCGGAATTTTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  2:877821/150‑1 (MQ=255)
                                                                                                              aGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg                                         >  1:4510203/1‑150 (MQ=255)
                                                                                                               gTCCGTGTGGGCGCGGCGTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        <  1:25697/150‑1 (MQ=255)
                                                                                                               gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        >  1:6537054/1‑150 (MQ=255)
                                                                                                               gTCCATGTGGGAGCGGCCTTGTGTCGCGATAGGATGGCAAAGCCCTCCCTGACTACCGCTCTCCCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAGCAACAAATTCgg                                        >  2:6779874/1‑149 (MQ=38)
                                                                                                                  cctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCaa                                     >  2:1219469‑M2/27‑150 (MQ=255)
                                                                                                                   ctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAAc                                    <  2:3661883‑M2/125‑1 (MQ=255)
                                                                                                                   ctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAAc                                    >  2:4859026‑M2/26‑150 (MQ=255)
                                                                                                                              gccttgggtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGa                         <  1:4305520‑M2/136‑1 (MQ=255)
                                                                                                                               ccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGAc                        <  2:4889756‑M2/137‑1 (MQ=255)
                                                                                                                                cttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACt                       >  1:4598731‑M2/13‑150 (MQ=255)
                                                                                                                                  tgtgtcgcgaTTGGAGGGCGAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:6293231‑M2/140‑1 (MQ=255)
                                                                                                                                  tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:2029786‑M2/140‑1 (MQ=255)
                                                                                                                                  tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCCCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  1:5499780‑M2/140‑1 (MQ=255)
                                                                                                                                   gtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTggg                    <  2:5185091‑M2/141‑1 (MQ=255)
                                                                                                                                      tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  2:3361301‑M2/7‑150 (MQ=255)
                                                                                                                                      tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:6619376‑M2/7‑150 (MQ=255)
                                                                                                                                      tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:1710606‑M2/7‑150 (MQ=255)
                                                                                                                                      tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:293796‑M2/7‑150 (MQ=255)
                                                                                                                                       cgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGTCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCAc                <  1:1509109‑M2/145‑1 (MQ=255)
                                                                                                                                         ttcagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACgg              >  1:3321766/6‑150 (MQ=255)
                                                                                                                                          tcagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             >  1:2286206/5‑150 (MQ=255)
                                                                                                                                          gaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             <  1:4910729‑M2/148‑1 (MQ=255)
                                                                                                                                          gaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             >  1:5101290‑M2/3‑150 (MQ=255)
                                                                                                                                            ttGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           <  1:885399/150‑1 (MQ=255)
                                                                                                                                            ttGGAGGGCAAAGCCCTCCCTGAATACCGCTCGACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           <  2:1533469/150‑1 (MQ=255)
                                                                                                                                            agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  2:4769353/3‑150 (MQ=255)
                                                                                                                                            agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  1:2460531/3‑150 (MQ=255)
                                                                                                                                            agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  2:6513801/3‑150 (MQ=255)
                                                                                                                                             gggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGtt          >  2:893555/2‑150 (MQ=255)
                                                                                                                                              ggAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTa         <  1:4168515/150‑1 (MQ=255)
                                                                                                                                               gAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTAc        >  2:274946/1‑150 (MQ=255)
                                                                                                                                                aGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt       <  1:2359468/150‑1 (MQ=255)
                                                                                                                                                 gggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGCGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTc      <  1:830743/150‑1 (MQ=255)
                                                                                                                                                  ggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                <  2:3703575/123‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898486‑4898784

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: