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,3020GT96.1% 539.2 / 12.2 180intergenic (+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 (4/3);  new base T (93/80);  total (97/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 182 total aligned reads displayed.

GAAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAG  >  CP009273/4045157‑4045446
                                                                                                                                                 |                                                                                                                                                
gaagaaGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGt                                                                                                                                                <  2:380106/148‑1 (MQ=255)
  agaagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTg                                                                                                                                              <  2:1927990/148‑1 (MQ=255)
   gaagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCtgttgt                                                                                                                                             >  2:743738/1‑148 (MQ=255)
   gaagaGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCtgttgt                                                                                                                                             >  1:1906039/1‑148 (MQ=255)
     agagGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcg                                                                                                                                           >  2:3508891/1‑148 (MQ=255)
       agGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTcgcg                                                                                                                                         >  1:1985849/1‑148 (MQ=255)
        ggAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCgcgg                                                                                                                                        >  1:784029/1‑148 (MQ=255)
        ggAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCgcgg                                                                                                                                        >  1:783722/1‑148 (MQ=255)
           agaagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATAGGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGa                                                                                                                                     >  2:3544715/1‑148 (MQ=255)
            gaagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGAc                                                                                                                                    >  2:519437/1‑148 (MQ=255)
             aagaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACa                                                                                                                                   >  2:1837143/1‑148 (MQ=255)
               gaCGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACaca                                                                                                                                 >  2:2588630/1‑146 (MQ=255)
                 cGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAg                                                                                                                               <  2:450458/148‑1 (MQ=255)
                  gagaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGt                                                                                                                              <  1:2356165/148‑1 (MQ=255)
                  gagaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGt                                                                                                                              <  1:1309837/148‑1 (MQ=255)
                    gaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTgg                                                                                                                            >  1:2662169/1‑148 (MQ=255)
                    gaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTgg                                                                                                                            >  2:2669813/1‑148 (MQ=255)
                    gaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTgg                                                                                                                            <  2:784029/148‑1 (MQ=255)
                    gaAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTgg                                                                                                                            >  2:381819/1‑148 (MQ=255)
                            aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCaa                                                                                                                    >  2:2227850/1‑148 (MQ=255)
                            aagaCATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCaa                                                                                                                    <  1:3267726/148‑1 (MQ=255)
                                cATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCgg                                                                                                                >  1:3401821/1‑148 (MQ=255)
                                 atgatgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGc                                                                                                               >  1:164523/1‑148 (MQ=255)
                                    atgCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtga                                                                                                            <  1:1376976/148‑1 (MQ=255)
                                      gCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                          >  2:1510317/1‑148 (MQ=255)
                                      gCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatg                                                                                                          >  2:3272410/1‑148 (MQ=255)
                                       cGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCtgatga                                                                                                         >  2:215343/1‑148 (MQ=255)
                                          ctgctgCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGAccc                                                                                                      >  1:2642951/1‑148 (MQ=255)
                                             ctgcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGaa                                                                                                   <  2:2907584/148‑1 (MQ=255)
                                               gcGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAAc                                                                                                 <  2:2622595/148‑1 (MQ=255)
                                                   ggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                             >  1:2253499/1‑148 (MQ=255)
                                                   ggCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGc                                                                                             >  1:3539204/1‑148 (MQ=255)
                                                     cAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCAt                                                                                           <  1:2164981/148‑1 (MQ=255)
                                                     cAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCAt                                                                                           >  1:1356119/1‑148 (MQ=255)
                                                        cTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                        >  1:2814158/1‑148 (MQ=255)
                                                        cTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                        >  2:671395/1‑148 (MQ=255)
                                                        cTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccg                                                                                        <  1:3507090/148‑1 (MQ=255)
                                                         tAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgc                                                                                       <  1:436661/148‑1 (MQ=255)
                                                           aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgccg                                                                                     <  2:1025303/148‑1 (MQ=255)
                                                           aCGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATccgccg                                                                                     <  1:1682970/148‑1 (MQ=255)
                                                             ggATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGt                                                                                   <  2:3052756/148‑1 (MQ=255)
                                                                ttGCCGCCGTGGGCCGTCCCGTTCTGCTGATCACCCTTCCGGTTATATGTTCTCTGGTGTGGTTAGTCGTTAAACTCCAGCTGTTGTCGCGGCGACAAAAGGGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACg                                                                                <  1:2300551/148‑1 (MQ=255)
                                                                 tGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGc                                                                               >  2:568511/1‑148 (MQ=255)
                                                                  gccgccGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGcc                                                                              >  2:1121038/1‑148 (MQ=255)
                                                                  gccgccGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGcc                                                                              >  2:1121160/1‑148 (MQ=255)
                                                                             ccTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCag                                                                   <  1:3391112/148‑1 (MQ=255)
                                                                                     ttCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCa                                                           <  2:247053/148‑1 (MQ=255)
                                                                                     ttCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCa                                                           <  2:3001443/148‑1 (MQ=255)
                                                                                     ttCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCa                                                           <  2:3001971/148‑1 (MQ=255)
                                                                                      tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTGCAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                          >  1:932168/1‑148 (MQ=255)
                                                                                      tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                          >  2:3501941/1‑148 (MQ=255)
                                                                                      tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCAt                                                          >  1:2709096/1‑148 (MQ=255)
                                                                                          ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                      >  1:1111545/1‑148 (MQ=255)
                                                                                          ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGa                                                      <  1:3197331/148‑1 (MQ=255)
                                                                                            gATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAg                                                    <  1:3416458/148‑1 (MQ=255)
                                                                                              tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  2:1704355/148‑1 (MQ=255)
                                                                                              tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  1:381819/148‑1 (MQ=255)
                                                                                              tAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                  <  1:961453/148‑1 (MQ=255)
                                                                                                aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                >  2:1918071/1‑148 (MQ=255)
                                                                                                aCCCTTCCGGTTATATGTTATCGGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGGGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                <  1:179210/148‑1 (MQ=255)
                                                                                                 cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                               >  1:1666500/1‑148 (MQ=255)
                                                                                                 cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                               >  1:2528326/1‑148 (MQ=255)
                                                                                                 cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACACGCTGAAGACCCGAAACGGGCATACGCCGGAACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAgag                                               >  1:1666484/1‑146 (MQ=255)
                                                                                                   cTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAg                                             <  2:1490785/148‑1 (MQ=255)
                                                                                                      ccGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCAt                                          <  2:2135795/148‑1 (MQ=255)
                                                                                                        ggTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGt                                        <  1:538347/148‑1 (MQ=255)
                                                                                                          ttATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCt                                      >  2:3427642/1‑148 (MQ=255)
                                                                                                          ttATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCt                                      >  2:614342/1‑148 (MQ=255)
                                                                                                          ttATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCt                                      >  1:2560323/1‑148 (MQ=255)
                                                                                                            atatGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGt                                    <  2:2814158/148‑1 (MQ=255)
                                                                                                            atatGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGt                                    <  1:1510317/148‑1 (MQ=255)
                                                                                                            atatGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGt                                    <  1:3272410/148‑1 (MQ=255)
                                                                                                             tatGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTa                                   >  2:2287689/1‑148 (MQ=255)
                                                                                                              atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTAc                                  <  2:2253499/148‑1 (MQ=255)
                                                                                                              atGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTAc                                  >  1:2543758/1‑148 (MQ=255)
                                                                                                               tGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGGGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACa                                 <  1:3546052/148‑1 (MQ=255)
                                                                                                                   aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                             >  1:2000769/1‑148 (MQ=255)
                                                                                                                     cTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                           <  1:3009210/148‑1 (MQ=255)
                                                                                                                     cTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                           >  2:381247/1‑148 (MQ=255)
                                                                                                                     cTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCAGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAAt                           <  1:1124828/148‑1 (MQ=255)
                                                                                                                       ggTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCg                         >  2:1545600/1‑148 (MQ=255)
                                                                                                                         tgtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGAc                       <  2:1765677/148‑1 (MQ=255)
                                                                                                                            ggTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                    >  2:707887/1‑148 (MQ=255)
                                                                                                                            ggTTATTCGTTAAACTACAGCTGGTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                    <  2:2604759/148‑1 (MQ=255)
                                                                                                                             gTTTTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTggg                   <  1:1184829/148‑1 (MQ=255)
                                                                                                                              ttATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGa                  <  2:92067/148‑1 (MQ=255)
                                                                                                                               tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                 >  2:395549/1‑148 (MQ=255)
                                                                                                                                aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                <  1:1462298/148‑1 (MQ=255)
                                                                                                                                  tCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTg              <  1:183821/148‑1 (MQ=255)
                                                                                                                                    gTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGc            >  2:3118230/1‑148 (MQ=255)
                                                                                                                                     ttAAACTCCAGCTGTGGTCGCGGCGACAAACGGGGCTGCGCAACCGGCTCATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGCGACGCCATCGGACGGAGTGAGCATGTCTGTACAGAAAATCGACTGGGATCTGGcc           <  1:688488/148‑1 (MQ=255)
                                                                                                                                     ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           <  2:2509642/148‑1 (MQ=255)
                                                                                                                                     ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           >  2:137404/1‑148 (MQ=255)
                                                                                                                                     atAAACTACAGCTGGTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc           <  2:618243/147‑1 (MQ=255)
                                                                                                                                      tAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc          >  1:3114856/1‑148 (MQ=255)
                                                                                                                                        aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGGACAGCAAATCGACTGGGATCTGGCCCTg        <  2:1906039/148‑1 (MQ=255)
                                                                                                                                        aaCTACAGCTGTTGTCGCGGCGAAAAATGTGGCTGCGCAACCGGATTATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCAATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg        <  1:871657/148‑1 (MQ=255)
                                                                                                                                          cTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt      >  1:1075653/1‑148 (MQ=255)
                                                                                                                                           tACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATc     >  2:1847449/1‑148 (MQ=255)
                                                                                                                                              aGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAg  >  2:1896098/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAAGAAGAGGAAGAAGACGAGAAGCAAGAAGACATGATGCGTCTGCTGCGGGGCAACTAACGGATTGCCGCCGTGGGCCTTCCCGTTCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAG  >  CP009273/4045157‑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: