Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,330,501 A→G 6.4% intergenic (+26/‑281) PP_5SE → / → PP_1157 5S ribosomal RNA/acetolactate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,330,5010AG 6.4% 527.2 / 16.3 171intergenic (+26/‑281)PP_5SE/PP_11575S ribosomal RNA/acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base A (80/80);  new base G (9/2);  total (89/82)
Fisher's exact test for biased strand distribution p-value = 5.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 174 total aligned reads displayed.

TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAACG  >  NC_002947/1330387‑1330650
                                                                                                                  |                                                                                                                                                     
 gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                   <  1:5178614/150‑1 (MQ=31)
  atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTg                                                                                                                  <  1:4119913/150‑1 (MQ=32)
    cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGa                                                                                                                <  2:1199080/150‑1 (MQ=33)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               >  1:1212538/1‑150 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               >  1:1212546/1‑150 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               <  2:450048/150‑1 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               <  1:2161362/150‑1 (MQ=34)
       atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCt                                                                                                             >  1:3733308/1‑150 (MQ=34)
       atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCt                                                                                                             >  1:4635489/1‑150 (MQ=34)
               cTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCgaga                                                                                                     <  2:28999/150‑1 (MQ=38)
                tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAATAAATTGGACCTGCTCGagag                                                                                                    >  2:4969181/1‑150 (MQ=255)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCCGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  >  2:2989221/1‑150 (MQ=255)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               >  1:1581629/1‑137 (MQ=34)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               <  2:1581629/137‑1 (MQ=34)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  >  1:4144551/1‑150 (MQ=40)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  <  2:394539/150‑1 (MQ=40)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  <  2:548511/150‑1 (MQ=40)
                   gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCgg                                                                                                 >  1:890491/1‑150 (MQ=255)
                   gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCgg                                                                                                 >  1:890642/1‑150 (MQ=255)
                   gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCgg                                                                                                 <  1:1933361/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  1:5413539/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  2:2494098/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  1:1505372/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  1:1505338/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  1:3903665/150‑1 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                              <  1:4853163/150‑1 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                              >  1:3469995/1‑150 (MQ=255)
                        gAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                            >  1:2986475/1‑150 (MQ=255)
                        gAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                            >  2:1978480/1‑150 (MQ=255)
                         aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCCTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATgg                                                                                           >  2:4093928/1‑150 (MQ=255)
                           aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCa                                                                                         <  1:1985145/150‑1 (MQ=255)
                              aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCg                                                                                      <  1:4031270/150‑1 (MQ=255)
                                gCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGc                                                                                    >  1:4974574/1‑150 (MQ=255)
                                  aTCGCCGATTGTAGTGTTGTGTTTCCCCATGTGAGAGTATTTCATCGTCAAGATTCATTTCGCAAAAACCCTATCTGCGCACTCAGGTAGGGTTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCTCCGGCtt                                                                                  <  2:5142483/150‑1 (MQ=17)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCCCAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                  >  1:3504694/1‑150 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                  <  1:4831657/150‑1 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                  >  1:1323739/1‑150 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                             <  2:4144551/150‑1 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTATTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                             <  1:1839227/150‑1 (MQ=255)
                                         aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                           <  2:1323739/150‑1 (MQ=255)
                                               gtgtGGGGGTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGATTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgc                                                                     <  1:5384965/150‑1 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                    >  1:3926337/1‑150 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                    >  1:4731466/1‑150 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   <  2:2350151/150‑1 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   <  2:2898838/150‑1 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   >  2:2712966/1‑150 (MQ=255)
                                                  tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCACGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                  >  2:1511125/1‑150 (MQ=255)
                                                       tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                             <  2:4478975/150‑1 (MQ=255)
                                                       tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                             <  2:5006399/150‑1 (MQ=255)
                                                          ccccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCg                                                          <  1:2411886/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                                <  1:2407863/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                                <  1:170347/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                                <  2:3138149/150‑1 (MQ=255)
                                                                     agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTATTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACagag                                               <  1:5001501/150‑1 (MQ=255)
                                                                         ggTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                           >  2:4345326/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                          >  2:1562691/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                          >  2:4897888/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                          >  2:4351385/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                          >  2:4351539/1‑150 (MQ=255)
                                                                            cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                           <  2:1762708/115‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCACAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAGGGCCCCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  2:1149329/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  2:3165932/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  2:3890621/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  1:1089875/150‑1 (MQ=255)
                                                                                gTCAAGATTCATGTCGCAACACCCCGATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAACAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  2:4971404/150‑1 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                                <  2:3101325/137‑1 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                                >  1:3101325/1‑137 (MQ=255)
                                                                                  cAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGc                                  >  1:2534755/1‑150 (MQ=255)
                                                                                    aGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCa                                <  1:1258668/150‑1 (MQ=255)
                                                                                     gATTCATTTCGCAAAACCCCTCTCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAg                               <  2:3435694/150‑1 (MQ=255)
                                                                                     gATTCATTTCGCAAAACCCCGATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAg                               <  1:3630667/150‑1 (MQ=255)
                                                                                      gttCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                              >  1:968782/2‑150 (MQ=255)
                                                                                      aTTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                              >  1:1042058/1‑150 (MQ=255)
                                                                                      aTTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                              >  2:2116771/1‑150 (MQ=255)
                                                                                      aTTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                              <  2:1537886/150‑1 (MQ=255)
                                                                                       ttCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGc                             >  2:2606336/1‑150 (MQ=255)
                                                                                         cATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTa                           <  2:2138781/150‑1 (MQ=255)
                                                                                         cATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTa                           <  2:2184458/150‑1 (MQ=255)
                                                                                          aTTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAg                          >  1:1806494/1‑150 (MQ=255)
                                                                                           tttCGCAACACCCCTCTCGGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGa                         <  2:3184772/150‑1 (MQ=255)
                                                                                           tttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGa                         >  1:329138/1‑150 (MQ=255)
                                                                                              gacAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATg                      >  2:5148277/3‑150 (MQ=255)
                                                                                               gCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATgg                     <  1:4799011/150‑1 (MQ=255)
                                                                                               gCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATgg                     >  1:2890228/1‑150 (MQ=255)
                                                                                                 aaaaCCCCTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc                   <  2:327920/150‑1 (MQ=255)
                                                                                                  aaaCCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCt                  <  1:5187144/150‑1 (MQ=255)
                                                                                                     ccccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCaa               >  2:2505127/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTTTACCCGc                                                         >  1:1048907/1‑107 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  2:3558267/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:214686/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:3171955/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:2701468/1‑150 (MQ=255)
                                                                                                        cTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgc            >  1:4994269/1‑150 (MQ=255)
                                                                                                         tATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgca           >  2:5308358/1‑150 (MQ=255)
                                                                                                          aTCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTCCAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcag          >  2:4372650/1‑150 (MQ=255)
                                                                                                           tCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcagc         >  2:4788046/1‑150 (MQ=255)
                                                                                                           tCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcagc         >  1:696852/1‑150 (MQ=255)
                                                                                                                gcacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTaa    >  1:2874088/1‑150 (MQ=255)
                                                                                                                 cacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACCGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAAc   <  2:4494398/150‑1 (MQ=255)
                                                                                                                  acgcaGGTAGGGGTTTTGTCTTTGTGGGCACAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAACg  >  2:3745974/1‑150 (MQ=255)
                                                                                                                  |                                                                                                                                                     
TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAACG  >  NC_002947/1330387‑1330650

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: