Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,330,500 C→T 6.4% intergenic (+25/‑282) 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,5000CT 6.4% 521.7 / 16.3 171intergenic (+25/‑282)PP_5SE/PP_11575S ribosomal RNA/acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base C (80/80);  new base T (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 = 1.00e+00

Only 100 of 173 total aligned reads displayed.

TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAAC  >  NC_002947/1330387‑1330649
                                                                                                                 |                                                                                                                                                     
tGATCCCAGCCCGAACTCAGGAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAATCCCCTCTCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAt                                                                                                                   <  2:2886299/150‑1 (MQ=2)
 gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCACAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                  <  1:4888957/150‑1 (MQ=14)
 gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                  <  1:5178614/150‑1 (MQ=31)
   tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTgg                                                                                                                <  2:3738773/150‑1 (MQ=32)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              <  1:2161362/150‑1 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              <  2:450048/150‑1 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              >  1:1212546/1‑150 (MQ=34)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              >  1:1212538/1‑150 (MQ=34)
       atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCt                                                                                                            >  1:4635489/1‑150 (MQ=34)
       atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCt                                                                                                            >  1:3733299/1‑150 (MQ=34)
               cTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCgaga                                                                                                    <  2:28999/150‑1 (MQ=38)
                tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAATAAATTGGACCTGCTCGagag                                                                                                   >  2:4969181/1‑150 (MQ=255)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              >  1:1581629/1‑137 (MQ=34)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                              <  2:1581629/137‑1 (MQ=34)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                 <  2:548511/150‑1 (MQ=40)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                 >  1:4144551/1‑150 (MQ=40)
                   gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCgg                                                                                                >  1:890617/1‑150 (MQ=255)
                   gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCgg                                                                                                <  1:1933361/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                               <  1:5413539/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                               <  1:1505372/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                               <  2:3125214/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                               <  1:1505338/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                               <  1:3903665/150‑1 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                             >  1:3469995/1‑150 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                             <  1:4853163/150‑1 (MQ=255)
                         aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCCTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATgg                                                                                          >  2:4093928/1‑150 (MQ=255)
                           aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCa                                                                                        <  1:1985145/150‑1 (MQ=255)
                              aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCg                                                                                     <  1:4031270/150‑1 (MQ=255)
                                  aTCGCCGATTGTAGTGTTGTGTTTCCCCATGTGAGAGTATTTCATCGTCAAGATTCATTTCGCAAAAACCCTATCTGCGCACTCAGGTAGGGTTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCTCCGGCtt                                                                                 <  2:5142483/150‑1 (MQ=17)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCCCAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                 >  1:3504694/1‑150 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                 >  2:2263000/1‑150 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                 <  1:4831657/150‑1 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                 >  1:1323739/1‑150 (MQ=255)
                                   tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTc                                                                                >  1:4773695/1‑150 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                            <  2:3041312/150‑1 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                            <  2:4144551/150‑1 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTATTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                            <  1:1839227/150‑1 (MQ=255)
                                         aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                          <  1:5209520/150‑1 (MQ=255)
                                         aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                          <  2:1323739/150‑1 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                   >  1:4731466/1‑150 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                   >  1:3926337/1‑150 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                   >  1:5178895/1‑150 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                  <  2:2350211/150‑1 (MQ=255)
                                                  tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                 >  2:1954115/1‑150 (MQ=255)
                                                  tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCACGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                 >  2:1511125/1‑150 (MQ=255)
                                                    gggTTTCACCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTGTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTa                                                               >  2:2541508/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)
                                                              aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTcc                                                     <  2:3468052/150‑1 (MQ=255)
                                                              aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTcc                                                     <  2:1166274/150‑1 (MQ=255)
                                                                gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTa                                                   <  2:2620964/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                               <  2:2174512/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                               <  1:2407863/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                               <  1:170347/150‑1 (MQ=255)
                                                                    gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                               >  1:2040331/1‑150 (MQ=255)
                                                                     agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTATTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACagag                                              <  1:5001501/150‑1 (MQ=255)
                                                                         ggTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                          >  2:4345326/1‑150 (MQ=255)
                                                                         ggTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATg                                          >  2:4006484/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                         >  2:4351539/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                         >  2:4351469/1‑150 (MQ=255)
                                                                          gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATgg                                         >  2:4897888/1‑150 (MQ=255)
                                                                            cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                          <  2:1762708/115‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                   <  2:3165932/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                   <  1:1089875/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                   <  2:3166140/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                   <  2:3890621/150‑1 (MQ=255)
                                                                                gTCAAGATTCATGTCGCAACACCCCGATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAACAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                   <  2:4971404/150‑1 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                           >  1:2957397/1‑93 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                           <  2:2957397/93‑1 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTc                                                      >  1:4105059/1‑130 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACaga                                               >  1:3101325/1‑137 (MQ=255)
                                                                                    aGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCa                               <  1:1258668/150‑1 (MQ=255)
                                                                                     gATTCATTTCGCAAAACCCCTCTCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAg                              <  2:3435694/150‑1 (MQ=255)
                                                                                     gATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAg                              <  1:851839/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)
                                                                                      aTTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAACTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                             >  2:2116723/1‑150 (MQ=255)
                                                                                         cATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTa                          <  2:2138781/150‑1 (MQ=255)
                                                                                          aTTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAg                         >  1:1806494/1‑150 (MQ=255)
                                                                                             tCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAAt                      <  2:3223459/150‑1 (MQ=255)
                                                                                              gacAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATg                     >  2:5148277/3‑150 (MQ=255)
                                                                                               gCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATgg                    >  1:2890228/1‑150 (MQ=255)
                                                                                                 aaaaCCCCTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc                  >  1:327920/1‑150 (MQ=255)
                                                                                                     ccccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCaa              >  2:2505127/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTTTACCCGc                                                        <  2:1048907/107‑1 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTTTACCCGc                                                        >  1:1048907/1‑107 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc             >  1:3171955/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGGAGGGGTTTTGTCTTTGTGGGCACAAAATTGGACCTGCTCGAGAGCGGAGAGGGCACCGGCTTCGCCGGTGTTCGCGGGTCCACCCGCTCCTACAGAGGCTGGTGCGACGCCCGGCTCGACGGGCCTCAAc             >  2:1543173/1‑150 (MQ=37)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc             >  1:2701468/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc             >  1:2532302/1‑150 (MQ=255)
                                                                                                        cTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgc           >  1:4994249/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)
                                                                                                            cTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcagca       <  2:1535520/150‑1 (MQ=255)
                                                                                                                gcacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTaa   >  1:2874088/1‑150 (MQ=255)
                                                                                                                 cacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACCGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAAc  <  2:4494398/150‑1 (MQ=255)
                                                                                                                 |                                                                                                                                                     
TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAAC  >  NC_002947/1330387‑1330649

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: