Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,514,591 G→T 7.2% V79V (GTC→GTA yjhH ← putative lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,514,5910GT 7.2% 935.1 / 32.0 250V79V (GTC→GTAyjhHputative lyase/synthase
Reads supporting (aligned to +/- strand):  ref base G (124/108);  new base T (11/7);  total (135/115)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 254 total aligned reads displayed.

GTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGG  >  CP009273/4514442‑4514739
                                                                                                                                                     |                                                                                                                                                    
gTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCg                                                                                                                                                      >  1:1031223/1‑150 (MQ=255)
gTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCg                                                                                                                                                      >  2:4397314/1‑150 (MQ=255)
   aCGCTCCGGGCGATCTGCTGGTAATAGTCGTCAAGATGTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACg                                                                                                                                                   <  2:3184918/150‑1 (MQ=255)
         cGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCcaatc                                                                                                                                             <  2:2999663/150‑1 (MQ=255)
          gggCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCcaatca                                                                                                                                            <  1:1445001/150‑1 (MQ=255)
               aTCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATAcc                                                                                                                                       >  2:2043691/1‑150 (MQ=255)
                tCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCg                                                                                                                                      <  1:3890147/150‑1 (MQ=255)
                    ctgGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCAc                                                                                                                                  >  2:2436365/1‑150 (MQ=255)
                          aTAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGccc                                                                                                                            <  2:1594081/150‑1 (MQ=255)
                            aGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCTACGCCAATCAATACCGGCACTCGCCCGt                                                                                                                          >  1:2064965/1‑150 (MQ=255)
                            aGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGt                                                                                                                          >  1:218490/1‑150 (MQ=255)
                            aGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGt                                                                                                                          >  1:218479/1‑150 (MQ=255)
                             gtcgtcAAGATTTCGTTGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTc                                                                                                                         <  2:3868730/150‑1 (MQ=255)
                             gtcgtcAAGATTTCGTTGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTc                                                                                                                         <  2:1558736/150‑1 (MQ=255)
                              tcgtcAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCg                                                                                                                        <  2:2000487/150‑1 (MQ=255)
                               cgtcAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGTCGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGa                                                                                                                       <  2:4237998/150‑1 (MQ=255)
                                  cAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACaa                                                                                                                    <  2:4091245/150‑1 (MQ=255)
                                  cAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCATCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACaa                                                                                                                    >  2:2908926/1‑150 (MQ=255)
                                     gATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATgg                                                                                                                 >  1:1250158/1‑150 (MQ=255)
                                          cGTGGTGCGACTTTCCAGTTGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACa                                                                                                            <  1:1776036/150‑1 (MQ=255)
                                             ggtgCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGct                                                                                                         >  1:3562487/1‑150 (MQ=255)
                                               tgCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttc                                                                                                       >  1:85008/1‑150 (MQ=255)
                                               tgCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttc                                                                                                       >  1:1236418/1‑150 (MQ=255)
                                                gCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttct                                                                                                      <  1:4395205/150‑1 (MQ=255)
                                                 cGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCAGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCTACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttctt                                                                                                     >  1:4475444/1‑150 (MQ=255)
                                                  gACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTAGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCttcttc                                                                                                    >  2:2582814/1‑150 (MQ=255)
                                                  gACTTTCCAGTAGTAGGGGTAGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCttcttc                                                                                                    >  1:2706598/1‑150 (MQ=255)
                                                   aCTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCg                                                                                                   <  2:1153373/150‑1 (MQ=255)
                                                    cTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCgg                                                                                                  >  1:2523457/1‑150 (MQ=255)
                                                    cTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCgg                                                                                                  <  1:210407/150‑1 (MQ=255)
                                                     tttCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGc                                                                                                 >  2:3103964/1‑150 (MQ=255)
                                                         cAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGt                                                                                             >  2:2403743/1‑150 (MQ=255)
                                                            tagtagGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGcc                                                                                          >  1:3516493/1‑150 (MQ=255)
                                                               tagGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATg                                                                                       <  1:388338/150‑1 (MQ=255)
                                                                agGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATgc                                                                                      >  2:4328840/1‑150 (MQ=255)
                                                                 ggggTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAATGCCATgcg                                                                                     >  2:2954239/1‑150 (MQ=255)
                                                                   ggTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATAACGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCt                                                                                   >  2:591442/1‑150 (MQ=255)
                                                                    gTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTg                                                                                  <  1:4115758/150‑1 (MQ=255)
                                                                    gTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTg                                                                                  <  1:1800284/150‑1 (MQ=255)
                                                                        aTGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCt                                                                              <  1:1345142/150‑1 (MQ=255)
                                                                          ggCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGt                                                                            >  2:3753665/1‑150 (MQ=255)
                                                                          ggCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGt                                                                            >  2:3754237/1‑150 (MQ=255)
                                                                           gCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTa                                                                           <  1:2845594/150‑1 (MQ=255)
                                                                             gacgaTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTAtt                                                                         >  2:1732324/1‑150 (MQ=255)
                                                                             gacgaTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTAtt                                                                         <  1:295169/150‑1 (MQ=255)
                                                                                gaTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCAt                                                                      >  2:3872239/1‑150 (MQ=255)
                                                                                 aTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCAtt                                                                     >  1:2829693/1‑150 (MQ=255)
                                                                                 aTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCAtt                                                                     >  1:475560/1‑150 (MQ=255)
                                                                                  tACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCAttt                                                                    <  2:1373955/150‑1 (MQ=255)
                                                                                  tACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCAttt                                                                    <  1:272807/150‑1 (MQ=255)
                                                                                    ccATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTgg                                                                  >  2:1963224/1‑150 (MQ=255)
                                                                                     cATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACACTGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGc                                                                 >  2:3877050/1‑150 (MQ=255)
                                                                                       tCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTa                                                               <  1:3581094/150‑1 (MQ=255)
                                                                                        cAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTaa                                                              >  2:3286766/1‑150 (MQ=255)
                                                                                         aGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTaaa                                                             >  1:415194/1‑150 (MQ=255)
                                                                                          gcgcCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAAt                                                            <  1:2216901/150‑1 (MQ=255)
                                                                                           cgcCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAAtt                                                           <  1:4352528/150‑1 (MQ=255)
                                                                                            gcCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTc                                                          >  1:4021950/1‑150 (MQ=255)
                                                                                               gTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTcacc                                                       >  1:3074805/1‑150 (MQ=255)
                                                                                                tAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTcacca                                                      <  2:2335345/150‑1 (MQ=255)
                                                                                                tAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTcacca                                                      >  1:1106669/1‑150 (MQ=255)
                                                                                                  ggCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCaccacc                                                    <  1:1821079/150‑1 (MQ=255)
                                                                                                     ttGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGt                                                 <  1:725424/150‑1 (MQ=255)
                                                                                                     ttGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGt                                                 <  2:1031223/150‑1 (MQ=255)
                                                                                                     ttGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGt                                                 >  1:2987066/1‑150 (MQ=255)
                                                                                                      tGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTa                                                <  2:4528320/150‑1 (MQ=255)
                                                                                                        cgcATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTCTTCATTTGGCTAAATTCACCACCGGTAcc                                              >  2:3315629/1‑150 (MQ=255)
                                                                                                           aTGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAg                                           >  2:3940431/1‑150 (MQ=255)
                                                                                                            tgctgcGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGa                                          >  2:3377098/1‑150 (MQ=255)
                                                                                                            tgctgcGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGa                                          >  1:3195111/1‑150 (MQ=255)
                                                                                                                gcgcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATaaa                                      <  1:805249/150‑1 (MQ=255)
                                                                                                                 cgcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATaaaa                                     <  1:2582814/150‑1 (MQ=255)
                                                                                                                 cgcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATaaaa                                     >  1:816326/1‑150 (MQ=255)
                                                                                                                  gcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCTACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAAc                                    <  2:2061234/150‑1 (MQ=255)
                                                                                                                  gcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAAc                                    >  1:1640494/1‑150 (MQ=255)
                                                                                                                    cAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAg                                  <  1:948648/150‑1 (MQ=255)
                                                                                                                     aGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGc                                 >  1:4132504/1‑150 (MQ=255)
                                                                                                                       tttGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGccc                               <  1:2332090/150‑1 (MQ=255)
                                                                                                                          gACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTc                            >  2:480830/1‑150 (MQ=255)
                                                                                                                          gACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAAAGGTTACAGCTTCTTCGGCGAGAGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTc                            >  2:4453069/1‑150 (MQ=255)
                                                                                                                            ccGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGa                          <  1:1952399/150‑1 (MQ=255)
                                                                                                                                   gTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCttt                   >  2:1471077/1‑150 (MQ=255)
                                                                                                                                     cAGTGGAAGGGGAACCTACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTttat                 >  2:3604664/1‑150 (MQ=255)
                                                                                                                                     cAGTGGAAGGGGAACCTACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTttat                 >  2:3604715/1‑150 (MQ=255)
                                                                                                                                     cAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTttat                 <  1:292035/150‑1 (MQ=255)
                                                                                                                                     cAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTttat                 <  1:858140/150‑1 (MQ=255)
                                                                                                                                       gTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTttatta               >  2:3781141/1‑150 (MQ=255)
                                                                                                                                          gAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATc            >  2:785803/1‑150 (MQ=255)
                                                                                                                                           aaGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCa           >  1:1777011/1‑150 (MQ=255)
                                                                                                                                           aaGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCa           <  1:983566/150‑1 (MQ=255)
                                                                                                                                           aaGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCa           >  1:1518451/1‑150 (MQ=255)
                                                                                                                                              gggAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGa        >  1:3452275/1‑150 (MQ=255)
                                                                                                                                              gggAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGa        <  1:3872192/150‑1 (MQ=255)
                                                                                                                                                gAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCTTCGTCCCCTTTATTAATCAGGAAg      <  1:4447601/150‑1 (MQ=255)
                                                                                                                                                gAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAg      <  1:2043691/150‑1 (MQ=255)
                                                                                                                                                gAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAg      <  1:2043536/150‑1 (MQ=255)
                                                                                                                                                gAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAg      <  1:3789805/150‑1 (MQ=255)
                                                                                                                                                 aaCCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGt     <  1:968844/150‑1 (MQ=255)
                                                                                                                                                  aCCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTc    <  2:3251721/150‑1 (MQ=255)
                                                                                                                                                    cGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCgg  >  1:339623/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGG  >  CP009273/4514442‑4514739

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: