Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,045,302 G→T 100% intergenic (+86/‑103) yihI → / → hemN activator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,045,3020GT98.1% 713.9 / 3.3 227intergenic (+86/‑103)yihI/hemNactivator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (134/88);  total (139/88)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 230 total aligned reads displayed.

AGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAG  >  CP009273/4045159‑4045446
                                                                                                                                               |                                                                                                                                                
 gaagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCtgttgt                                                                                                                                             <  2:925896/148‑1 (MQ=255)
 gaagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCtgttgt                                                                                                                                             <  1:1288781/148‑1 (MQ=255)
     agGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcgcg                                                                                                                                         >  1:3386688/1‑148 (MQ=255)
        aagaagACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCgcggcg                                                                                                                                      <  2:3408064/148‑1 (MQ=255)
           aagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACa                                                                                                                                   >  1:1377410/1‑148 (MQ=255)
            agaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaa                                                                                                                                  >  1:1553393/1‑148 (MQ=255)
             gaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaaa                                                                                                                                 >  1:1436691/1‑148 (MQ=255)
             gaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaaa                                                                                                                                 >  2:2412416/1‑148 (MQ=255)
             gaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTAGTTAAAATAAAGCTGTTGTCGCGGCGACaaa                                                                                                                                 >  2:599249/1‑148 (MQ=255)
                   aaGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGc                                                                                                                           >  1:889870/1‑148 (MQ=255)
                    aGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCt                                                                                                                          >  1:533442/1‑148 (MQ=255)
                    aGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCt                                                                                                                          >  1:2456331/1‑148 (MQ=255)
                     gCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATCTGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTg                                                                                                                         >  2:2185602/1‑148 (MQ=255)
                      cAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTgc                                                                                                                        >  1:2208009/1‑148 (MQ=255)
                      cAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTgc                                                                                                                        >  1:2103756/1‑148 (MQ=255)
                      cAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTGATATGTTATCTCGTGTGGTTCTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTgc                                                                                                                        >  1:3671740/1‑148 (MQ=255)
                         gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTACACTACAGCTGTTGTCGCGGCGACACAAGTGGCTGCGCa                                                                                                                     >  2:2338385/1‑148 (MQ=255)
                         gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCa                                                                                                                     >  1:3197485/1‑148 (MQ=255)
                         gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCa                                                                                                                     >  2:2804935/1‑148 (MQ=255)
                         gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGGTCTAGGGTATCTGGTGTGGGTGGTGGGTGAGCTACAGCGGTTGTCGCGGGGGCGTGAGTGGCTgcgcg                                                                                                                     >  1:3585911/1‑147 (MQ=255)
                         gaagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCCTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCa                                                                                                                     >  2:891343/1‑148 (MQ=255)
                          aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCaa                                                                                                                    >  1:1715381/1‑148 (MQ=255)
                            gaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAAcc                                                                                                                  >  2:1513104/1‑148 (MQ=255)
                             aCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCg                                                                                                                 >  1:223739/1‑148 (MQ=255)
                              cATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCTTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCgg                                                                                                                >  1:2551444/1‑148 (MQ=255)
                               atgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGc                                                                                                               >  2:531985/1‑148 (MQ=255)
                               atgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGc                                                                                                               >  1:2471997/1‑148 (MQ=255)
                                tgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCt                                                                                                              >  2:634214/1‑148 (MQ=255)
                                  atgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtga                                                                                                            >  2:275232/1‑148 (MQ=255)
                                  atgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtga                                                                                                            >  2:2202413/1‑148 (MQ=255)
                                   tgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgat                                                                                                           >  2:2158497/1‑148 (MQ=255)
                                    gCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                          >  2:2841253/1‑148 (MQ=255)
                                    gCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                          >  2:2838863/1‑148 (MQ=255)
                                         tgctgcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCg                                                                                                     >  2:2250770/1‑148 (MQ=255)
                                          gctgcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGa                                                                                                    >  1:728434/1‑148 (MQ=255)
                                            tgcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGaaa                                                                                                  <  1:429964/148‑1 (MQ=255)
                                                 ggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                             >  1:2560975/1‑148 (MQ=255)
                                                 ggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                             >  2:518136/1‑148 (MQ=255)
                                                 ggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                             >  1:3802019/1‑148 (MQ=255)
                                                         aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgccg                                                                                     >  1:1830271/1‑148 (MQ=255)
                                                           ggATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGt                                                                                   >  1:2167260/1‑148 (MQ=255)
                                                           ggATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGt                                                                                   <  2:2716395/148‑1 (MQ=255)
                                                                 ccgccgTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCg                                                                             >  1:1499932/1‑148 (MQ=255)
                                                                     cgTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgc                                                                         <  1:303554/148‑1 (MQ=255)
                                                                     cgTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgc                                                                         <  1:303769/148‑1 (MQ=255)
                                                                        gggCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAgccgc                                                                      <  2:1906217/148‑1 (MQ=255)
                                                                           ccTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCag                                                                   >  2:3690283/1‑148 (MQ=255)
                                                                             ttCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCagag                                                                 >  2:2007222/1‑148 (MQ=255)
                                                                                    tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCAt                                                          >  2:1666993/1‑148 (MQ=255)
                                                                                    tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCAt                                                          >  2:3212602/1‑148 (MQ=255)
                                                                                        ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                      <  2:2903285/148‑1 (MQ=255)
                                                                                            tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  2:3365991/148‑1 (MQ=255)
                                                                                            tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  2:121208/148‑1 (MQ=255)
                                                                                            tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  2:3365712/148‑1 (MQ=255)
                                                                                             aaCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAg                                                 <  1:1071753/148‑1 (MQ=255)
                                                                                              aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                >  2:456733/1‑148 (MQ=255)
                                                                                               cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                               <  1:2083125/148‑1 (MQ=255)
                                                                                               cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                               >  1:3672920/1‑148 (MQ=255)
                                                                                                ccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGa                                              >  1:298755/1‑148 (MQ=255)
                                                                                                  ttCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGc                                            <  2:1140489/148‑1 (MQ=255)
                                                                                                   tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCa                                           <  2:2280410/148‑1 (MQ=255)
                                                                                                   tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCa                                           >  2:3010458/1‑148 (MQ=255)
                                                                                                     cGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATg                                         <  1:2239525/148‑1 (MQ=255)
                                                                                                     cGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATg                                         <  1:2187822/148‑1 (MQ=255)
                                                                                                      agTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                        <  1:3301863/147‑1 (MQ=255)
                                                                                                            atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTAc                                  <  2:1830271/148‑1 (MQ=255)
                                                                                                            atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTAc                                  >  1:1386986/1‑148 (MQ=255)
                                                                                                            atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTAc                                  >  1:754252/1‑148 (MQ=255)
                                                                                                              gTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAg                                >  2:557830/1‑148 (MQ=255)
                                                                                                                 aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             <  2:2278349/148‑1 (MQ=255)
                                                                                                                 aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             <  2:2471997/148‑1 (MQ=255)
                                                                                                                 aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             <  2:1304389/148‑1 (MQ=255)
                                                                                                                 aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             <  2:1580129/148‑1 (MQ=255)
                                                                                                                  tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTTCGCCGTAGCCGCCAGAGACTCAATCGTAAGGAGTGAGCATGTCTGTACAGCaaa                            <  2:3495251/148‑1 (MQ=255)
                                                                                                                  tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                            <  2:1430738/148‑1 (MQ=255)
                                                                                                                  tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                            <  2:1818336/148‑1 (MQ=255)
                                                                                                                  tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                            <  1:2531590/148‑1 (MQ=255)
                                                                                                                  tCTGGTGTGGTTATTCGTTAAACTACAGCTGGTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                            <  2:2883106/148‑1 (MQ=255)
                                                                                                                   cTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                           >  2:3721711/1‑148 (MQ=255)
                                                                                                                        gtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACt                      >  1:2326414/1‑148 (MQ=255)
                                                                                                                          ggTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                    >  1:963861/1‑148 (MQ=255)
                                                                                                                             tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                 >  2:177914/1‑148 (MQ=255)
                                                                                                                             tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                 >  2:177330/1‑148 (MQ=255)
                                                                                                                              aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                <  2:195554/148‑1 (MQ=255)
                                                                                                                              aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                <  2:1083490/148‑1 (MQ=255)
                                                                                                                               ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             <  1:3217625/134‑1 (MQ=255)
                                                                                                                                tCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTg              >  2:3767414/1‑148 (MQ=255)
                                                                                                                                 cGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTgg             >  1:2857657/1‑148 (MQ=255)
                                                                                                                                   ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           <  1:1274962/148‑1 (MQ=255)
                                                                                                                                   ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           <  2:407432/148‑1 (MQ=255)
                                                                                                                                   ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           <  2:889056/148‑1 (MQ=255)
                                                                                                                                    tAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc          <  2:3266265/148‑1 (MQ=255)
                                                                                                                                      aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg        <  2:3591513/148‑1 (MQ=255)
                                                                                                                                       aCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa       <  2:3641775/148‑1 (MQ=255)
                                                                                                                                       aCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa       <  1:2412416/148‑1 (MQ=255)
                                                                                                                                        cTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt      >  1:3770682/1‑148 (MQ=255)
                                                                                                                                        cTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt      >  1:1292021/1‑148 (MQ=255)
                                                                                                                                         tACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATc     >  2:28113/1‑148 (MQ=255)
                                                                                                                                          aCAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATcc    >  1:1948930/1‑148 (MQ=255)
                                                                                                                                          aCAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATcc    >  1:678817/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAG  >  CP009273/4045159‑4045446

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