Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,753,424 G→T 57.8% G361* (GGA→TGA)  yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,753,4240GT57.8% 26.2 / 300.0 185G361* (GGA→TGA) yfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (30/48);  new base T (37/70);  total (67/118)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 185 total aligned reads displayed.

TTCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATT  >  CP009273/2753298‑2753571
                                                                                                                              |                                                                                                                                                   
ttCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCa                                                                                                                              <  1:2698108/150‑1 (MQ=31)
ttCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCa                                                                                                                              <  1:2698017/150‑1 (MQ=31)
  cTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGt                                                                                                                            <  2:53726/150‑1 (MQ=32)
  cTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGt                                                                                                                            <  2:1690359/150‑1 (MQ=32)
   tACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTa                                                                                                                           >  1:527903/1‑150 (MQ=18)
          aGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAagcag                                                                                                                    <  1:2114515/150‑1 (MQ=255)
          aGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAagcag                                                                                                                    >  2:28992/1‑150 (MQ=255)
            gTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAgcagca                                                                                                                  >  1:1828409/1‑150 (MQ=255)
              cGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGt                                                                                                                >  2:1713568/1‑150 (MQ=255)
              cGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGt                                                                                                                <  2:1957987/150‑1 (MQ=255)
              cGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGt                                                                                                                <  2:1958002/150‑1 (MQ=255)
               gAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTg                                                                                                               <  1:2279527/150‑1 (MQ=255)
               gAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTg                                                                                                               >  2:1841353/1‑150 (MQ=255)
                aGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTgg                                                                                                              <  1:4255098/150‑1 (MQ=255)
                 gTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGga                                                                                                             >  1:560284/1‑150 (MQ=255)
                 gTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGga                                                                                                             >  1:2336486/1‑150 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAATGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  2:3557062/150‑1 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  2:3557030/150‑1 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  1:1140972/150‑1 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  1:3131749/150‑1 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  2:3703963/150‑1 (MQ=255)
                        aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  1:1140950/150‑1 (MQ=255)
                                 aTTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGtat                                                                                             <  2:4471221/150‑1 (MQ=255)
                                      ccTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTga                                                                                        <  1:1568714/150‑1 (MQ=255)
                                      ccTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTga                                                                                        >  2:2440767/1‑150 (MQ=255)
                                       cttcttTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaa                                                                                       <  1:2264245/150‑1 (MQ=255)
                                        ttcttTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaag                                                                                      <  2:3892990/150‑1 (MQ=255)
                                        ttcttTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaag                                                                                      <  2:2491067/150‑1 (MQ=255)
                                          cttTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaagaa                                                                                    >  1:932031/1‑150 (MQ=255)
                                            ttATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATa                                                                                  <  2:2855062/150‑1 (MQ=255)
                                            ttATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATa                                                                                  <  2:410680/150‑1 (MQ=255)
                                             tATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATAc                                                                                 >  1:920188/1‑150 (MQ=255)
                                               tGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACAt                                                                               >  1:4061010/1‑150 (MQ=255)
                                               tGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACAt                                                                               >  1:2429587/1‑150 (MQ=255)
                                                gATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATc                                                                              <  2:258101/150‑1 (MQ=255)
                                                 aTTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGATGCTTGCTGCCTTACTTCATTACTTCAGAGGTGATGGAGGCGATATATCTGATACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACAAcc                                                                             >  1:1382907/1‑150 (MQ=255)
                                                 aTTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATcc                                                                             >  1:3259677/1‑150 (MQ=255)
                                                  ttttAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCg                                                                            <  2:2523212/150‑1 (MQ=255)
                                                    ttAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCt                                                                          >  2:4403395/1‑150 (MQ=255)
                                                    ttAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCt                                                                          >  2:4403422/1‑150 (MQ=255)
                                                    ttAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCt                                                                          >  2:2073520/1‑150 (MQ=255)
                                                     tAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCtt                                                                         >  2:3947004/1‑150 (MQ=255)
                                                     tAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCtt                                                                         <  1:1813423/150‑1 (MQ=255)
                                                     tAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCtt                                                                         >  2:1618062/1‑150 (MQ=255)
                                                     tAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTAATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCtt                                                                         <  2:1560602/150‑1 (MQ=255)
                                                       aGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGt                                                                       >  2:1870682/1‑150 (MQ=255)
                                                        ggATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGtt                                                                      <  1:3544233/150‑1 (MQ=255)
                                                        ggATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGtt                                                                      <  2:1343826/150‑1 (MQ=255)
                                                         gATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTtc                                                                     >  1:2319858/1‑150 (MQ=255)
                                                           ttATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTtctc                                                                   >  2:3317352/1‑150 (MQ=255)
                                                           ttATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTtctc                                                                   <  1:4544911/150‑1 (MQ=255)
                                                             aTGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTa                                                                 >  2:2542318/1‑150 (MQ=255)
                                                             aTGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTa                                                                 >  1:2544349/1‑150 (MQ=255)
                                                              tGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaa                                                                <  2:4104879/150‑1 (MQ=255)
                                                                    tAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaaga                                                          <  2:4606499/150‑1 (MQ=255)
                                                                    tAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaaga                                                          >  2:4130221/1‑150 (MQ=255)
                                                                      aaaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaagaag                                                        <  1:3804767/150‑1 (MQ=255)
                                                                      aaaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaagaag                                                        >  1:1212632/1‑150 (MQ=255)
                                                                       aaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAagaaga                                                       <  1:4236201/150‑1 (MQ=255)
                                                                         tGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGt                                                     <  2:1739782/150‑1 (MQ=255)
                                                                         tGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGt                                                     <  2:2418702/150‑1 (MQ=255)
                                                                          ggCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGtt                                                    >  1:393154/1‑150 (MQ=255)
                                                                           gCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGttt                                                   <  1:150217/150‑1 (MQ=255)
                                                                             ggAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTc                                                 >  1:3236005/1‑150 (MQ=255)
                                                                               agaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCtt                                               <  1:4154340/150‑1 (MQ=255)
                                                                                gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              <  1:2079253/150‑1 (MQ=255)
                                                                                gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              <  1:519692/150‑1 (MQ=255)
                                                                                gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              <  2:16419/150‑1 (MQ=255)
                                                                                gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              >  2:3914785/1‑150 (MQ=255)
                                                                                gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              >  2:2136563/1‑150 (MQ=255)
                                                                                  aCATGGCGAGGCTTGCTGCCTTACTTCTTTACTTCAGCGTTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCATCAGTGGTGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAg                                            <  1:4618043/150‑1 (MQ=37)
                                                                                  aCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAg                                            >  1:4169322/1‑150 (MQ=255)
                                                                                    aTGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGAt                                          <  1:477611/150‑1 (MQ=255)
                                                                                    aTGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGAt                                          >  1:739590/1‑150 (MQ=255)
                                                                                    aTGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGAt                                          <  2:1069030/150‑1 (MQ=255)
                                                                                      ggCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGt                                        <  2:3014289/150‑1 (MQ=255)
                                                                                      ggCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGt                                        >  2:2283779/1‑150 (MQ=255)
                                                                                         gAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTc                                     >  1:3074911/1‑150 (MQ=255)
                                                                                            gCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGa                                  <  1:639046/150‑1 (MQ=255)
                                                                                             cTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGaa                                 >  2:2277211/1‑150 (MQ=255)
                                                                                             cTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGaa                                 <  1:4528731/150‑1 (MQ=255)
                                                                                              ttGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                <  1:4115199/150‑1 (MQ=255)
                                                                                              ttGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                >  1:3578553/1‑150 (MQ=255)
                                                                                               tgctgcCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGc                               <  2:2192961/150‑1 (MQ=255)
                                                                                               tgctgcCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGc                               <  2:3856497/150‑1 (MQ=255)
                                                                                                gctgcCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCa                              >  2:942176/1‑150 (MQ=255)
                                                                                                 ctgcCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAg                             <  2:2346363/150‑1 (MQ=255)
                                                                                                       tacttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGc                       >  2:1443109/1‑150 (MQ=255)
                                                                                                           tcattacttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTa                   <  2:4284157/150‑1 (MQ=255)
                                                                                                             attacttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAtt                 >  1:258662/1‑150 (MQ=255)
                                                                                                             attacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAtt                 <  2:3339898/150‑1 (MQ=255)
                                                                                                               tacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgt               <  1:3092826/150‑1 (MQ=255)
                                                                                                                acttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgtt              >  1:4185497/1‑150 (MQ=255)
                                                                                                                 cttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgttg             >  2:4221292/1‑150 (MQ=255)
                                                                                                                  ttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTgg            <  1:60731/150‑1 (MQ=255)
                                                                                                                  ttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTgg            >  2:1312560/1‑150 (MQ=255)
                                                                                                                     aGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATa         <  2:3823753/150‑1 (MQ=255)
                                                                                                                         gTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAg     >  1:1965867/1‑150 (MQ=255)
                                                                                                                           gATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGAt   <  2:560284/150‑1 (MQ=255)
                                                                                                                            aTGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGAtt  >  2:216583/1‑150 (MQ=255)
                                                                                                                              |                                                                                                                                                   
TTCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATT  >  CP009273/2753298‑2753571

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: