Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,654,099 C→A 19.3% A264A (GCG→GCT iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,654,0990CA19.3% 794.1 / 148.1 275A264A (GCG→GCTiscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (131/91);  new base A (27/26);  total (158/117)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 280 total aligned reads displayed.

TCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACA  >  CP009273/2653951‑2654244
                                                                                                                                                    |                                                                                                                                                 
tCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTc                                                                                                                                                   >  1:2650003/1‑149 (MQ=255)
tCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTc                                                                                                                                                   >  1:1232717/1‑149 (MQ=255)
 cAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCg                                                                                                                                                  >  1:1481283/1‑149 (MQ=255)
    cGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCa                                                                                                                                               >  2:925071/1‑149 (MQ=255)
       cTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtct                                                                                                                                            >  1:1726524/1‑149 (MQ=255)
        tCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc                                                                                                                                           <  1:4192349/149‑1 (MQ=255)
         cACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATctct                                                                                                                                          >  2:386015/1‑149 (MQ=255)
          aCCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt                                                                                                                                         >  2:1769743/1‑149 (MQ=255)
          aCCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtctt                                                                                                                                         >  2:501093/1‑149 (MQ=255)
           ccTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCtcttc                                                                                                                                        >  2:1141109/1‑149 (MQ=255)
           ccTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCtcttc                                                                                                                                        >  1:4435571/1‑149 (MQ=255)
           ccTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCtcttc                                                                                                                                        >  1:1801102/1‑149 (MQ=255)
             ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt                                                                                                                                      >  1:622813/1‑149 (MQ=255)
             ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt                                                                                                                                      >  1:4141830/1‑149 (MQ=255)
             ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTcttctt                                                                                                                                      <  1:1805392/149‑1 (MQ=255)
             ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTcttctt                                                                                                                                      >  1:4068229/1‑149 (MQ=255)
              tCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCttt                                                                                                                                     >  1:1672568/1‑149 (MQ=255)
              tCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCttt                                                                                                                                     >  2:1865044/1‑149 (MQ=255)
              tCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCttt                                                                                                                                     <  1:2708124/149‑1 (MQ=255)
              tCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCttt                                                                                                                                     <  1:366512/149‑1 (MQ=255)
                aaCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtg                                                                                                                                   >  1:2424882/1‑149 (MQ=255)
                 aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc                                                                                                                                  <  2:2293432/149‑1 (MQ=255)
                 aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc                                                                                                                                  <  2:3439677/149‑1 (MQ=255)
                 aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc                                                                                                                                  <  2:1398183/149‑1 (MQ=255)
                   gTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcga                                                                                                                                >  1:4475747/1‑149 (MQ=255)
                       ttGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTtgcgatgcg                                                                                                                            <  2:1374440/149‑1 (MQ=255)
                           aGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAg                                                                                                                        >  1:1730067/1‑149 (MQ=255)
                            gCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAgg                                                                                                                       >  1:3097555/1‑149 (MQ=255)
                               gACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCt                                                                                                                    <  1:870708/149‑1 (MQ=255)
                                aCGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTc                                                                                                                   >  1:3236644/1‑149 (MQ=255)
                                      agaATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCAt                                                                                                             <  1:1417648/149‑1 (MQ=255)
                                       gaATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCAtt                                                                                                            >  1:672427/1‑149 (MQ=255)
                                       gaATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCAtt                                                                                                            >  2:3006447/1‑149 (MQ=255)
                                        aaTGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTc                                                                                                           <  2:989840/149‑1 (MQ=255)
                                         aTGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTcc                                                                                                          <  2:1205556/149‑1 (MQ=255)
                                         aTGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTcc                                                                                                          <  1:3391762/149‑1 (MQ=255)
                                           gTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCcga                                                                                                        <  2:288151/149‑1 (MQ=255)
                                              cGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCcgacga                                                                                                     <  1:2959368/149‑1 (MQ=255)
                                                 cgcACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCt                                                                                                  >  1:353934/1‑149 (MQ=255)
                                                    aCCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGt                                                                                               >  1:2199098/1‑149 (MQ=255)
                                                    aCCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGt                                                                                               >  1:2199274/1‑149 (MQ=255)
                                                     ccGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTg                                                                                              >  2:4140881/1‑149 (MQ=255)
                                                     ccGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTg                                                                                              >  1:2758913/1‑149 (MQ=255)
                                                      cGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACTATCTGGTGa                                                                                             >  1:3526429/1‑149 (MQ=255)
                                                       gtgtTCCAGGTCACCTTTCAGGTAAACTACTACGATATCTTTGATGCCGTTCAACAGACGGATAAGCAGGCCGCTCAGACGTTCCATCTCTGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTACGACTATCTGGTGaa                                                                                            >  2:237118/1‑149 (MQ=255)
                                                        tgtTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAAc                                                                                           <  2:3689263/149‑1 (MQ=255)
                                                         gtTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACa                                                                                          >  2:806983/1‑149 (MQ=255)
                                                           tCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACTATCTGGTGAACAgg                                                                                        <  2:3486189/149‑1 (MQ=255)
                                                            ccAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACTATCTGGTGAACAGGc                                                                                       >  2:4496359/1‑149 (MQ=255)
                                                            ccAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGc                                                                                       <  1:4228483/149‑1 (MQ=255)
                                                             cAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCa                                                                                      <  1:1266524/149‑1 (MQ=255)
                                                               ggTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCaga                                                                                    >  1:277512/1‑149 (MQ=255)
                                                                 tCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGt                                                                                  <  1:3027631/149‑1 (MQ=255)
                                                                  cACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTg                                                                                 >  1:746121/1‑149 (MQ=255)
                                                                   aCCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGc                                                                                <  2:4458563/149‑1 (MQ=255)
                                                                   aCCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGc                                                                                >  1:3251743/1‑149 (MQ=255)
                                                                       ttCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGa                                                                            >  1:1117476/1‑149 (MQ=255)
                                                                          aGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACg                                                                         >  1:4518971/1‑149 (MQ=255)
                                                                           ggTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGc                                                                        <  2:857641/149‑1 (MQ=255)
                                                                             tAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCAt                                                                      <  2:1312963/149‑1 (MQ=255)
                                                                             tAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCAt                                                                      <  2:3668946/149‑1 (MQ=255)
                                                                             tAAACTTCTTCGATATCTTTGATGCCGGTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCAt                                                                      <  2:265319/149‑1 (MQ=255)
                                                                                   tcttcGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACcgcg                                                                >  2:1571308/1‑149 (MQ=255)
                                                                                      tcGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTc                                                             <  2:3409859/149‑1 (MQ=255)
                                                                                      tcGATATCTTTGATGCCGGTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTc                                                             <  1:1944949/149‑1 (MQ=255)
                                                                                          tatCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGa                                                         >  2:1710001/1‑149 (MQ=255)
                                                                                          tatCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGa                                                         <  2:4082258/149‑1 (MQ=255)
                                                                                          gatCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGa                                                         <  2:2474740/148‑1 (MQ=255)
                                                                                           atCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAc                                                        <  2:1952098/149‑1 (MQ=255)
                                                                                            tCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAcc                                                       >  1:3449065/1‑149 (MQ=255)
                                                                                              tttGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACTATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAccgc                                                     >  2:437546/1‑149 (MQ=255)
                                                                                                  aTGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAccgccgcc                                                 >  2:203619/1‑149 (MQ=255)
                                                                                                   tGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAccgccgccg                                                >  2:3697844/1‑149 (MQ=255)
                                                                                                      cGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGc                                             <  2:778901/149‑1 (MQ=255)
                                                                                                       gTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCa                                            >  2:2258829/1‑149 (MQ=255)
                                                                                                        ttCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCAt                                           >  1:1315334/1‑149 (MQ=255)
                                                                                                         tCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCAtt                                          <  1:268062/149‑1 (MQ=255)
                                                                                                         tCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCAtt                                          <  2:4091379/149‑1 (MQ=255)
                                                                                                           cacaGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTg                                        >  1:3744475/1‑149 (MQ=255)
                                                                                                              aGACGGGTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTgcgc                                     <  2:3069982/149‑1 (MQ=255)
                                                                                                                  ggTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCg                                 <  1:1533886/149‑1 (MQ=255)
                                                                                                                   gTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGa                                <  2:707347/149‑1 (MQ=255)
                                                                                                                      aCGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGc                             >  1:2468142/1‑149 (MQ=255)
                                                                                                                          aGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTAc                         <  1:822291/149‑1 (MQ=255)
                                                                                                                          aGGCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTACATTTGCGCTTCGATGCGTAc                         >  1:4161934/1‑149 (MQ=255)
                                                                                                                           ggCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTAcg                        >  1:1097589/1‑149 (MQ=255)
                                                                                                                            gCCGCGCAGACGTTCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTAcgc                       <  1:4347602/149‑1 (MQ=255)
                                                                                                                                  cAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGtt                 >  2:1991249/1‑149 (MQ=255)
                                                                                                                                   aGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGttt                >  1:4613554/1‑149 (MQ=255)
                                                                                                                                    gACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTa               >  1:602203/1‑149 (MQ=255)
                                                                                                                                       gTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTacga            >  2:4285538/1‑149 (MQ=255)
                                                                                                                                         tCCATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTacgacg          >  1:1871253/1‑149 (MQ=255)
                                                                                                                                          ccATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCCTTTGCGCTTCGATGCGTCCGCGCGGTTTACGACGt         >  2:2899808/1‑149 (MQ=255)
                                                                                                                                           cATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGCCGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGta        <  2:3953610/149‑1 (MQ=255)
                                                                                                                                           cATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACTATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGta        <  2:3150534/149‑1 (MQ=255)
                                                                                                                                           cATCTCGGTAGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGta        <  2:34080/149‑1 (MQ=255)
                                                                                                                                             tctcGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGtaca      >  2:4425993/1‑149 (MQ=255)
                                                                                                                                             tctcGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGtaca      >  2:1613501/1‑149 (MQ=255)
                                                                                                                                                 ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGTACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGtacataca  <  2:517495/149‑1 (MQ=255)
                                                                                                                                                 ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGtacataca  <  1:501093/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACA  >  CP009273/2653951‑2654244

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: