Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,335,450 A→G 6.6% 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,335,4500AG 6.6% 1288.9 / 51.6 421intergenic (+26/‑281)PP_5SE/PP_11575S ribosomal RNA/acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base A (173/220);  new base G (26/2);  total (199/222)
Fisher's exact test for biased strand distribution p-value = 1.81e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 432 total aligned reads displayed.

TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAACG  >  NC_002947/1335336‑1335599
                                                                                                                  |                                                                                                                                                     
 gATCCCATCCCGAACTCAGTAGTTAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                   <  1:1232568/150‑1 (MQ=14)
 gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                   <  2:8741675/150‑1 (MQ=31)
 gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAAtt                                                                                                                   <  2:2071718/150‑1 (MQ=31)
  atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTATGGGTTTTGACTTTGTGGGCAAAAAATTg                                                                                                                  <  1:6474422/150‑1 (MQ=17)
   tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTgg                                                                                                                 <  1:2593732/150‑1 (MQ=32)
    cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGa                                                                                                                <  1:1345589/150‑1 (MQ=33)
    cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGa                                                                                                                >  1:223682/1‑150 (MQ=33)
     ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGAc                                                                                                               <  2:192305/150‑1 (MQ=34)
       atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCt                                                                                                             >  1:4532486/1‑150 (MQ=34)
        tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTg                                                                                                            >  1:6336174/1‑150 (MQ=34)
        tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTg                                                                                                            <  1:3476001/150‑1 (MQ=34)
           cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTc                                                                                                         >  1:1117350/1‑150 (MQ=35)
             aaCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCga                                                                                                       >  2:7235883/1‑150 (MQ=37)
                tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGagag                                                                                                    >  1:6051945/1‑150 (MQ=39)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  >  2:3390305/1‑150 (MQ=40)
                  agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTATTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCg                                                                                                  <  2:4492940/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  2:1641059/150‑1 (MQ=255)
                    tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGga                                                                                                <  2:3632017/150‑1 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                              >  1:4298056/1‑150 (MQ=255)
                      gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGgaga                                                                                              >  1:8115468/1‑150 (MQ=255)
                       tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAt                                                                                             >  1:5644622/1‑150 (MQ=255)
                       tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAt                                                                                             >  2:313549/1‑150 (MQ=255)
                       tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAt                                                                                             >  2:8367475/1‑150 (MQ=255)
                       tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGAt                                                                                             >  2:8546088/1‑150 (MQ=255)
                        gAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATg                                                                                            <  2:4536940/150‑1 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCtt                                                                                  >  1:741169/1‑150 (MQ=255)
                                  aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATATTCAAGATTCATTTCGCAAAAAACAAATCTGAGAACTAAGGTAGGGGATTTGTCTTTTTGGGCAAAAAATTGGACCTGCTCGAGAGCGTAGATGGCACCGGCtt                                                                                  >  2:2946859/1‑150 (MQ=17)
                                     gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGc                                                                               <  1:6175372/150‑1 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                             <  1:7120361/150‑1 (MQ=255)
                                       cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCg                                                                             <  2:6810371/150‑1 (MQ=255)
                                         aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgt                                                                           <  1:101029/150‑1 (MQ=255)
                                          tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgtg                                                                          <  2:392437/150‑1 (MQ=255)
                                          tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGgtg                                                                          >  1:6561128/1‑150 (MQ=255)
                                            gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGtt                                                                        >  1:5311156/1‑150 (MQ=255)
                                            gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGtt                                                                        >  2:6097083/1‑150 (MQ=255)
                                              aGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcg                                                                      <  2:3526735/150‑1 (MQ=255)
                                                tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTcgcg                                                                    <  2:6965504/150‑1 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   <  1:7032590/150‑1 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   >  1:7674649/1‑150 (MQ=255)
                                                 gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCgg                                                                   <  2:5101521/150‑1 (MQ=255)
                                                  tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                  >  1:2769641/1‑150 (MQ=255)
                                                  tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCggg                                                                  <  1:3944977/150‑1 (MQ=255)
                                                   ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGt                                                                 <  1:3168389/150‑1 (MQ=255)
                                                    gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTa                                                                >  1:156867/1‑150 (MQ=255)
                                                       tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                             <  2:7921050/150‑1 (MQ=255)
                                                       tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAc                                                             <  1:8481666/150‑1 (MQ=255)
                                                        ttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAAcc                                                            <  2:6388680/150‑1 (MQ=255)
                                                          ccccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCg                                                          >  1:2530891/1‑150 (MQ=255)
                                                           cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGc                                                         <  1:4800206/150‑1 (MQ=255)
                                                            ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCt                                                        <  2:7081510/150‑1 (MQ=255)
                                                            ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCt                                                        >  1:8639641/1‑150 (MQ=255)
                                                             cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTc                                                       <  1:5882967/150‑1 (MQ=255)
                                                              aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTcc                                                      <  2:7670376/150‑1 (MQ=255)
                                                                gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGATCCTa                                                    <  2:2236023/150‑1 (MQ=255)
                                                                 tgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTAc                                                   >  2:5912768/1‑150 (MQ=255)
                                                                  gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACa                                                  <  1:2446239/150‑1 (MQ=255)
                                                                  gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACa                                                  >  2:7353013/1‑150 (MQ=255)
                                                                           tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGt                                         <  2:4279407/150‑1 (MQ=255)
                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGaa                                    <  1:690264/150‑1 (MQ=255)
                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAg                                   >  1:3575035/1‑150 (MQ=255)
                                                                                   aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGcc                                 >  1:7970914/1‑150 (MQ=255)
                                                                                   aaGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGcc                                 >  2:5793505/1‑150 (MQ=255)
                                                                                    aGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGGGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCa                                <  2:2726160/150‑1 (MQ=255)
                                                                                      aTTCATTTCGCAACACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAgg                              <  1:5770639/150‑1 (MQ=255)
                                                                                       ttCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTCCAGAGGATGGTGCGAAGCCAGGc                             >  2:5524557/1‑150 (MQ=255)
                                                                                       ttCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGc                             >  2:5313832/1‑150 (MQ=255)
                                                                                          aTTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAg                          >  1:4010867/1‑150 (MQ=255)
                                                                                           tttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGa                         >  1:2037986/1‑150 (MQ=255)
                                                                                            ttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGaa                        >  2:6413960/1‑150 (MQ=255)
                                                                                            ttCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGaa                        <  2:3827940/150‑1 (MQ=255)
                                                                                             tCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAAt                       <  2:8006851/150‑1 (MQ=255)
                                                                                              cGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATg                      <  1:5114193/150‑1 (MQ=255)
                                                                                              cGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATg                      <  1:1715541/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  2:1727030/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  1:8279326/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  1:3311418/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  1:907253/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  2:2503915/150‑1 (MQ=255)
                                                                                                cAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGc                    <  2:6953576/150‑1 (MQ=255)
                                                                                                 aaaaCCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc                   <  1:730547/150‑1 (MQ=255)
                                                                                                 aaaaCCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc                   >  1:6692593/1‑150 (MQ=255)
                                                                                                 aaaaCCCCTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGcc                   >  1:7779009/1‑150 (MQ=255)
                                                                                                     ccccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTc                 <  1:2841988/148‑1 (MQ=255)
                                                                                                     ccccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCaa               >  2:1534616/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:4442421/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  2:1113948/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              <  1:8407294/150‑1 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  2:4765496/1‑150 (MQ=255)
                                                                                                      cccTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:5466593/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:1849500/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:3061340/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:4183142/1‑150 (MQ=255)
                                                                                                      cccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAAc              >  1:4973946/1‑150 (MQ=255)
                                                                                                       ccTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACg             >  2:8129470/1‑150 (MQ=255)
                                                                                                        cTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgc            >  1:6295673/1‑150 (MQ=255)
                                                                                                        cTACCTGCGTGCGCAGATAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgc            >  1:3096936/1‑150 (MQ=255)
                                                                                                           tCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACgcagc         <  2:4010959/150‑1 (MQ=255)
                                                                                                             tGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCAt       >  2:2996665/1‑150 (MQ=255)
                                                                                                               cgcacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTa     >  2:2312344/1‑150 (MQ=255)
                                                                                                                 cacgcaGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAAc   <  1:6786512/150‑1 (MQ=255)
                                                                                                                  |                                                                                                                                                     
TGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTGTCTTTGTGGGCAAAAAATTGGACCTGCTCGAGAGCGGAGATGGCACCGGCTTCGCCGGTGTTCGCGGGTAAACCCGCTCCTACAGAGGATGGTGCGAAGCCAGGCTAGAATGGCCTCAACGCAGCATTAACG  >  NC_002947/1335336‑1335599

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: