Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,882,054 C→A 10.6% D469Y (GAT→TAT)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,882,0540CA10.6% 580.8 / 32.4 170D469Y (GAT→TAT) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base C (75/77);  new base A (12/6);  total (87/83)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 171 total aligned reads displayed.

ACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTA  >  CP009273/2881907‑2882200
                                                                                                                                                   |                                                                                                                                                  
aCATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATc                                                                                                                                                    >  1:486851/1‑148 (MQ=255)
 cATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                                                                                                                                                   >  2:323145/1‑148 (MQ=255)
 cATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                                                                                                                                                   >  2:323159/1‑148 (MQ=255)
   tCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCt                                                                                                                                                 >  2:645258/1‑148 (MQ=255)
    ccGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTg                                                                                                                                                >  2:1034689/1‑148 (MQ=255)
      gTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac                                                                                                                                              >  1:1378213/1‑148 (MQ=255)
      gTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGac                                                                                                                                              >  1:3018601/1‑148 (MQ=255)
        gggATACGTGTCCAAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGTCTGTTACACGCATCACGATATGCTCATCGCTGacac                                                                                                                                            <  1:2736456/148‑1 (MQ=255)
         ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACAcc                                                                                                                                           <  1:229201/148‑1 (MQ=255)
         ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACAcc                                                                                                                                           <  1:2947027/148‑1 (MQ=255)
          gATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCa                                                                                                                                          >  1:1738729/1‑148 (MQ=255)
          gATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCa                                                                                                                                          <  2:266190/148‑1 (MQ=255)
            tACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg                                                                                                                                        <  2:1600162/148‑1 (MQ=255)
            tACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATg                                                                                                                                        >  1:505447/1‑148 (MQ=255)
             aCGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGt                                                                                                                                       >  1:1899221/1‑148 (MQ=255)
              cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGtt                                                                                                                                      <  1:106486/148‑1 (MQ=255)
               gtgtCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGttt                                                                                                                                     <  2:1768792/148‑1 (MQ=255)
                tgtCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTc                                                                                                                                    >  2:1183778/1‑148 (MQ=255)
                tgtCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTc                                                                                                                                    >  1:1590640/1‑148 (MQ=255)
                   cccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGcc                                                                                                                                 >  2:1493497/1‑148 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               <  1:74674/148‑1 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               <  1:1608748/148‑1 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               <  1:2265091/148‑1 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               <  1:598700/148‑1 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               <  1:74781/148‑1 (MQ=255)
                     cAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCAt                                                                                                                               >  1:329785/1‑148 (MQ=255)
                        tGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTaa                                                                                                                            >  2:2794917/1‑148 (MQ=255)
                         gCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTaaa                                                                                                                           <  2:450219/148‑1 (MQ=255)
                         gCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTaaa                                                                                                                           <  1:35700/148‑1 (MQ=255)
                          cGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAAt                                                                                                                          <  2:2645146/148‑1 (MQ=255)
                           ggTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAAtt                                                                                                                         <  2:1514429/148‑1 (MQ=255)
                             ttGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTAt                                                                                                                       <  2:1828572/148‑1 (MQ=255)
                                 cgccAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAt                                                                                                                   >  2:1559907/1‑148 (MQ=255)
                                    cAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtgtt                                                                                                                <  1:993030/148‑1 (MQ=255)
                                    cAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAtgtt                                                                                                                >  2:2958957/1‑148 (MQ=255)
                                       gATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTc                                                                                                             >  1:1020426/1‑148 (MQ=255)
                                        aTGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCg                                                                                                            >  1:1646930/1‑148 (MQ=255)
                                             ggTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATaaa                                                                                                       <  2:2302581/148‑1 (MQ=255)
                                              gTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATaaaa                                                                                                      <  1:3290326/148‑1 (MQ=255)
                                                tGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGa                                                                                                    >  1:1948769/1‑148 (MQ=255)
                                                   aGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACgg                                                                                                 <  2:1566488/148‑1 (MQ=255)
                                                   aGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACgg                                                                                                 <  2:3291719/148‑1 (MQ=255)
                                                    gCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGc                                                                                                >  2:499619/1‑148 (MQ=255)
                                                     cGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATAGATGTTGTCGATAAAAGACGGCa                                                                                               >  1:655557/1‑148 (MQ=255)
                                                       aCCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAgg                                                                                             <  1:2445746/148‑1 (MQ=255)
                                                         ccGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGaa                                                                                           <  2:2676095/148‑1 (MQ=255)
                                                              gcTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGct                                                                                      <  1:2285567/148‑1 (MQ=255)
                                                               cTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                                                     <  2:2739987/148‑1 (MQ=255)
                                                               cTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGctc                                                                                     <  1:1898339/148‑1 (MQ=255)
                                                                         caGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGcc                                                                           >  1:834005/1‑148 (MQ=255)
                                                                           ggCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCAt                                                                         <  1:3082969/148‑1 (MQ=255)
                                                                                aCTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCa                                                                    >  1:2583136/1‑148 (MQ=255)
                                                                                aCTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCa                                                                    >  1:929124/1‑148 (MQ=255)
                                                                                        cAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCaa                                                            <  1:2197316/148‑1 (MQ=255)
                                                                                          gCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGt                                                          >  2:2241034/1‑148 (MQ=255)
                                                                                             tCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCgg                                                       >  1:2360790/1‑148 (MQ=255)
                                                                                              cgcgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCACGTCggg                                                      >  2:1438256/1‑148 (MQ=255)
                                                                                              cgcgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                                                      >  2:730383/1‑148 (MQ=255)
                                                                                              cgcgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCggg                                                      >  2:730309/1‑148 (MQ=255)
                                                                                                cgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGaa                                                    <  1:26020/148‑1 (MQ=255)
                                                                                                cgcgACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGaa                                                    <  2:3120510/148‑1 (MQ=255)
                                                                                                   gACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGa                                                 <  2:1009359/148‑1 (MQ=255)
                                                                                                        caACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGc                                            >  1:2755740/1‑148 (MQ=255)
                                                                                                         aaCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCTCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa                                           <  2:667784/148‑1 (MQ=255)
                                                                                                         aaCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa                                           >  2:2342107/1‑148 (MQ=255)
                                                                                                         aaCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCa                                           <  2:667906/148‑1 (MQ=255)
                                                                                                           ccGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAg                                         >  2:2339971/1‑148 (MQ=255)
                                                                                                           ccGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAg                                         <  1:668262/148‑1 (MQ=255)
                                                                                                              ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTTTCGATAAAATACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                      >  1:929096/1‑148 (MQ=255)
                                                                                                              ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                      <  2:1436724/148‑1 (MQ=255)
                                                                                                              ttCGGGCAGCCTGTTACACGCATCACGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCa                                      >  1:914004/1‑148 (MQ=255)
                                                                                                                 gggCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCg                                   <  2:1755796/148‑1 (MQ=255)
                                                                                                                  ggCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGa                                  >  1:2525883/1‑148 (MQ=255)
                                                                                                                  ggCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGa                                  >  1:2525950/1‑148 (MQ=255)
                                                                                                                      gCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGttt                              >  2:3074008/1‑148 (MQ=255)
                                                                                                                      gCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGttt                              >  1:1788715/1‑148 (MQ=255)
                                                                                                                       ccTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtttt                             >  2:3072444/1‑148 (MQ=255)
                                                                                                                       ccTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtttt                             >  1:2825450/1‑148 (MQ=255)
                                                                                                                       ccTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtttt                             >  2:1690939/1‑148 (MQ=255)
                                                                                                                       ccTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGtttt                             >  2:1690892/1‑148 (MQ=255)
                                                                                                                          gTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCAc                          >  1:2929941/1‑148 (MQ=255)
                                                                                                                            tACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGc                        >  1:456179/1‑148 (MQ=255)
                                                                                                                                    tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                >  1:3189253/1‑148 (MQ=255)
                                                                                                                                    tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                >  1:1548318/1‑148 (MQ=255)
                                                                                                                                    tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                <  2:1590640/148‑1 (MQ=255)
                                                                                                                                     cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg               <  1:1493497/148‑1 (MQ=255)
                                                                                                                                     cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg               >  2:2839207/1‑148 (MQ=255)
                                                                                                                                      aCGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGa              >  1:1139452/1‑148 (MQ=255)
                                                                                                                                      aCGATATGCTCATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGa              >  1:1139451/1‑148 (MQ=255)
                                                                                                                                        gATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACg            >  1:1558821/1‑148 (MQ=255)
                                                                                                                                         atatGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACgg           >  2:1333401/1‑148 (MQ=255)
                                                                                                                                           atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca         <  1:628417/148‑1 (MQ=255)
                                                                                                                                           atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGca         <  1:138860/148‑1 (MQ=255)
                                                                                                                                             gCTCATCGCTGACACCATGTTTCACCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt       >  1:2210545/1‑148 (MQ=255)
                                                                                                                                              cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc      <  1:2653692/148‑1 (MQ=255)
                                                                                                                                              cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc      <  1:2653703/148‑1 (MQ=255)
                                                                                                                                               tCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACTCGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattca     <  1:1034689/148‑1 (MQ=255)
                                                                                                                                                cATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcat    >  1:1568640/1‑148 (MQ=255)
                                                                                                                                                cATAGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcat    <  1:1238701/148‑1 (MQ=255)
                                                                                                                                                  tCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa  <  2:550929/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTA  >  CP009273/2881907‑2882200

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: