Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,321,444 T→A 5.2% D60V (GAT→GTT)  yfaP ← DUF2135 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,321,4440TA 5.2% 1160.0 / 38.2 314D60V (GAT→GTT) yfaPDUF2135 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base T (158/139);  new base A (5/12);  total (163/151)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

AGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGCCAC  >  CP009273/2321296‑2321592
                                                                                                                                                    |                                                                                                                                                    
aGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGat                                                                                                                                                      <  1:332933/149‑1 (MQ=255)
aGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGat                                                                                                                                                      >  1:1190646/1‑149 (MQ=255)
aGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGat                                                                                                                                                      >  1:4619557/1‑149 (MQ=255)
 gTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatc                                                                                                                                                     <  1:3071433/149‑1 (MQ=255)
  tGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatca                                                                                                                                                    >  1:690197/1‑149 (MQ=255)
  tGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatca                                                                                                                                                    <  2:5084735/149‑1 (MQ=255)
     aaaTATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCt                                                                                                                                                 <  2:2057002/149‑1 (MQ=255)
     aaaTATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCt                                                                                                                                                 <  2:3702314/149‑1 (MQ=255)
     aaaTATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCt                                                                                                                                                 >  1:3800342/1‑149 (MQ=255)
       atatAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGc                                                                                                                                               >  1:2423323/1‑149 (MQ=255)
       atatAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGc                                                                                                                                               >  2:5201036/1‑149 (MQ=255)
        tataAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGcc                                                                                                                                              <  1:4241661/149‑1 (MQ=255)
        tataAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGcc                                                                                                                                              <  2:583967/149‑1 (MQ=255)
          taAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCAt                                                                                                                                            >  1:3238667/1‑149 (MQ=255)
          taAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCAt                                                                                                                                            <  2:4238578/149‑1 (MQ=255)
           aaGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATg                                                                                                                                           <  1:1872402/149‑1 (MQ=255)
            aGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGt                                                                                                                                          >  2:3075934/1‑149 (MQ=255)
            aGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGt                                                                                                                                          >  2:3122953/1‑149 (MQ=255)
             gggCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATACTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGtt                                                                                                                                         >  1:1403699/1‑149 (MQ=255)
               gCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTga                                                                                                                                       <  2:1069342/149‑1 (MQ=255)
                cGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgac                                                                                                                                      <  1:2229489/149‑1 (MQ=255)
                 gTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacg                                                                                                                                     >  1:176706/1‑149 (MQ=255)
                 gTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacg                                                                                                                                     <  1:4000816/149‑1 (MQ=255)
                  tGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacgg                                                                                                                                    >  1:2831134/1‑149 (MQ=255)
                   gCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacgga                                                                                                                                   <  2:3337210/149‑1 (MQ=255)
                      aaTGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGgacga                                                                                                                                >  2:443475/1‑149 (MQ=255)
                       aTGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGgacgac                                                                                                                               <  1:3525355/149‑1 (MQ=255)
                         gAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGc                                                                                                                             >  1:4815031/1‑149 (MQ=255)
                         gAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGc                                                                                                                             <  2:2848895/149‑1 (MQ=255)
                              aTCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGt                                                                                                                        <  1:5189458/149‑1 (MQ=255)
                              aTCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGt                                                                                                                        >  1:3518103/1‑149 (MQ=255)
                                cGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTaa                                                                                                                      <  2:1754289/149‑1 (MQ=255)
                                 ggATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAAt                                                                                                                     <  1:1244502/149‑1 (MQ=255)
                                 ggATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAAt                                                                                                                     >  2:605185/1‑149 (MQ=255)
                                        aTTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGa                                                                                                              <  2:379304/149‑1 (MQ=255)
                                        aTTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGa                                                                                                              >  2:646172/1‑149 (MQ=255)
                                         ttCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGAt                                                                                                             >  1:4921286/1‑149 (MQ=255)
                                          tctcTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGAtt                                                                                                            <  2:1286974/149‑1 (MQ=255)
                                                     tACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTg                                                                                                 >  1:832410/1‑149 (MQ=255)
                                                        tGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCa                                                                                              <  1:2295431/149‑1 (MQ=255)
                                                          gCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATc                                                                                            >  1:3379319/1‑149 (MQ=255)
                                                             ccGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTc                                                                                         >  2:1864977/1‑149 (MQ=255)
                                                             ccGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTc                                                                                         >  2:4784039/1‑149 (MQ=255)
                                                               gTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGc                                                                                       >  1:236415/1‑149 (MQ=255)
                                                               gTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGc                                                                                       <  1:3677603/149‑1 (MQ=255)
                                                               gTTGACCACCACCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGc                                                                                       >  1:499747/1‑149 (MQ=255)
                                                                      accaGCCGGCCCTGCTGAACTTTACCCGCCGCAGCAAAATACTTTATTTTGCCGAGGAACTGCGCTGAGATATTTTGAACATCTGCCATGTAGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGTCATCTTCTCCCTCTGc                                                                                >  1:4390958/1‑149 (MQ=255)
                                                                            cGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATg                                                                          >  1:1348395/1‑149 (MQ=255)
                                                                             ggCCCTGCTGAACTTTACCCGCCGCAGCCTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTGGACTGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGACCACTGCCGGATGc                                                                         >  1:1891641/1‑149 (MQ=255)
                                                                             ggCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGc                                                                         <  1:3105448/149‑1 (MQ=255)
                                                                               cccTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCa                                                                       <  1:1685968/149‑1 (MQ=255)
                                                                                  tgctgAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTACTTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGcc                                                                    >  2:4570122/1‑149 (MQ=255)
                                                                                    ctgAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTg                                                                  <  1:2124344/149‑1 (MQ=255)
                                                                                      gAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTgag                                                                >  1:4984936/1‑149 (MQ=255)
                                                                                       aaCTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTgaga                                                               >  2:2338028/1‑149 (MQ=255)
                                                                                         cTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGc                                                             >  1:3250495/1‑149 (MQ=255)
                                                                                          tttACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCg                                                            >  1:2850231/1‑149 (MQ=255)
                                                                                            tACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGt                                                          <  2:865142/149‑1 (MQ=255)
                                                                                             aCCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTg                                                         >  2:932371/1‑149 (MQ=255)
                                                                                                    gcagcaTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAAt                                                  >  2:2201476/1‑149 (MQ=255)
                                                                                                     cagcaTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAAtc                                                 >  1:3572944/1‑149 (MQ=255)
                                                                                                      agcaTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAAtct                                                <  2:1178955/149‑1 (MQ=255)
                                                                                                           aaTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACt                                           >  2:2390902/1‑149 (MQ=255)
                                                                                                           aaTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACt                                           <  2:2450826/149‑1 (MQ=255)
                                                                                                            aTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACtt                                          <  1:257640/149‑1 (MQ=255)
                                                                                                               cTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCt                                       >  1:2999770/1‑149 (MQ=255)
                                                                                                                tttATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTg                                      >  1:2807230/1‑149 (MQ=255)
                                                                                                                 ttATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGc                                     >  1:1656363/1‑149 (MQ=255)
                                                                                                                 ttATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGc                                     >  2:3188467/1‑149 (MQ=255)
                                                                                                                    ttttGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGcaca                                  >  2:1332470/1‑149 (MQ=255)
                                                                                                                       tGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCt                               <  1:2718277/149‑1 (MQ=255)
                                                                                                                         ccGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTc                             >  1:2534103/1‑149 (MQ=255)
                                                                                                                          cgcgGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCa                            >  1:4514473/1‑149 (MQ=255)
                                                                                                                           gcgGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCAc                           >  2:1852884/1‑149 (MQ=255)
                                                                                                                           gcgGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCAc                           <  1:102969/149‑1 (MQ=255)
                                                                                                                             ggATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtg                         >  1:3871011/1‑149 (MQ=255)
                                                                                                                              gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                        >  2:1523126/1‑149 (MQ=255)
                                                                                                                              gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                        <  2:3240182/149‑1 (MQ=255)
                                                                                                                              gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                        <  2:330961/149‑1 (MQ=255)
                                                                                                                              gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                        >  2:4128466/1‑149 (MQ=255)
                                                                                                                              gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                        >  1:962531/1‑149 (MQ=255)
                                                                                                                               aTCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgtg                       >  1:426050/1‑149 (MQ=255)
                                                                                                                               aTCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgtg                       <  2:5040729/149‑1 (MQ=255)
                                                                                                                               aTCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgtg                       <  2:4911730/149‑1 (MQ=255)
                                                                                                                                tCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGa                      >  1:2012208/1‑149 (MQ=255)
                                                                                                                                tCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGa                      >  1:1606845/1‑149 (MQ=255)
                                                                                                                                 cTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAg                     <  2:269387/149‑1 (MQ=255)
                                                                                                                                 cTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAg                     <  2:569741/149‑1 (MQ=255)
                                                                                                                                      cTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACa                <  2:744075/149‑1 (MQ=255)
                                                                                                                                       tGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCTGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAg               >  2:3958565/1‑149 (MQ=255)
                                                                                                                                        gagaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCTGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAgg              >  2:3413967/1‑149 (MQ=255)
                                                                                                                                        gagaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAgg              >  1:3757764/1‑149 (MQ=255)
                                                                                                                                        gagaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAgg              >  2:3513944/1‑149 (MQ=255)
                                                                                                                                        gagaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAgg              <  1:1107011/149‑1 (MQ=255)
                                                                                                                                         agaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGc             <  2:4745306/149‑1 (MQ=255)
                                                                                                                                            tatTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGaa          >  2:4490575/1‑149 (MQ=255)
                                                                                                                                              ttttGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAg        >  1:1302058/1‑149 (MQ=255)
                                                                                                                                                ttGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAgcg      <  1:443475/149‑1 (MQ=255)
                                                                                                                                                 tGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAgcgc     <  2:4063864/149‑1 (MQ=255)
                                                                                                                                                  gATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGcc    <  1:2953355/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGCCAC  >  CP009273/2321296‑2321592

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: ‑