Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 703,038 G→A 5.2% intergenic (+32/‑227) 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,0380GA 5.2% 453.8 / 10.3 153intergenic (+32/‑227)PP_5SD/PP_05985S ribosomal RNA/paraquat‑inducible protein A
Reads supporting (aligned to +/- strand):  ref base G (74/71);  new base A (5/3);  total (79/74)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 156 total aligned reads displayed.

ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGCACA  >  NC_002947/702912‑703185
                                                                                                                              |                                                                                                                                                   
accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGa                                                                                                                              <  1:5440231/150‑1 (MQ=11)
  caccTGATCCCATCCCGAACTCAATAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaa                                                                                                                            <  2:4285533/150‑1 (MQ=17)
   accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaa                                                                                                                           >  2:3237736/1‑150 (MQ=32)
    ccTGATCCCATCCCGAACTCAGTAGTGAAAAGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTATGGGTTTTGTCTTTGGTGGGaaaaa                                                                                                                          >  1:169332/1‑150 (MQ=17)
      tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                                                                        >  1:31613/1‑150 (MQ=32)
        atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                                                                      <  2:2466395/150‑1 (MQ=33)
        atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                                                                      >  2:2615074/1‑150 (MQ=33)
         tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTg                                                                                                                     <  1:1052932/150‑1 (MQ=33)
         tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTg                                                                                                                     <  1:5343496/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)
              tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACt                                                                                                                >  1:4532099/1‑150 (MQ=33)
                  gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAGGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                                                                            <  1:260585/150‑1 (MQ=18)
                  gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                                                                            <  1:2954270/150‑1 (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)
                                aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                                                                              <  1:4437579/150‑1 (MQ=33)
                                aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                                                                              <  2:5125233/150‑1 (MQ=33)
                                 aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAg                                                                                             >  2:2380224/1‑150 (MQ=33)
                                 aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAg                                                                                             >  2:2380226/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)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                                           <  1:1573333/150‑3 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                                           >  2:1573333/1‑148 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                                                                           <  1:1573234/150‑3 (MQ=18)
                                                   tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGt                                                                           >  1:911062/1‑150 (MQ=33)
                                                     gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTg                                                                         <  1:5224317/150‑1 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                                      >  2:1976950/1‑150 (MQ=33)
                                                        tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                                                                      >  1:2064142/1‑150 (MQ=33)
                                                         ggggTTTCCGCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCTCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTGTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCgg                                                                     >  2:3517098/1‑150 (MQ=16)
                                                          gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAATCt                                                                                                                      >  2:4505774/1‑100 (MQ=18)
                                                                 cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCTCCGGCGTTGCCGGTGTTTGCg                                                             >  2:2293838/1‑150 (MQ=25)
                                                                 cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGACACCGGCAACGCCGgtgccatct                                                             <  1:5177217/150‑7 (MQ=11)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                            <  1:2199495/150‑1 (MQ=35)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                            <  1:3478408/150‑1 (MQ=35)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                            <  2:76859/150‑1 (MQ=35)
                                                                  ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg                                                            <  1:1298131/150‑1 (MQ=35)
                                                                   cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCggg                                                           <  2:495844/150‑1 (MQ=35)
                                                                     tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa                                                         >  2:140103/1‑150 (MQ=35)
                                                                     tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGAAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGTGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa                                                         >  2:140038/1‑150 (MQ=25)
                                                                        gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGTCACCGGCGTTGCCGGTGTTTGCGGGTAAAc                                                      <  1:4526834/150‑1 (MQ=25)
                                                                         agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAcc                                                     <  2:2030737/150‑1 (MQ=35)
                                                                            gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGc                                                  <  2:2547260/150‑1 (MQ=35)
                                                                             tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt                                                 <  2:4360223/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                                         <  2:5248830/150‑1 (MQ=37)
                                                                                     cGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                                         <  1:2324651/150‑1 (MQ=37)
                                                                                         aaGATTCATTTCGCAAATCCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTTTTTGCGGGTAAACCCGCTCCTACAAGGGAt                                     <  2:745062/150‑1 (MQ=255)
                                                                                         aaGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGAt                                     <  1:5336925/150‑1 (MQ=40)
                                                                                         aaGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGAt                                     >  2:3350715/1‑150 (MQ=40)
                                                                                          aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATg                                    >  1:1185305/1‑150 (MQ=40)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                                   >  2:1813848/1‑150 (MQ=255)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                                   >  2:308501/1‑150 (MQ=255)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                                   >  2:2309898/1‑150 (MQ=255)
                                                                                           gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATgg                                   <  2:2251855/150‑1 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                                 <  1:1010570/150‑1 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                                 >  1:1140336/1‑150 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                                 >  1:817681/1‑150 (MQ=255)
                                                                                             ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGGTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGt                                 <  2:15168/150‑1 (MQ=255)
                                                                                              tCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTg                                <  1:5426214/150‑1 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             >  2:1404353/1‑150 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             <  1:2309898/150‑1 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             >  1:2738515/1‑150 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             >  1:5273139/1‑150 (MQ=255)
                                                                                                 tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTGTGGGGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGCGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACCAGGGATGGGTGCac                             >  2:1796236/1‑149 (MQ=255)
                                                                                                   agacAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACg                           >  2:4506325/4‑150 (MQ=17)
                                                                                                     gCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGcc                         >  2:2887104/1‑150 (MQ=255)
                                                                                                     gCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACGGCTCGAATCACTGCGGGAGCGGCGATGGCACCGGCGGTGCCGGTGTTTGCGGGTAAACCCGCTCCTCCAAGGGATGGGGGCAACGcc                         >  2:2887277/1‑150 (MQ=14)
                                                                                                      cAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCca                        <  1:3517098/150‑1 (MQ=255)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                                         >  1:272087/1‑129 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                                         <  2:272148/129‑1 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                                         <  2:272114/129‑1 (MQ=11)
                                                                                                          aCCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAg                                         >  1:272148/1‑129 (MQ=11)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                                  >  1:2689497/1‑102 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                                  >  1:2689466/1‑102 (MQ=35)
                                                                                                            cccTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGgtgt                                                                  <  2:2689497/102‑1 (MQ=35)
                                                                                                              cTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAgggg                >  1:4979935/1‑150 (MQ=255)
                                                                                                                aTCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGATCCTACAAGGGATGGGTGCAACGCCAGGCAGGGgcg              >  1:3194212/1‑150 (MQ=255)
                                                                                                                  cTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTAGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCa            >  1:487549/1‑150 (MQ=255)
                                                                                                                   tGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGACAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCAc           >  2:5063971/1‑150 (MQ=255)
                                                                                                                      gcATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCAcccc        <  1:1229675/150‑1 (MQ=255)
                                                                                                                          gCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGca    >  1:2945023/1‑150 (MQ=255)
                                                                                                                            gggTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGcaca  >  2:3870385/2‑150 (MQ=255)
                                                                                                                            aGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTa                                             >  2:573983/1‑107 (MQ=38)
                                                                                                                            aGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGcaca  >  1:1474109/1‑150 (MQ=255)
                                                                                                                            aGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACCGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGcaca  >  2:118814/1‑150 (MQ=255)
                                                                                                                              gTAGGTGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             <  1:4387793/121‑1 (MQ=255)
                                                                                                                              gTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCaa                             >  2:4387793/1‑121 (MQ=255)
                                                                                                                              |                                                                                                                                                   
ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCTCCTACAAGGGATGGGTGCAACGCCAGGCAGGGGCGCACCCCCGCACA  >  NC_002947/702912‑703185

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: