Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,363,550 T→C intergenic (+206/‑109) RPPX_19410 → / → RPPX_19415 alkyl hydroperoxide reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,363,5500TC78.0% 720.3 / 205.4 362intergenic (+206/‑109)RPPX_19410/RPPX_19415alkyl hydroperoxide reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/59);  new base C (135/142);  total (158/204)
Fisher's exact test for biased strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 362 total aligned reads displayed.

CTGTGAAGATCCCCGTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGTTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCTAGTCCC  >  CP009974/4363403‑4363696
                                                                                                                                                    |                                                                                                                                                   
cTGTGAAGATCCCCGTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGCTCGGTGGTAGACACACATCAGAGGGTgg                                                                                                                                                  >  2:2375511/1‑151 (MQ=255)
      aGATCCCCGTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAg                                                                                                                                           <  1:4029152/151‑1 (MQ=255)
        gccgttgTAGGAGCGGGTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGgc                                                                                                                                         >  1:5276767/7‑151 (MQ=17)
        aTCCCCGTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAgg                                                                                                                                          <  1:5407580/151‑1 (MQ=255)
         tcccgttAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGgcc                                                                                                                                        <  2:819964/145‑1 (MQ=255)
         ccgttgTAGGAGCGGGTTTACCCGCGAAGAGGCCGGGGCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCCCCTACACATATCGGGGGTAGACACACATCGGGCGCTAGGCACAGgcg                                                                                                                                        >  2:4524801/6‑150 (MQ=14)
          cgttgTAGGAGCGGGTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGGAAACGCGCTCCTACACAGAGCGGTGGTAGGCACATATCGGAGGCTGGCCACAGgccg                                                                                                                                       >  2:1060407/5‑151 (MQ=14)
            gtggaGGGGGGGGTTGACCCGCGATGAGGCCGGGCCTGCGGGCCAGTCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATGCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGgccgcc                                                                                                                                     <  2:5638870/147‑1 (MQ=17)
             tgTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGccgccg                                                                                                                                    <  1:1038589/150‑1 (MQ=255)
               aaGGGGCGGGTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTa                                                                                                                                  <  2:59767/150‑1 (MQ=255)
                gggggCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGgccgttgta                                                                                                                                  <  2:3794332/150‑6 (MQ=18)
                  gggCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGCCCCAGGCCGCTGTAgg                                                                                                                                >  1:2475386/1‑151 (MQ=18)
                  gggCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGgccgttgtagg                                                                                                                                <  2:5779701/151‑8 (MQ=18)
                  gggCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGgccgttgtagg                                                                                                                                <  2:2465316/151‑8 (MQ=18)
                  gagCGGGTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGa                                                                                                                               <  1:140069/149‑1 (MQ=255)
                   ggCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGa                                                                                                                               <  1:76795/151‑1 (MQ=18)
                     cGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGc                                                                                                                             <  1:5969247/151‑1 (MQ=18)
                      gggTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCg                                                                                                                            <  2:4179703/151‑1 (MQ=18)
                      gggTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTGGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCg                                                                                                                            <  2:1550048/151‑1 (MQ=17)
                       ggTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCggg                                                                                                                          <  1:2046050/151‑1 (MQ=255)
                        gTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGt                                                                                                                         <  1:592380/151‑1 (MQ=255)
                        gTTTACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGt                                                                                                                         >  1:4398330/1‑151 (MQ=255)
                         ttCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGCTgg                                                                                                                                                  >  1:1531556/1‑126 (MQ=255)
                          ttaCCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGttt                                                                                                                       <  1:4839426/149‑1 (MQ=255)
                          ttaCCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGAGGTAGACACACATCGGAGGCTGGGCACAGGCCGCCGTCGCAGCGGGttt                                                                                                                       >  1:4642423/3‑151 (MQ=17)
                           cACCGGAGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGttt                                                                                                                       <  2:2020008/151‑1 (MQ=17)
                            aCCCGCGAAGAGGCGGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTa                                                                                                                      <  1:1005742/151‑1 (MQ=17)
                            aCCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGGTTGGCCCCAGGCGGGTTAAGGGGCGGGTTTa                                                                                                                      >  2:1724060/1‑151 (MQ=17)
                               cgcgAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTAccc                                                                                                                   <  1:5732894/151‑1 (MQ=17)
                                gcgAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTACCcg                                                                                                                  <  2:5885235/151‑1 (MQ=17)
                                 cgAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTGGCACACAGCCCGTTGTAGGAGCGGGCTTAccccc                                                                                                                 >  2:6529056/1‑149 (MQ=17)
                                  gAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTACCcgcg                                                                                                                <  2:3352491/151‑1 (MQ=17)
                                    agagGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGaaa                                                                                                             >  2:6019592/1‑151 (MQ=255)
                                       ggCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTACCCGCGaagag                                                                                                           <  2:3279911/151‑4 (MQ=17)
                                       ggCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTACCCGCGaagag                                                                                                           <  1:5591682/151‑4 (MQ=17)
                                         ccGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTCCACAGATCGGCGGTAGACACACATCGGAGGCTGGACACACGCCGCCGTAGGAGCGGGCTTACCCGCGAAAGGGcc                                                                                                        >  2:4771834/1‑151 (MQ=255)
                                         ccGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGATCGGGTTTACCCGCGAAAGGGcc                                                                                                        >  1:729358/1‑151 (MQ=255)
                                         ccGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                                                        >  2:6437784/1‑151 (MQ=255)
                                         ccGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                                                        <  1:3163234/151‑1 (MQ=255)
                                            ggctcgctgGCCAATCTCCCCCGGTCGTTCACCGTTCCTGCCGCTACACTCTTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGCCACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGt                                                                                                     <  2:142381/146‑1 (MQ=17)
                                            ggCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGCTGGACCCCGGCCGCCGGAGGAGGGGGGTTACCCGCCGAAGGGCCggg                                                                                                     >  2:4387748/1‑150 (MQ=18)
                                             gCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTg                                                                                                    >  1:4947701/1‑151 (MQ=255)
                                             gCCTGCGGGCTGATCGCCGCCAGGCGTTCACCGTTGCTGCCGGTACCGTAGTCGCGGGTAAACCCGCTCCTACACGGGGCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTg                                                                                                    <  2:305884/151‑1 (MQ=17)
                                              ccTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTgc                                                                                                   >  2:3179708/1‑151 (MQ=255)
                                              ccTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCCGAGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgggc                                                                                                    >  1:1951147/1‑149 (MQ=14)
                                               ctgctgGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTgca                                                                                                  >  1:4824011/1‑151 (MQ=255)
                                               ctgctgCCCAATCTCCCCCAGTCGTTCACCGTTGCTGCCGCTACCGCATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACAAACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTgca                                                                                                  <  1:1861959/151‑1 (MQ=18)
                                                  ctgGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTgcaagc                                                                                               >  1:1334762/1‑151 (MQ=255)
                                                  ctgGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGCATTCGCGGGTACACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTgcaagc                                                                                               <  1:3527574/151‑1 (MQ=255)
                                                       cAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGGCACCAGCCCGCGGAAGAGGGGGGTTACCCGCGGAAGGGGCGGGGCcagcaaagc                                                                                          >  2:1670868/1‑151 (MQ=18)
                                                        aaTCTCCACCAGTCGTTCGCCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCaagcaaagca                                                                                         <  1:2154685/151‑1 (MQ=255)
                                                        aaTCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCaagcaaagca                                                                                         <  1:6070937/151‑1 (MQ=255)
                                                        aaTCGCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCaagcaaagca                                                                                         <  1:332605/151‑1 (MQ=255)
                                                          tctcCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAgg                                                                                       >  1:4738931/1‑151 (MQ=255)
                                                              caccaGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGAccc                                                                                   >  2:3984856/1‑150 (MQ=255)
                                                                ccaGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGa                                                                                 <  1:6349099/151‑1 (MQ=255)
                                                                 caGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAGCCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGaa                                                                                <  2:2799487/151‑1 (MQ=255)
                                                                   gTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAgg                                                                              <  2:2411996/151‑1 (MQ=255)
                                                                          aCCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTAt                                                                       >  1:1771860/1‑151 (MQ=255)
                                                                             gTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGACCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGCAGCGGGGTTACCCGCGCAAGGGCCGGTGCAAGCAAAGCAGGCGCTGAAGGTGCCTATGCt                                                                    >  2:1389101/1‑151 (MQ=21)
                                                                              ttGCTGCCGCGCCCGTATTGGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCtt                                                                   <  2:1280387/151‑1 (MQ=255)
                                                                                  ggccggcaCCGTATGCGCGGGTAAACCCGCTCCTACACAGAGCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTg                                                               <  2:5283928/144‑1 (MQ=21)
                                                                                    ccGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCa                                                             <  1:6385743/151‑1 (MQ=255)
                                                                                     cGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAg                                                            <  1:237941/151‑1 (MQ=255)
                                                                                       cTACCGTATTCGCGGGTAAACCCGCACCTACGCAGGTGGGTGGTAGACACACATAGGCGGCTGGCAACAGGCCGGCGTAGGAGGGGGTTGACGCGCGAAAGGGCCGGTGCAAGCACAGGAGGGCCGGACAGTGTCTAGGCTTGGCGCAGAc                                                          >  2:631205/1‑151 (MQ=12)
                                                                                         gccGTAGTCGCGGGCAAACCCCCTACTACACCGATCGGTAGTAGACACCCATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACct                                                        <  1:5807794/150‑1 (MQ=25)
                                                                                         aCCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACct                                                        >  2:740478/1‑151 (MQ=255)
                                                                                          ccGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACctc                                                       >  1:550335/1‑151 (MQ=255)
                                                                                          ccGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACctc                                                       <  1:2952703/151‑1 (MQ=255)
                                                                                           cGTATTCGCGGGTAAGCCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACctct                                                      <  1:2817192/151‑1 (MQ=255)
                                                                                            gTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCtctc                                                     >  2:3430309/1‑151 (MQ=255)
                                                                                             tATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCa                                                    <  1:920322/151‑1 (MQ=255)
                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGATCGGTTACAGACACACATCGGCGCCTAGACACAGGCCGCCGCAGGAGGGGGTTTACCCGCCGAAGGACCCGTGCAACCACATCAGCGCCTGAAGGTGTCGTTGCTTACTCCTGAccccgcac                                                   >  2:4944049/1‑145 (MQ=9)
                                                                                               tgcgcgGGTAAACCGGCTCGTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACGCAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAgg                                                  <  2:4688186/149‑1 (MQ=255)
                                                                                                   cgGGTAAACCCGCCCCTACACAGATCGGTGGTAGACACACATCGGACGCGGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAACGGGTCTAGGCTTAGTCCAGACCCCTCAGGaacc                                              >  2:5664900/1‑149 (MQ=18)
                                                                                                   cgGGGAAACCCGCTCCTACACAGGTCGGGGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATc                                              <  1:6541517/151‑1 (MQ=255)
                                                                                                    gggTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATca                                             >  2:3955909/1‑151 (MQ=255)
                                                                                                    gggTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATca                                             >  1:2781421/1‑151 (MQ=255)
                                                                                                     ggTAACCCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCATGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATcac                                            <  2:6126681/151‑1 (MQ=255)
                                                                                                     ggTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATcac                                            <  1:3866252/151‑1 (MQ=255)
                                                                                                       tAAGCCCGCTCGTACACAGATCGGAGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCacac                                          <  1:6263440/151‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACACAGAGCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCacac                                          <  2:1872438/151‑1 (MQ=255)
                                                                                                          aCCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCa                                       >  2:3724654/1‑151 (MQ=255)
                                                                                                           cccGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCGCACCCAt                                      >  2:4331316/1‑151 (MQ=255)
                                                                                                             cGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACCCCCATGa                                    <  1:4143938/151‑1 (MQ=255)
                                                                                                              gCTCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGAc                                   <  2:885873/151‑1 (MQ=255)
                                                                                                                tCCTACACAGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGAccc                                 >  2:253309/1‑151 (MQ=255)
                                                                                                                    aaacaGATCGGTGGTAGAACCCCATCGGAGGCTGGAGACAGGCAGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGAt                             <  2:2641348/149‑1 (MQ=38)
                                                                                                                     cacaGATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATc                            >  2:222424/1‑151 (MQ=255)
                                                                                                                         gATCGGTGGTAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGcc                        >  2:5398510/1‑151 (MQ=255)
                                                                                                                          ttCGGTGGTAGACACACACCGGGGGCGGGACACAGGCCGCCGTAGGGGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCcg                       <  2:1196189/150‑1 (MQ=38)
                                                                                                                             ggtggtAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACGTGAAGGTGTCTATCCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGa                    >  2:6268144/1‑151 (MQ=255)
                                                                                                                                  tAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCt               >  1:6460884/1‑151 (MQ=255)
                                                                                                                                  tAGACACACATCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCt               >  1:3062646/1‑151 (MQ=255)
                                                                                                                                      cacccaTCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCtggtg           <  2:5485696/151‑1 (MQ=255)
                                                                                                                                         acaTCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCt        >  1:1658488/1‑151 (MQ=255)
                                                                                                                                           aTCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCTAg      >  1:1191722/1‑151 (MQ=255)
                                                                                                                                           aTCGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGCTGGCTAg      >  2:6197791/1‑151 (MQ=255)
                                                                                                                                             cGGAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCTAGTc    <  1:1724060/151‑1 (MQ=255)
                                                                                                                                              ggAGGCTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCTAGTcc   >  1:5153717/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGTGAAGATCCCCGTAGGGGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCAGTCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACAGATCGGTGGTAGACACACATCGGAGGTTGGACACAGGCCGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTGCAAGCAAAGCAGGACCTGAAGGTGTCTATGCTTACTGCAGACCTCTCAGGAATCACACCCATGACCCTGATCAGCCGCGAACCCTGGTGGCTAGTCCC  >  CP009974/4363403‑4363696

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