Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,869,943 A→T 15.5% Y124N (TAC→AAC)  yidA ← sugar phosphate phosphatase; substrates include erythrose 4‑P and mannose 1‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,869,9430AT15.5% 492.4 / 72.7 179Y124N (TAC→AAC) yidAsugar phosphate phosphatase; substrates include erythrose 4‑P and mannose 1‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base A (82/69);  new base T (16/12);  total (98/81)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

CGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTG  >  CP009273/3869800‑3870090
                                                                                                                                               |                                                                                                                                                   
cGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTg                                                                                                                                                 <  1:927731/148‑1 (MQ=255)
       cgcgcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCAc                                                                                                                                          <  2:1035648/148‑1 (MQ=255)
        gcgcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACg                                                                                                                                         <  1:268904/148‑1 (MQ=255)
        gcgcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACg                                                                                                                                         <  1:268935/148‑1 (MQ=255)
        gcgcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACg                                                                                                                                         <  1:503515/148‑1 (MQ=255)
          gcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGt                                                                                                                                       <  2:1202132/148‑1 (MQ=255)
             aTAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTgg                                                                                                                                    >  2:1602928/1‑148 (MQ=255)
              tAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGc                                                                                                                                   >  2:2764007/1‑148 (MQ=255)
                 ccTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGgt                                                                                                                                <  1:2468980/148‑1 (MQ=255)
                 ccTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGgt                                                                                                                                <  1:2761388/148‑1 (MQ=255)
                 ccTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGgt                                                                                                                                >  1:3011993/1‑148 (MQ=255)
                   tGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGgtgt                                                                                                                              <  1:2767371/148‑1 (MQ=255)
                   tGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGgtgt                                                                                                                              >  2:1270690/1‑148 (MQ=255)
                    ggTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTa                                                                                                                             >  1:1801083/1‑148 (MQ=255)
                     gTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTAc                                                                                                                            >  2:3275130/1‑148 (MQ=255)
                     gTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTAc                                                                                                                            >  2:3275140/1‑148 (MQ=255)
                       cGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAg                                                                                                                          >  1:603579/1‑148 (MQ=255)
                       cGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAg                                                                                                                          <  1:1979500/148‑1 (MQ=255)
                          ggATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCgt                                                                                                                       <  1:1938894/148‑1 (MQ=255)
                           gATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCgtg                                                                                                                      >  1:1690958/1‑148 (MQ=255)
                            aTGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCgtgg                                                                                                                     >  1:2475315/1‑148 (MQ=255)
                                cGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCg                                                                                                                 <  2:2361577/148‑1 (MQ=255)
                                 gggTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCgg                                                                                                                <  1:1090339/148‑1 (MQ=255)
                                  ggTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGt                                                                                                               >  1:2440445/1‑148 (MQ=255)
                                   gTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTc                                                                                                              <  1:25758/148‑1 (MQ=255)
                                   gTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTc                                                                                                              <  2:419940/148‑1 (MQ=255)
                                     tcatcaATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCa                                                                                                            <  1:2362503/148‑1 (MQ=255)
                                       atcaatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAgg                                                                                                          >  1:1253201/1‑148 (MQ=255)
                                       atcaatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAgg                                                                                                          >  1:1253526/1‑148 (MQ=255)
                                       atcaatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAgg                                                                                                          >  1:1914286/1‑148 (MQ=255)
                                         caatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGc                                                                                                        >  2:745442/1‑148 (MQ=255)
                                            tcatcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTg                                                                                                     >  1:1905046/1‑148 (MQ=255)
                                             catcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTgg                                                                                                    >  2:2419469/1‑148 (MQ=255)
                                               tcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaa                                                                                                  >  1:1032444/1‑148 (MQ=255)
                                               tcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaa                                                                                                  >  1:1032434/1‑148 (MQ=255)
                                               tcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaa                                                                                                  <  2:3018596/148‑1 (MQ=255)
                                               tcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaa                                                                                                  >  2:92649/1‑148 (MQ=255)
                                                catcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaaa                                                                                                 >  1:2551575/1‑148 (MQ=255)
                                                catcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGaaa                                                                                                 >  1:1169488/1‑148 (MQ=255)
                                                 atcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAAt                                                                                                >  2:83343/1‑148 (MQ=255)
                                                  tcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATg                                                                                               >  2:1498860/1‑148 (MQ=255)
                                                    aCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgag                                                                                             <  2:3183542/148‑1 (MQ=255)
                                                    aCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgag                                                                                             <  1:776649/148‑1 (MQ=255)
                                                     cTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgaga                                                                                            >  2:2369894/1‑148 (MQ=255)
                                                      tttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGaa                                                                                           <  2:848836/148‑1 (MQ=255)
                                                      tttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGaa                                                                                           <  2:1039905/148‑1 (MQ=255)
                                                       ttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAAc                                                                                          >  2:2682703/1‑148 (MQ=255)
                                                       ttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAAc                                                                                          >  1:628598/1‑148 (MQ=255)
                                                       ttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGAGGAATGGTGTCAACGAAGGATACATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAAc                                                                                          >  1:1655253/1‑148 (MQ=255)
                                                         cAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCg                                                                                        <  2:772359/148‑1 (MQ=255)
                                                         cAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCg                                                                                        <  2:537864/148‑1 (MQ=255)
                                                          aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGa                                                                                       >  1:2006682/1‑148 (MQ=255)
                                                          aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGa                                                                                       >  1:263875/1‑148 (MQ=255)
                                                           ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGAc                                                                                      >  2:1066173/1‑148 (MQ=255)
                                                           ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGAc                                                                                      >  2:1314886/1‑148 (MQ=255)
                                                                 gggTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgcg                                                                                >  1:168787/1‑148 (MQ=255)
                                                                   gTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCgag                                                                              <  1:98394/148‑1 (MQ=255)
                                                                   gTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCgag                                                                              <  1:2419469/148‑1 (MQ=255)
                                                                           gTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGtttt                                                                      >  1:2433339/1‑148 (MQ=255)
                                                                             ccATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATGTCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTc                                                                    >  1:1952459/1‑148 (MQ=255)
                                                                              cATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTcc                                                                   >  2:336874/1‑148 (MQ=255)
                                                                                ttttCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAg                                                                 >  1:816744/1‑148 (MQ=255)
                                                                                 tttCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAgg                                                                >  2:1712013/1‑148 (MQ=255)
                                                                                    ctcCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGaaa                                                             >  1:1311591/1‑148 (MQ=255)
                                                                                    ctcCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGaaa                                                             <  2:2412038/148‑1 (MQ=255)
                                                                                     tcCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAAc                                                            >  2:531414/1‑148 (MQ=255)
                                                                                      ccGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACg                                                           <  1:1452144/148‑1 (MQ=255)
                                                                                       cGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGa                                                          <  1:2306411/148‑1 (MQ=255)
                                                                                        gCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGAt                                                         <  1:104662/148‑1 (MQ=255)
                                                                                        gCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGAt                                                         <  1:1983842/148‑1 (MQ=255)
                                                                                         cTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGaga                                                        >  2:1838150/1‑146 (MQ=255)
                                                                                         cTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATa                                                        >  2:2355376/1‑148 (MQ=255)
                                                                                         cTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATa                                                        >  2:2588589/1‑148 (MQ=255)
                                                                                         cTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGAAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAACTGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATa                                                        >  2:301533/1‑148 (MQ=255)
                                                                                           tCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAgt                                                      >  2:2287584/1‑148 (MQ=255)
                                                                                            cGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAgtc                                                     <  1:1192322/148‑1 (MQ=255)
                                                                                              cAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAgtcgt                                                   <  1:1517133/148‑1 (MQ=255)
                                                                                                  acacCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATa                                               <  2:221665/148‑1 (MQ=255)
                                                                                                    acCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTTTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGc                                             <  1:2958131/148‑1 (MQ=255)
                                                                                                       aGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTga                                          <  1:467016/148‑1 (MQ=255)
                                                                                                        gCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTgag                                         <  2:401482/148‑1 (MQ=255)
                                                                                                        gCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTgag                                         <  1:1168884/148‑1 (MQ=255)
                                                                                                             aTGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAg                                    >  1:836965/1‑148 (MQ=255)
                                                                                                              tggtggCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGt                                   <  1:2314868/148‑1 (MQ=255)
                                                                                                               ggtggCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGtt                                  <  1:2511665/148‑1 (MQ=255)
                                                                                                                   gCAACGAAGGATTCATGCACCGTGTAGTTGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTgcg                              >  1:1787249/1‑148 (MQ=255)
                                                                                                                    cAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTgcgc                             <  2:643034/148‑1 (MQ=255)
                                                                                                                    cAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTgcgc                             <  1:2839057/148‑1 (MQ=255)
                                                                                                                      aCGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCa                           >  2:284331/1‑148 (MQ=255)
                                                                                                                       cGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCAc                          >  1:219688/1‑148 (MQ=255)
                                                                                                                         aaGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACgg                        >  2:334575/1‑148 (MQ=255)
                                                                                                                          aGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTACAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGt                       >  2:3006936/1‑148 (MQ=255)
                                                                                                                           ggATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTg                      >  1:1995414/1‑148 (MQ=255)
                                                                                                                            gATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGc                     >  2:2130742/1‑148 (MQ=255)
                                                                                                                                 aTGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCa                <  1:3231120/148‑1 (MQ=255)
                                                                                                                                      ccGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGc           <  2:2004621/148‑1 (MQ=255)
                                                                                                                                        gtgtAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCgg         >  1:515163/1‑148 (MQ=255)
                                                                                                                                             gtagCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTc    <  2:3160074/148‑1 (MQ=255)
                                                                                                                                             gtagCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTc    <  1:1606472/148‑1 (MQ=255)
                                                                                                                                               agCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTg  >  2:574858/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAAACGATAGTCGTCATAGCTGAGAGCAGTTTGCGCCACGGTGCTACCATCAGCGGCCTTCTG  >  CP009273/3869800‑3870090

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: