Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 703,025 T→C 5.6% intergenic (+19/‑240) PP_5SD → / → PP_0598 5S ribosomal RNA/paraquat‑inducible protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947703,0250TC 5.6% 417.6 / 10.1 142intergenic (+19/‑240)PP_5SD/PP_05985S ribosomal RNA/paraquat‑inducible protein A
Reads supporting (aligned to +/- strand):  ref base T (67/67);  new base C (5/3);  total (72/70)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 142 total aligned reads displayed.

ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCG  >  NC_002947/702912‑703173
                                                                                                                 |                                                                                                                                                    
  caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaa                                                                                                                <  2:517560/150‑1 (MQ=32)
  caccTGATCCCATCCCGAACTCAATAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaa                                                                                                                <  2:4285533/150‑1 (MQ=17)
   accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaa                                                                                                               >  2:3237736/1‑150 (MQ=32)
    ccTGATCCCATCCCGAACTCAGTAGTGAAAAGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTATGGGTTTTGTCTTTGGTGGGaaaaa                                                                                                              >  1:169332/1‑150 (MQ=17)
      tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                                                            >  1:3375543/1‑150 (MQ=32)
       gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAAc                                                                                                           <  2:504060/150‑1 (MQ=32)
        atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                                                          <  2:2556147/150‑1 (MQ=33)
        atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                                                          <  2:2466395/150‑1 (MQ=33)
         tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTg                                                                                                         <  1:1052932/150‑1 (MQ=33)
           ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                                                       >  2:2324651/1‑150 (MQ=33)
           ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                                                       >  1:2030737/1‑150 (MQ=33)
           ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                                                       >  1:12532/1‑150 (MQ=18)
            catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGaa                                                                                                      >  2:1098830/1‑150 (MQ=33)
            catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGaa                                                                                                      >  1:4338326/1‑150 (MQ=18)
             atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAAc                                                                                                     >  1:947141/1‑150 (MQ=18)
                 cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGTGGTTTCACCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCAGCGCATGAAGGTAGGGGTTTTGTATTTGGTGGGAAAAAAACTGGAAATGCt                                                                                                 >  1:5320719/1‑150 (MQ=11)
                  gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAGGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                                                                <  1:260585/150‑1 (MQ=18)
                  gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                                                                <  1:2954270/150‑1 (MQ=33)
                    aCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGa                                                                                              <  1:982464/150‑1 (MQ=33)
                    aCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGa                                                                                              <  1:353626/150‑1 (MQ=33)
                      tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAAt                                                                                            >  1:2251855/1‑150 (MQ=33)
                      tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAAt                                                                                            >  1:225085/1‑150 (MQ=33)
                        agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCa                                                                                          <  1:1834159/150‑1 (MQ=33)
                          tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACt                                                                                        <  2:3859376/150‑1 (MQ=33)
                             tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGCAAAAAACTGGAACTGCTCGAATCACTGCg                                                                                     <  1:4218846/150‑1 (MQ=14)
                             tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCg                                                                                     <  2:3144628/150‑1 (MQ=33)
                             tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCg                                                                                     >  2:2346423/1‑150 (MQ=33)
                                aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                                                                  <  2:2396260/150‑1 (MQ=33)
                                aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                                                                  <  2:5125233/150‑1 (MQ=33)
                                aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                                                                  <  1:4437579/150‑1 (MQ=33)
                                 aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAg                                                                                 >  2:2380224/1‑150 (MQ=33)
                                     tGCATCGCCGATGGGAGTGGGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCAGGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGga                                                                             <  2:1976259/150‑1 (MQ=16)
                                        aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                                                            <  2:4367266/116‑1 (MQ=32)
                                        aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGAt                                                                          <  1:4122702/150‑1 (MQ=33)
                                         tCGGCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTGTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                                                         >  1:2521700/1‑150 (MQ=18)
                                         tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                                                         >  1:2547260/1‑150 (MQ=33)
                                         tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                                                         <  2:672102/150‑1 (MQ=33)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                               >  2:1573335/1‑148 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                               <  1:1573333/150‑3 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                               <  1:1573234/150‑3 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                               >  2:1573234/1‑148 (MQ=18)
                                                     gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTg                                                             <  1:2190286/150‑1 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                          >  2:1977031/1‑150 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                          >  1:2318648/1‑150 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                          >  1:2064142/1‑150 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                          >  2:3478408/1‑150 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTCCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                          >  1:4780890/1‑150 (MQ=18)
                                                         ggggTTTCCGCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCTCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTGTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCgg                                                         >  2:3517098/1‑150 (MQ=16)
                                                          gggTTTCCCCATGTGAGAGTAGGTCATCGTCAATATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                                                          <  1:4505774/100‑1 (MQ=25)
                                                                 cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCTCCGGCGTTGCCGGTGTTTGCg                                                 >  2:2293838/1‑150 (MQ=25)
                                                                 cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGACACCGGCAACGCCGgtgccatct                                                 <  1:5177505/150‑7 (MQ=11)
                                                                 cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGACACCGGCAACGCCGgtgccatct                                                 <  1:5177217/150‑7 (MQ=11)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTTGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                <  1:1327179/150‑1 (MQ=25)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                <  1:909790/150‑1 (MQ=35)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                <  1:1298131/150‑1 (MQ=35)
                                                                     tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa                                             >  2:140103/1‑150 (MQ=35)
                                                                     tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGAAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGTGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa                                             >  2:140038/1‑150 (MQ=25)
                                                                      gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTaa                                            >  2:2865183/1‑150 (MQ=35)
                                                                        gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGTCACCGGCGTTGCCGGTGTTTGCGGGTAAAc                                          <  1:4526834/150‑1 (MQ=25)
                                                                         agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAcc                                         <  2:911062/150‑1 (MQ=35)
                                                                         agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAcc                                         <  2:2030737/150‑1 (MQ=35)
                                                                            gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGc                                      <  2:3872234/150‑1 (MQ=35)
                                                                             tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt                                     <  2:4360225/150‑1 (MQ=35)
                                                                             tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt                                     <  2:2592802/150‑1 (MQ=35)
                                                                                gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGTGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCt                                  >  2:921910/1‑150 (MQ=25)
                                                                                  cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTAc                                >  1:1244382/1‑150 (MQ=35)
                                                                                     cGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             <  1:2324651/150‑1 (MQ=37)
                                                                                         aaGATTCATTTCGCAAATCCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTTTTTGCGGGTAAACCCGCTCCTACAAGGGAt                         <  2:745062/150‑1 (MQ=255)
                                                                                         aaGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGAt                         <  1:5336925/150‑1 (MQ=40)
                                                                                          aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATg                        >  1:1185305/1‑150 (MQ=40)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                       >  2:2309898/1‑150 (MQ=255)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                       <  2:3218691/150‑1 (MQ=255)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                       <  2:2251855/150‑1 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                     >  1:817681/1‑150 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                     <  1:1010570/150‑1 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGGTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                     <  2:15168/150‑1 (MQ=255)
                                                                                              tCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTg                    <  1:5426214/150‑1 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                 >  1:5273139/1‑150 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                 >  2:1404353/1‑150 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTGTGGGGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGCGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACCAGGGATGGGTGCac                 >  2:1796236/1‑149 (MQ=255)
                                                                                                   agacAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACg               >  2:4506325/4‑150 (MQ=17)
                                                                                                     gCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACGGCTCGAATCACTGCGGGAGCGGCGATGGCACCGGCGGTGCCGGTGTTTGCGGGTAAACCCGCTCCTCCAAGGGATGGGGGCAACGcc             >  2:2887277/1‑150 (MQ=14)
                                                                                                      cAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCca            <  1:3517098/150‑1 (MQ=255)
                                                                                                       aaaaCCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCcag           >  2:1843957/1‑150 (MQ=255)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             >  1:272087/1‑129 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             >  1:272148/1‑129 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             <  2:272087/129‑1 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             <  2:272148/129‑1 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                             >  1:272114/1‑129 (MQ=11)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                      >  1:2689497/1‑102 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                      <  2:2689497/102‑1 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                      >  1:2689466/1‑102 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                      <  2:2689466/102‑1 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTGTTGTCTTTGGTGGGAAAAAAACTGGAACTGCGCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACACGGGATGGGTGCAACGCcaggcagg      >  2:4870597/1‑150 (MQ=255)
                                                                                                              cTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAgggg    >  1:1128542/1‑150 (MQ=255)
                                                                                                              cTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAgggg    >  1:4979935/1‑150 (MQ=255)
                                                                                                               tATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAccc                                        >  1:298313/1‑113 (MQ=38)
                                                                                                               tATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAccc                                        <  2:298313/113‑1 (MQ=38)
                                                                                                                aTCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGgcg  >  1:5381900/1‑150 (MQ=255)
                                                                                                                aTCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGATCCTACAAGGGATGGGTGCAACGCCAGGCAGGGgcg  >  1:3194212/1‑150 (MQ=255)
                                                                                                                 |                                                                                                                                                    
ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCG  >  NC_002947/702912‑703173

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