Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,002,710 T→C 9.5% intergenic (+210/‑1476) PP_4407 → / → PP_4409 gamma‑aminobutyrate transporter/site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,002,7100TC 9.5% 274.5 / 13.8 95intergenic (+210/‑1476)PP_4407/PP_4409gamma‑aminobutyrate transporter/site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base T (44/42);  new base C (7/2);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ACAAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCG  >  NC_002947/5002561‑5002859
                                                                                                                                                     |                                                                                                                                                     
acacGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCt                                                                                                                                                       <  2:3687100/150‑1 (MQ=255)
 cAAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTCTCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTg                                                                                                                                                      <  1:5398106/150‑1 (MQ=255)
 cAAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTg                                                                                                                                                      >  2:2971427/1‑150 (MQ=255)
    gATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGAc                                                                                                                                                   <  1:4494248/150‑1 (MQ=255)
      tGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTc                                                                                                                                                 <  1:3624602/150‑1 (MQ=255)
          aaTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTACCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGc                                                                                                                                             <  1:1806708/150‑1 (MQ=255)
           aTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCa                                                                                                                                            <  2:4930386/150‑1 (MQ=255)
               cGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAg                                                                                                                                        >  2:3595353/1‑150 (MQ=255)
                gagaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGc                                                                                                                                       >  2:2054328/1‑150 (MQ=255)
                 agaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCt                                                                                                                                      <  2:1289245/150‑1 (MQ=255)
                   aaTGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGc                                                                                                                                    >  1:2691710/1‑150 (MQ=255)
                             caggcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGc                                                                                                                          >  2:2737486/1‑150 (MQ=255)
                             caggcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGc                                                                                                                          <  2:5192063/150‑1 (MQ=255)
                             caggcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGc                                                                                                                          <  1:4464005/150‑1 (MQ=255)
                               ggcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGa                                                                                                                        <  2:4829350/150‑1 (MQ=255)
                                 caggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACt                                                                                                                      >  1:3326915/1‑150 (MQ=255)
                                           aCCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTg                                                                                                            >  1:2069212/1‑150 (MQ=255)
                                            cccTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTgg                                                                                                           <  2:4414743/150‑1 (MQ=255)
                                             ccTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGc                                                                                                          <  2:1308206/150‑1 (MQ=255)
                                             ccTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGc                                                                                                          >  1:947000/1‑150 (MQ=255)
                                               ttATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCtg                                                                                                        <  2:4501224/150‑1 (MQ=255)
                                                  tCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTgtgt                                                                                                     <  2:673495/150‑1 (MQ=255)
                                                  tCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTgtgt                                                                                                     <  2:673364/150‑1 (MQ=255)
                                                   cTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGtt                                                                                                    <  2:2025678/150‑1 (MQ=255)
                                                     ggggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTaa                                                                                                  <  2:1654486/150‑1 (MQ=255)
                                                      ggggGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAAc                                                                                                 <  1:706493/150‑1 (MQ=255)
                                                      gggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAAc                                                                                                 <  1:705689/150‑1 (MQ=255)
                                                          gCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCt                                                                                             >  2:4962447/1‑150 (MQ=255)
                                                           cctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTc                                                                                            >  2:3457346/1‑150 (MQ=255)
                                                           cctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTc                                                                                            >  2:3457535/1‑150 (MQ=255)
                                                            ctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                           >  1:1238764/1‑150 (MQ=255)
                                                             tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTccc                                                                                          >  2:2446214/1‑150 (MQ=255)
                                                             tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTccc                                                                                          >  2:3953031/1‑150 (MQ=255)
                                                                tgcctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATc                                                                                       >  2:914030/1‑150 (MQ=255)
                                                                 gcctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGATTTTTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATcc                                                                                      <  1:970115/150‑1 (MQ=255)
                                                                  cctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCt                                                                                     >  2:5349681/1‑150 (MQ=255)
                                                                  cctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCt                                                                                     <  1:5362675/150‑1 (MQ=255)
                                                                     gATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTgc                                                                                  <  2:1869947/150‑1 (MQ=255)
                                                                        cAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGcgcc                                                                               >  2:5017879/1‑150 (MQ=255)
                                                                        cAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGcgcc                                                                               >  2:5017863/1‑150 (MQ=255)
                                                                                cGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTa                                                                       <  2:75984/150‑1 (MQ=255)
                                                                                 gtggtgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAg                                                                      <  2:2420938/150‑1 (MQ=255)
                                                                                   ggtgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGt                                                                    <  1:2253248/150‑1 (MQ=255)
                                                                                    gtgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGtt                                                                   >  1:1802817/1‑150 (MQ=255)
                                                                                      gCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                 <  1:4777876/150‑1 (MQ=255)
                                                                                      gCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                 <  1:4777834/150‑1 (MQ=255)
                                                                                          tATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTAAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTAc                                                             <  1:5339797/150‑1 (MQ=255)
                                                                                            tGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGa                                                           >  1:2553689/1‑150 (MQ=255)
                                                                                            tGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGa                                                           >  1:2628660/1‑150 (MQ=255)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCTGCAGGAGTCCAGACTCCTGc                                                                                                                                             >  2:461863/1‑64 (MQ=37)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCTGCAGGAGTCCAGACTCCTGc                                                                                                                                             <  1:461863/64‑1 (MQ=37)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATc                                                       >  1:3614794/1‑150 (MQ=255)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATc                                                       <  2:1443894/150‑1 (MQ=255)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATc                                                       <  2:1443904/150‑1 (MQ=255)
                                                                                                tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATc                                                       >  2:4865496/1‑150 (MQ=255)
                                                                                                   gATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCg                                                    >  1:571465/1‑150 (MQ=255)
                                                                                                    aTTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCga                                                   <  1:4909459/150‑1 (MQ=255)
                                                                                                      tttCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCgaga                                                 <  1:4622073/150‑1 (MQ=255)
                                                                                                          cGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACt                                             >  1:1299981/1‑150 (MQ=255)
                                                                                                          cGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACt                                             <  2:2250911/150‑1 (MQ=255)
                                                                                                          cGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACt                                             >  1:1299720/1‑150 (MQ=255)
                                                                                                           gTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTg                                            >  1:3667809/1‑150 (MQ=255)
                                                                                                            tAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGa                                           <  1:3448655/150‑1 (MQ=255)
                                                                                                            tAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCAATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCAGTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGa                                           >  1:5044531/1‑150 (MQ=255)
                                                                                                                  agtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGAc                                     >  1:987725/1‑150 (MQ=255)
                                                                                                                    taCTTCTTCGTATTGTTGTTTTCCGAAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTAAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGc                                   <  1:927618/150‑1 (MQ=255)
                                                                                                                    taCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGc                                   >  2:3310548/1‑150 (MQ=255)
                                                                                                                    taCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGc                                   >  1:3840407/1‑150 (MQ=255)
                                                                                                                    taCTTCTTCGTAGAGTTGTTTTACGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATACTTGCGCCGCTCCCTAGGATCTGTACGAAATCCCGAGAAACTGAATCGACGc                                   >  2:3310555/1‑150 (MQ=255)
                                                                                                                      cttcttCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGccc                                 <  1:2637828/150‑1 (MQ=255)
                                                                                                                       ttcttcGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCga                                                   >  2:4447940/1‑131 (MQ=255)
                                                                                                                       ttcttcGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCga                                                   <  1:4447940/131‑1 (MQ=255)
                                                                                                                         cttcGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTc                              >  1:521328/1‑150 (MQ=255)
                                                                                                                          ttcGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCg                             >  1:3621082/1‑150 (MQ=255)
                                                                                                                             gTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGa                          >  2:1216862/1‑150 (MQ=255)
                                                                                                                              tAGTGTTGTTTTCCGCTGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAg                         <  1:4104711/150‑1 (MQ=255)
                                                                                                                              tAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAg                         >  2:8352/1‑150 (MQ=255)
                                                                                                                                  gttgttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGa                     <  2:1089716/150‑1 (MQ=255)
                                                                                                                                   ttgttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGaa                    >  1:822387/1‑150 (MQ=255)
                                                                                                                                      ttttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAAtt                 >  1:1536238/1‑150 (MQ=255)
                                                                                                                                      ttttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAAtt                 >  2:4546405/1‑150 (MQ=255)
                                                                                                                                        ttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt               >  2:1003324/1‑150 (MQ=255)
                                                                                                                                        ttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt               >  2:1585986/1‑150 (MQ=255)
                                                                                                                                        ttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt               >  1:3349325/1‑150 (MQ=255)
                                                                                                                                        ttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt               <  2:1412357/150‑1 (MQ=255)
                                                                                                                                        ttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt               <  1:2971427/150‑1 (MQ=255)
                                                                                                                                         tctgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                 <  2:4606043/96‑1 (MQ=39)
                                                                                                                                         tctgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                 >  1:4606043/4‑99 (MQ=39)
                                                                                                                                           tgCAGGAGTCCAGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTa            >  2:849872/2‑150 (MQ=255)
                                                                                                                                           cGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTa            <  2:521328/150‑1 (MQ=255)
                                                                                                                                            gCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCACATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAg           >  1:2376496/1‑150 (MQ=255)
                                                                                                                                            gCAGGAGTCCAGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAg           >  2:1482251/1‑150 (MQ=255)
                                                                                                                                             cAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAgg          >  1:3086201/1‑150 (MQ=255)
                                                                                                                                                   tCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTtctc    >  2:626884/1‑150 (MQ=255)
                                                                                                                                                   tCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTtctc    >  1:3000416/1‑150 (MQ=255)
                                                                                                                                                    cTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTcc   <  1:1289055/150‑1 (MQ=255)
                                                                                                                                                     tGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCg  >  2:1510047/1‑150 (MQ=255)
                                                                                                                                                     tGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCg  >  1:2706236/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCG  >  NC_002947/5002561‑5002859

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