Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 291,453 A→T 16.2% 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,4530AT16.2% 755.1 / 124.0 276L452* (TTA→TAA) intFCP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base A (129/102);  new base T (26/19);  total (155/121)
Fisher's exact test for biased strand distribution p-value = 8.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

GTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCA  >  CP009273/291304‑291602
                                                                                                                                                     |                                                                                                                                                     
gtgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                                                                                                                                       <  1:2950962/150‑1 (MQ=255)
 tgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTaa                                                                                                                                                      <  1:3479728/150‑1 (MQ=255)
   tAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGa                                                                                                                                                    <  1:3538139/150‑1 (MQ=255)
     gTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATc                                                                                                                                                  <  1:1400145/150‑1 (MQ=255)
       tGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCa                                                                                                                                                >  2:1534417/1‑150 (MQ=255)
        gCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCAt                                                                                                                                               >  1:4042097/1‑150 (MQ=255)
        gCTGTATCATTGACTAAACTGTGTATTACGATACATAAAAATCTATCTTGATCAATAGTTGATGATATTTTTTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCAt                                                                                                                                               >  2:543437/1‑150 (MQ=255)
         cTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATg                                                                                                                                              <  2:2117545/150‑1 (MQ=255)
          tGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGt                                                                                                                                             >  2:2978443/1‑150 (MQ=255)
           gTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtc                                                                                                                                            >  2:3900342/1‑150 (MQ=255)
           gTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtc                                                                                                                                            >  2:3227026/1‑150 (MQ=255)
            tATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtct                                                                                                                                           >  1:41336/1‑150 (MQ=255)
             aTCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtctc                                                                                                                                          >  1:646433/1‑150 (MQ=255)
               cATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACATGCTTCATTTAATATCCAtttataaa                                                                                                                                        >  1:2557067/1‑143 (MQ=255)
                aTTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAAt                                                                                                                                       >  1:1493941/1‑150 (MQ=255)
                  tGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAAttt                                                                                                                                     >  1:1325281/1‑150 (MQ=255)
                     cTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                                                                                                                                  >  2:3552121/1‑150 (MQ=255)
                      tAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtctc                                                                                                                                 <  2:3871543/150‑1 (MQ=255)
                         aCTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTg                                                                                                                              <  2:2477270/150‑1 (MQ=255)
                                aTTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCAtttt                                                                                                                       >  1:3722091/1‑150 (MQ=255)
                                  tACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCg                                                                                                                     <  2:3271437/150‑1 (MQ=255)
                                   aCGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGt                                                                                                                    <  1:3714664/150‑1 (MQ=255)
                                     gATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTATCTTGCCATTTTCGTaa                                                                                                                  >  2:2361916/1‑150 (MQ=255)
                                      aTACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAAc                                                                                                                 <  2:594057/150‑1 (MQ=255)
                                          aTGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGAt                                                                                                             >  1:1208879/1‑150 (MQ=255)
                                             aaaaTCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTa                                                                                                          >  2:2637116/1‑150 (MQ=255)
                                                atctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACg                                                                                                       >  2:4260073/1‑150 (MQ=255)
                                                 tctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACgg                                                                                                      >  1:4507742/1‑150 (MQ=255)
                                                 tctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACgg                                                                                                      <  1:2369419/150‑1 (MQ=255)
                                                   tatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGgac                                                                                                    <  1:3056999/150‑1 (MQ=255)
                                                    atctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGgacg                                                                                                   <  2:1430597/150‑1 (MQ=255)
                                                    atctTGATCAATAGTTGATGATATTTTGTAATATGTCGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGgacg                                                                                                   >  2:1193258/1‑150 (MQ=255)
                                                       ttGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGacgacg                                                                                                >  1:1008492/1‑150 (MQ=255)
                                                       ttGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGacgacg                                                                                                <  1:476678/150‑1 (MQ=255)
                                                        tGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACgg                                                                                               >  2:1841439/1‑150 (MQ=255)
                                                          aTCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCg                                                                                             >  2:2081183/1‑150 (MQ=255)
                                                           tCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                                                            <  1:2541852/150‑1 (MQ=255)
                                                             aaTAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                                                          >  1:550847/1‑150 (MQ=255)
                                                                aGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                                                       >  2:2193909/1‑150 (MQ=255)
                                                                aGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                                                       >  2:364527/1‑150 (MQ=255)
                                                                  ttGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtt                                                                                     <  1:515200/150‑1 (MQ=255)
                                                                   tgatgaTATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttt                                                                                    >  1:818323/1‑150 (MQ=255)
                                                                    gatgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                                                                   >  2:2261602/1‑150 (MQ=255)
                                                                     atgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttttt                                                                                  >  1:1767866/1‑150 (MQ=255)
                                                                      tgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                                                                 <  1:1631746/150‑1 (MQ=255)
                                                                       gatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCg                                                                                <  2:1787765/150‑1 (MQ=255)
                                                                       gatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCg                                                                                >  1:3996065/1‑150 (MQ=255)
                                                                         tatTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGc                                                                              >  1:1033925/1‑150 (MQ=255)
                                                                          atTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGAAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCa                                                                             >  2:1246207/1‑150 (MQ=255)
                                                                             ttGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCTGCAaga                                                                          <  2:4160616/150‑1 (MQ=255)
                                                                             ttGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                                                          >  1:2901170/1‑150 (MQ=255)
                                                                              tGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                                                         <  1:3175820/150‑1 (MQ=255)
                                                                                tAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                                                       <  1:179421/150‑1 (MQ=255)
                                                                                tAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                                                       <  1:3953144/150‑1 (MQ=255)
                                                                                tAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                                                       >  1:1523709/1‑150 (MQ=255)
                                                                                 aaTATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                                                      >  1:1452197/1‑150 (MQ=255)
                                                                                 aaTATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                                                      >  1:3772695/1‑150 (MQ=255)
                                                                                 aaTATGTTGATTTTTAATGGAATCAACGCATATCAACGATGTTTTTTATGGTACTACCTGATTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                                                      >  1:4592056/1‑150 (MQ=255)
                                                                                  atatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                                                     <  1:2184413/150‑1 (MQ=255)
                                                                                  atatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                                                     <  2:104057/150‑1 (MQ=255)
                                                                                   tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTg                                                                    <  2:4263723/150‑1 (MQ=255)
                                                                                   tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTg                                                                    >  2:53833/1‑150 (MQ=255)
                                                                                    atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                                                                   >  2:2907936/1‑150 (MQ=255)
                                                                                        tGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGttt                                                               <  1:2743205/150‑1 (MQ=255)
                                                                                         gATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtg                                                              >  2:4522719/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           >  2:4447368/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           <  1:2561454/150‑1 (MQ=255)
                                                                                             tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                                                          >  1:3412731/1‑150 (MQ=255)
                                                                                                 aTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                                                      <  2:1460606/150‑1 (MQ=255)
                                                                                                  tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                                                     >  2:3998850/1‑150 (MQ=255)
                                                                                                   ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                                                    <  2:891711/150‑1 (MQ=255)
                                                                                                    gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt                                                   >  2:407190/1‑150 (MQ=255)
                                                                                                    gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt                                                   >  1:88631/1‑150 (MQ=255)
                                                                                                      aTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt                                                 >  2:1362129/1‑150 (MQ=255)
                                                                                                        cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa                                               >  1:4418245/1‑150 (MQ=255)
                                                                                                        cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa                                               >  1:4322620/1‑150 (MQ=255)
                                                                                                         aaCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAg                                              >  2:1075707/1‑150 (MQ=255)
                                                                                                         aaCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAg                                              <  1:2833345/150‑1 (MQ=255)
                                                                                                         aaCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAg                                              >  1:3328583/1‑150 (MQ=255)
                                                                                                           cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca                                            <  1:3779071/150‑1 (MQ=255)
                                                                                                            gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaa                                           >  2:1115845/1‑150 (MQ=255)
                                                                                                            gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaa                                           >  1:2066901/1‑150 (MQ=255)
                                                                                                                atCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAAt                                       <  2:2557074/150‑1 (MQ=255)
                                                                                                                 tCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa                                      >  1:12244/1‑150 (MQ=255)
                                                                                                                 tCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa                                      >  1:3228685/1‑150 (MQ=255)
                                                                                                                  cAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATAc                                     <  1:312780/150‑1 (MQ=255)
                                                                                                                   aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATAcc                                    >  2:3799832/1‑150 (MQ=255)
                                                                                                                   aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATAcc                                    >  2:3086038/1‑150 (MQ=255)
                                                                                                                    aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCa                                   <  1:186213/150‑1 (MQ=255)
                                                                                                                     cGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAg                                  >  2:2737750/1‑150 (MQ=255)
                                                                                                                         gttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTggg                              <  2:4592528/150‑1 (MQ=255)
                                                                                                                             tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACt                          <  1:760298/150‑1 (MQ=255)
                                                                                                                                 tGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCaa                      >  1:4303832/1‑150 (MQ=255)
                                                                                                                                     aaTACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCa                  >  2:798825/1‑150 (MQ=255)
                                                                                                                                     aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCa                  <  1:3346056/150‑1 (MQ=255)
                                                                                                                                        aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTc               <  2:2925630/150‑1 (MQ=255)
                                                                                                                                          cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCgg             >  1:1966191/1‑150 (MQ=255)
                                                                                                                                              ttCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCaa         <  2:2181720/150‑1 (MQ=255)
                                                                                                                                                    tAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTAcc   >  1:2597663/1‑150 (MQ=255)
                                                                                                                                                     aaGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCa  <  1:2676870/150‑1 (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: