Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 291,453 A→T 29.1% L452* (TTA→TAA)  intF ← CP4‑6 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273291,4530AT29.1% 271.3 / 155.0 165L452* (TTA→TAA) intFCP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base A (74/43);  new base T (28/20);  total (102/63)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 171 total aligned reads displayed.

GTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCA  >  CP009273/291304‑291602
                                                                                                                                                     |                                                                                                                                                     
gtgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                                                                                                                                       <  1:1010523/150‑1 (MQ=255)
gtgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                                                                                                                                       <  1:2692532/150‑1 (MQ=255)
     gTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATc                                                                                                                                                  <  2:2479421/150‑1 (MQ=255)
             aTCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGtctc                                                                                                                                          <  1:2657417/150‑1 (MQ=255)
             aTCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGtctc                                                                                                                                          <  1:2657453/150‑1 (MQ=255)
                   gACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAAtttt                                                                                                                                    <  2:1218775/150‑1 (MQ=255)
                     cTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                                                                                                                                  >  1:3103382/1‑150 (MQ=255)
                      tAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTtctc                                                                                                                                 <  2:2992841/150‑1 (MQ=255)
                      tAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtctc                                                                                                                                 >  2:4544494/1‑150 (MQ=255)
                        aaCTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGt                                                                                                                               >  1:1591476/1‑150 (MQ=255)
                        aaCTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGt                                                                                                                               >  2:2638323/1‑150 (MQ=255)
                          cTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGc                                                                                                                             <  2:1726205/150‑1 (MQ=255)
                           tgtgTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGACTCAATTTTCTCGTGcc                                                                                                                            >  2:560027/1‑150 (MQ=255)
                              gtATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCAtt                                                                                                                         <  1:3518511/150‑1 (MQ=255)
                                aTTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCAtttt                                                                                                                       <  1:1275247/150‑1 (MQ=255)
                                       tACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACa                                                                                                                <  1:2547503/150‑1 (MQ=255)
                                       tACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACa                                                                                                                <  1:2547627/150‑1 (MQ=255)
                                        aCATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAg                                                                                                               >  2:4626341/1‑150 (MQ=255)
                                        aCATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAg                                                                                                               <  1:4587312/150‑1 (MQ=255)
                                         cATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACTCATATCAACGTTGTTTTTTATTGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTTCCATTTTAGTAACAGa                                                                                                              >  1:1330437/1‑150 (MQ=255)
                                         cATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTATCGTGCCATTTTCGTAACAGa                                                                                                              >  1:3029444/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            >  2:3374454/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            >  2:3374456/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATAACTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGAGCCATTTTCGTAACATATc                                                                                                            >  2:6538/1‑150 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           >  2:3497328/1‑150 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           >  1:1891358/1‑150 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           >  1:3120360/1‑150 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           >  1:730798/1‑150 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           >  1:1891016/1‑150 (MQ=255)
                                               aaTCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAAc                                                                                                        >  2:2582374/1‑150 (MQ=255)
                                                 tctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACgg                                                                                                      >  1:4229609/1‑150 (MQ=255)
                                                       ttGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGacgacg                                                                                                >  2:3016692/1‑150 (MQ=255)
                                                         gATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGc                                                                                              >  2:2686999/1‑150 (MQ=255)
                                                         gATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGc                                                                                              <  1:3304878/150‑1 (MQ=255)
                                                            cAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGAt                                                                                           >  2:4347177/1‑150 (MQ=255)
                                                            cAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGAt                                                                                           >  2:4345692/1‑150 (MQ=255)
                                                            cAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGAt                                                                                           <  1:2585061/150‑1 (MQ=255)
                                                             taTAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTTATACCTTCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                                                          <  1:4491265/149‑1 (MQ=255)
                                                             aaTAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                                                          <  1:4275769/150‑1 (MQ=255)
                                                                  ttGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtt                                                                                     >  1:3086813/1‑150 (MQ=255)
                                                                    gatgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                                                                   <  1:3632093/150‑1 (MQ=255)
                                                                    gatgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                                                                   <  1:1218049/150‑1 (MQ=255)
                                                                      tgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                                                                 >  2:3222630/1‑150 (MQ=255)
                                                                      tgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                                                                 <  1:4175707/150‑1 (MQ=255)
                                                                      tgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                                                                 <  1:4174591/150‑1 (MQ=255)
                                                                         tatTTTGTAATATATTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGc                                                                              >  2:3183989/1‑150 (MQ=255)
                                                                                  atatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                                                     >  2:1281843/1‑150 (MQ=255)
                                                                                    atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                                                                   >  1:1635338/1‑150 (MQ=255)
                                                                                    atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                                                                   >  2:2138975/1‑150 (MQ=255)
                                                                                      gTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGt                                                                 >  2:2411198/1‑150 (MQ=255)
                                                                                          aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                                                             <  1:1003502/150‑1 (MQ=255)
                                                                                          aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                                                             >  2:1179258/1‑150 (MQ=255)
                                                                                           tttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgtg                                                            >  1:3675991/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           >  1:3558728/1‑150 (MQ=255)
                                                                                             tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                                                          >  1:2480916/1‑150 (MQ=255)
                                                                                             tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                                                          >  1:2833448/1‑150 (MQ=255)
                                                                                              ttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTcc                                                         <  1:3319266/150‑1 (MQ=255)
                                                                                                aaTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                                                       <  1:534616/150‑1 (MQ=255)
                                                                                                 aTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                                                      <  1:635680/150‑1 (MQ=255)
                                                                                                  tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                                                     >  2:1722633/1‑150 (MQ=255)
                                                                                                   ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                                                    >  1:635177/1‑150 (MQ=255)
                                                                                                   ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                                                    >  1:2415283/1‑150 (MQ=255)
                                                                                                   ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                                                    >  2:1433060/1‑150 (MQ=255)
                                                                                                      aTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTTCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt                                                 <  1:270082/150‑1 (MQ=255)
                                                                                                       tCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTg                                                <  1:2802264/150‑1 (MQ=255)
                                                                                                        cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa                                               <  1:1678021/150‑1 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCCTTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTTAGcaaca                                         >  2:4093986/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  2:1007162/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  2:1005912/1‑150 (MQ=255)
                                                                                                                atCAACGTTGTTTTTTATGGTAATAACTGCTTCATTTTAGATCCATGTATCAATTTTATAGTGCCAATTTCGTAACAGATATAACCGAAGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCACATAATTTGAGcaacaaa                                       >  1:3600547/1‑149 (MQ=255)
                                                                                                                   aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCTGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAAAAATAcc                                    >  1:4440307/1‑150 (MQ=255)
                                                                                                                    aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCa                                   <  2:4635166/150‑1 (MQ=255)
                                                                                                                    aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCa                                   >  1:4541832/1‑150 (MQ=255)
                                                                                                                    aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCa                                   >  1:4538103/1‑150 (MQ=255)
                                                                                                                     cGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATCATTTGAGCAACAATACCAg                                  >  2:1877445/1‑150 (MQ=255)
                                                                                                                       ttgttTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTg                                >  1:4686634/1‑150 (MQ=255)
                                                                                                                        tgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTgg                               >  2:2813130/1‑150 (MQ=255)
                                                                                                                         gttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTggg                              <  1:845048/150‑1 (MQ=255)
                                                                                                                           tttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAACTTGAGCAACAATACCAGTGGGaa                            <  1:1113761/150‑1 (MQ=255)
                                                                                                                            ttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAAc                           >  2:3098866/1‑150 (MQ=255)
                                                                                                                              ttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACtt                         <  2:552896/150‑1 (MQ=255)
                                                                                                                               tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTc                        <  1:1198851/150‑1 (MQ=255)
                                                                                                                                aTGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCa                       >  2:3678673/1‑150 (MQ=255)
                                                                                                                                aTGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCa                       >  2:4364601/1‑150 (MQ=255)
                                                                                                                                 tGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCaa                      >  1:3667491/1‑150 (MQ=255)
                                                                                                                                    tAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAAccc                   >  2:838403/1‑150 (MQ=255)
                                                                                                                                     aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCa                  <  2:2868144/150‑1 (MQ=255)
                                                                                                                                     aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCa                  >  1:1964415/1‑150 (MQ=255)
                                                                                                                                     aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCa                  <  2:902977/150‑1 (MQ=255)
                                                                                                                                      aTACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCAc                 >  2:2873029/1‑150 (MQ=255)
                                                                                                                                         ccTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCg              >  1:383602/1‑150 (MQ=255)
                                                                                                                                             cTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCa          <  1:1544968/150‑1 (MQ=255)
                                                                                                                                               tCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCaaa        >  1:4296145/1‑150 (MQ=255)
                                                                                                                                                cATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAg       >  2:3435529/1‑150 (MQ=255)
                                                                                                                                                cATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGACCTTCAACCCACTCGGCCAAAg       >  1:413776/1‑150 (MQ=255)
                                                                                                                                                cATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAg       >  1:3201461/1‑150 (MQ=255)
                                                                                                                                                cATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAg       >  1:1549832/1‑150 (MQ=255)
                                                                                                                                                    tAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTAc  >  2:1414913/1‑149 (MQ=255)
                                                                                                                                                     aaGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCa  >  2:3787906/1‑150 (MQ=255)
                                                                                                                                                     aaGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCa  >  1:950473/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCA  >  CP009273/291304‑291602

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: