Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,753,424 G→T 19.6% 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,4240GT19.6% 677.0 / 116.1 240G361* (GGA→TGA) yfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (89/104);  new base T (24/23);  total (113/127)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 240 total aligned reads displayed.

ATTCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATT  >  CP009273/2753297‑2753571
                                                                                                                               |                                                                                                                                                   
 ttCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCa                                                                                                                              <  2:201179/150‑1 (MQ=14)
   cTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGt                                                                                                                            >  1:4157116/1‑150 (MQ=32)
   cTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGt                                                                                                                            >  2:2804166/1‑150 (MQ=32)
      cAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTaaa                                                                                                                         >  2:4374996/1‑150 (MQ=35)
        aCCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGTTTATGCTTCTAAAATGGCGGAGATCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAag                                                                                                                       <  2:3182701/150‑1 (MQ=255)
        aCCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAag                                                                                                                       <  2:1171302/150‑1 (MQ=37)
        aCCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAag                                                                                                                       <  2:3477177/150‑1 (MQ=37)
         ccAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAagc                                                                                                                      <  1:267740/150‑1 (MQ=38)
         ccAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAagc                                                                                                                      >  2:1482917/1‑150 (MQ=38)
         ccAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAagc                                                                                                                      >  2:3759667/1‑150 (MQ=38)
              tCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGcagcag                                                                                                                 <  2:1884344/150‑1 (MQ=255)
                gAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTg                                                                                                               >  1:1238920/1‑150 (MQ=255)
                  gTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGga                                                                                                             <  1:277420/150‑1 (MQ=255)
                    cAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGgaga                                                                                                           <  1:2010532/150‑1 (MQ=255)
                    cAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGgaga                                                                                                           <  1:31919/150‑1 (MQ=255)
                        aaTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGt                                                                                                       <  2:1381451/150‑1 (MQ=255)
                         aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  1:1406859/150‑1 (MQ=255)
                         aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      >  2:2988379/1‑150 (MQ=255)
                         aTGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTg                                                                                                      <  2:1541268/150‑1 (MQ=255)
                             gaATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCtt                                                                                                  <  2:3215319/150‑1 (MQ=255)
                                 gATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGta                                                                                              <  1:4549921/150‑1 (MQ=255)
                                 gATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGta                                                                                              >  1:854068/1‑150 (MQ=255)
                                  aTTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGtat                                                                                             >  2:2114617/1‑150 (MQ=255)
                                    tGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTatat                                                                                           <  1:1876792/150‑1 (MQ=255)
                                    tGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTatat                                                                                           >  2:799524/1‑150 (MQ=255)
                                      gCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTg                                                                                         >  1:3299460/1‑150 (MQ=255)
                                      gCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTg                                                                                         <  2:3714509/150‑1 (MQ=255)
                                       ccTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTga                                                                                        >  2:2410635/1‑150 (MQ=255)
                                       ccTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTga                                                                                        >  1:2258552/1‑150 (MQ=255)
                                         ttcttTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaag                                                                                      <  1:4114450/150‑1 (MQ=255)
                                            tttATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAAt                                                                                   >  2:692480/1‑150 (MQ=255)
                                             ttATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATa                                                                                  <  2:4138218/150‑1 (MQ=255)
                                             ttATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATa                                                                                  >  1:4012052/1‑150 (MQ=255)
                                              tATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATAc                                                                                 >  2:2259910/1‑150 (MQ=255)
                                               aTGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACa                                                                                >  2:1661587/1‑150 (MQ=255)
                                               aTGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACa                                                                                >  2:1445766/1‑150 (MQ=255)
                                               aTGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACa                                                                                >  2:4076821/1‑150 (MQ=255)
                                                tGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACAt                                                                               <  1:2964255/150‑1 (MQ=255)
                                                 gATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATc                                                                              >  1:1149088/1‑150 (MQ=255)
                                                 gATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATc                                                                              >  1:3892208/1‑150 (MQ=255)
                                                  aTTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATcc                                                                             <  2:3147032/150‑1 (MQ=255)
                                                   ttttAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCg                                                                            <  2:176756/150‑1 (MQ=255)
                                                       aaGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTg                                                                        <  1:2852225/150‑1 (MQ=255)
                                                        aGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGt                                                                       <  1:1552891/150‑1 (MQ=255)
                                                        aGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGt                                                                       <  2:2648881/150‑1 (MQ=255)
                                                         ggATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGtt                                                                      >  1:4490844/1‑150 (MQ=255)
                                                           aTTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTtct                                                                    <  2:1550491/150‑1 (MQ=255)
                                                            ttATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTtctc                                                                   <  2:3548694/150‑1 (MQ=255)
                                                             tATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTctct                                                                  <  1:4252502/150‑1 (MQ=255)
                                                             tATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTctct                                                                  >  2:4481395/1‑150 (MQ=255)
                                                             tATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTctct                                                                  <  1:1045847/150‑1 (MQ=255)
                                                              aTGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTa                                                                 >  2:3384994/1‑150 (MQ=255)
                                                                gCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaa                                                               >  1:919457/1‑150 (MQ=255)
                                                                 cTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaaa                                                              >  2:173219/1‑150 (MQ=255)
                                                                   tCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaaaaa                                                            >  2:599788/1‑150 (MQ=255)
                                                                   tCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaaaaa                                                            <  2:4110366/150‑1 (MQ=255)
                                                                     tAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaaga                                                          <  2:3401344/150‑1 (MQ=255)
                                                                      aaaaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaagaa                                                         >  1:2572724/1‑150 (MQ=255)
                                                                        aaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAagaaga                                                       <  2:4528124/150‑1 (MQ=255)
                                                                        aaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAagaaga                                                       <  1:2039275/150‑1 (MQ=255)
                                                                        aaTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAagaaga                                                       >  1:2622524/1‑150 (MQ=255)
                                                                         aTGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAagag                                                      <  1:3415295/150‑1 (MQ=255)
                                                                          tGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGt                                                     <  2:401552/150‑1 (MQ=255)
                                                                           ggCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGtt                                                    >  1:1797175/1‑150 (MQ=255)
                                                                           ggCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGtt                                                    <  2:2169856/150‑1 (MQ=255)
                                                                             cGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGtttt                                                  <  2:4579716/150‑1 (MQ=255)
                                                                              ggAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTc                                                 >  1:592255/1‑150 (MQ=255)
                                                                                agaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCtt                                               <  1:3497289/150‑1 (MQ=255)
                                                                                agaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCtt                                               <  1:235421/150‑1 (MQ=255)
                                                                                 gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                              <  1:2783634/150‑1 (MQ=255)
                                                                                  aaCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTa                                             <  1:3708574/150‑1 (MQ=255)
                                                                                  aaCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTa                                             <  1:3821254/150‑1 (MQ=255)
                                                                                  aaCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTa                                             >  1:3854968/1‑150 (MQ=255)
                                                                                  aaCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTa                                             <  1:2372787/150‑1 (MQ=255)
                                                                                     aTGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGAt                                          >  2:563032/1‑150 (MQ=255)
                                                                                         cGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGttt                                      >  2:753330/1‑150 (MQ=255)
                                                                                          gAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTc                                     <  2:3015387/150‑1 (MQ=255)
                                                                                            ggCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAg                                   >  1:4239087/1‑150 (MQ=255)
                                                                                            ggCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTCTCAg                                   <  1:2508750/150‑1 (MQ=255)
                                                                                               ttGCTGCCTTACTTCATTACTTCAGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                >  2:2243020/1‑150 (MQ=255)
                                                                                               ttGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                <  2:1691635/150‑1 (MQ=255)
                                                                                               ttGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                <  2:4157153/150‑1 (MQ=255)
                                                                                                tgctgcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGc                               >  1:1940586/1‑150 (MQ=255)
                                                                                                  ctgcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAg                             <  2:3635169/150‑1 (MQ=255)
                                                                                                  cggcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAg                             <  2:3868483/148‑1 (MQ=255)
                                                                                                     ccTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATg                          >  2:4434260/1‑150 (MQ=255)
                                                                                                      cTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGa                         >  2:4460855/1‑150 (MQ=255)
                                                                                                      cTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGa                         >  2:3039995/1‑150 (MQ=255)
                                                                                                          cttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCtt                     <  1:3132155/150‑1 (MQ=255)
                                                                                                            tcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTa                   >  2:3055134/1‑150 (MQ=255)
                                                                                                            tcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTa                   <  1:3479882/150‑1 (MQ=255)
                                                                                                             cattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAt                  >  2:1369353/1‑150 (MQ=255)
                                                                                                                 acttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgtt              <  1:617032/150‑1 (MQ=255)
                                                                                                                  cttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCCTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgttg             <  1:2443706/150‑1 (MQ=255)
                                                                                                                   ttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAATAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTgg            >  2:140507/1‑150 (MQ=255)
                                                                                                                   ttcaGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTgg            >  1:285367/1‑150 (MQ=255)
                                                                                                                     caGCGGTGATTGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGAt          >  1:559397/1‑150 (MQ=255)
                                                                                                                          gTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAg     >  2:1121215/1‑150 (MQ=255)
                                                                                                                            gATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGAt   >  2:4260036/1‑150 (MQ=255)
                                                                                                                             aTGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGAtt  <  2:1710674/150‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                   
ATTCTACAACCAGGTCGAGTCAGAAATGAGAATGATTGGCCTTCTTTATGATTTTAAGGATTATGCTTCTAAAATGGCGGAGAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATT  >  CP009273/2753297‑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: