Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 291,453 A→T 24.6% 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,4530AT24.6% 561.6 / 216.2 272L452* (TTA→TAA) intFCP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base A (134/71);  new base T (46/21);  total (180/92)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 279 total aligned reads displayed.

GTGTAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCA  >  CP009273/291304‑291602
                                                                                                                                                     |                                                                                                                                                     
gtgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                                                                                                                                       <  1:1745860/150‑1 (MQ=255)
gtgtAGTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTa                                                                                                                                                       >  1:4549118/1‑150 (MQ=255)
     gTTGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATc                                                                                                                                                  <  1:4286743/150‑1 (MQ=255)
      ttGCTGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATcc                                                                                                                                                 >  2:4202264/1‑150 (MQ=255)
          tGTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCAAGt                                                                                                                                             >  1:4359328/1‑150 (MQ=255)
           gTATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtc                                                                                                                                            <  1:3186400/150‑1 (MQ=255)
            tATCATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGtct                                                                                                                                           >  1:3487252/1‑150 (MQ=255)
               cATTGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCaa                                                                                                                                        >  2:3673426/1‑150 (MQ=255)
                 ttGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAAtt                                                                                                                                      >  2:818442/1‑150 (MQ=255)
                 ttGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAAtt                                                                                                                                      <  1:145825/150‑1 (MQ=255)
                  tGACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAAttt                                                                                                                                     >  1:1400335/1‑150 (MQ=255)
                   gACTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAAtttt                                                                                                                                    <  1:1799113/150‑1 (MQ=255)
                    aCTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAATATCCATGTCTCAATTTtc                                                                                                                                   >  2:4161881/1‑150 (MQ=255)
                     cTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                                                                                                                                  >  1:1450250/1‑150 (MQ=255)
                     cTAAACTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTtct                                                                                                                                  >  2:1511055/1‑150 (MQ=255)
                       aaaCTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCg                                                                                                                                <  2:1616658/150‑1 (MQ=255)
                          cTGTGTATTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGc                                                                                                                             <  2:2063516/150‑1 (MQ=255)
                                aTTACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCAtttt                                                                                                                       >  1:153631/1‑150 (MQ=255)
                                 ttACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAAAGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTc                                                                                                                      <  2:3503960/150‑1 (MQ=255)
                                  tACGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCg                                                                                                                     >  2:4025978/1‑150 (MQ=255)
                                   aCGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGt                                                                                                                    >  1:3083715/1‑150 (MQ=255)
                                    cGATACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTa                                                                                                                   >  1:4574451/1‑150 (MQ=255)
                                      aTACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAAc                                                                                                                 >  2:3787589/1‑150 (MQ=255)
                                       tACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACa                                                                                                                >  2:3763935/1‑150 (MQ=255)
                                       tACATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACa                                                                                                                >  1:667917/1‑150 (MQ=255)
                                         cATGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGa                                                                                                              >  2:2615070/1‑150 (MQ=255)
                                          aTGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGAt                                                                                                             >  2:2338194/1‑150 (MQ=255)
                                          aTGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGAt                                                                                                             >  2:1953196/1‑150 (MQ=255)
                                          aTGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGAt                                                                                                             >  1:4429694/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            <  1:1358305/150‑1 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            >  1:4141630/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            >  2:4603315/1‑150 (MQ=255)
                                           tGAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                                                            <  1:4793212/150‑1 (MQ=255)
                                            gAAAATCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCt                                                                                                           <  2:3338518/150‑1 (MQ=255)
                                              aaaTCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTaa                                                                                                         >  2:4049655/1‑150 (MQ=255)
                                              aaaTCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTaa                                                                                                         <  1:4901924/150‑1 (MQ=255)
                                              aaaTCTATCTTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTaa                                                                                                         <  2:532739/150‑1 (MQ=255)
                                                 tctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACgg                                                                                                      <  2:4708482/150‑1 (MQ=255)
                                                 tctatctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACgg                                                                                                      >  1:2820083/1‑150 (MQ=255)
                                                     tctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGgacga                                                                                                  <  2:2172931/150‑1 (MQ=255)
                                                      ctTGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGgacgac                                                                                                 >  1:2516688/1‑150 (MQ=255)
                                                       ttGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGacgacg                                                                                                >  1:2213776/1‑150 (MQ=255)
                                                        tGATCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACgg                                                                                               <  1:77145/150‑1 (MQ=255)
                                                           tCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                                                            <  2:326120/150‑1 (MQ=255)
                                                           tCAATAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                                                            <  1:2610731/150‑1 (MQ=255)
                                                              aTAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAg                                                                                         >  1:362596/1‑150 (MQ=255)
                                                               tAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCATCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgt                                                                                        <  1:1212354/150‑1 (MQ=255)
                                                               tAGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgt                                                                                        >  2:1075459/1‑150 (MQ=255)
                                                                aGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                                                       <  1:5099999/150‑1 (MQ=255)
                                                                aGTTGATGATATTTTGTAATATGTTGATTTTTAATGGAATCAAAGCATATCAAAGTTTTTTTTTATGGTAATACCTTCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATATTCGTAACATATCTAACGGACGACGGCGATAgag                                                                                       >  1:892253/1‑148 (MQ=255)
                                                                     atgatATTTTTTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttttt                                                                                  <  2:3923320/150‑1 (MQ=255)
                                                                      tgatATTTTGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGAGAGTGTTTTTc                                                                                 >  1:1033180/1‑150 (MQ=255)
                                                                            tttGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAag                                                                           >  2:2253503/1‑150 (MQ=255)
                                                                             ttGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                                                          <  2:1207907/150‑1 (MQ=255)
                                                                             ttGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                                                          >  1:3594748/1‑150 (MQ=255)
                                                                              tGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                                                         >  1:466612/1‑150 (MQ=255)
                                                                              tGTAATATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                                                         <  1:1516792/150‑1 (MQ=255)
                                                                                 aaTATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                                                      >  2:2340913/1‑150 (MQ=255)
                                                                                 aaTATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                                                      >  2:4884338/1‑150 (MQ=255)
                                                                                      gTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGt                                                                 <  1:3439999/150‑1 (MQ=255)
                                                                                        tGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGttt                                                               <  2:182045/150‑1 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           >  1:4134336/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           >  1:1606006/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           >  1:1652702/1‑150 (MQ=255)
                                                                                            ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                                                           <  1:1898389/150‑1 (MQ=255)
                                                                                             tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                                                          >  1:4328674/1‑150 (MQ=255)
                                                                                                aaTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                                                       >  2:2707487/1‑150 (MQ=255)
                                                                                                  tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                                                     >  1:1604327/1‑150 (MQ=255)
                                                                                                    gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt                                                   <  2:3229743/150‑1 (MQ=255)
                                                                                                     aaTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAtt                                                  <  1:2625497/150‑1 (MQ=255)
                                                                                                        cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa                                               >  1:1966733/1‑150 (MQ=255)
                                                                                                          aCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGc                                             <  1:1913833/150‑1 (MQ=255)
                                                                                                           cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca                                            >  1:2958926/1‑150 (MQ=255)
                                                                                                            gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaa                                           >  1:3129481/1‑150 (MQ=255)
                                                                                                            gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaa                                           >  2:1625350/1‑150 (MQ=255)
                                                                                                             cATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaac                                          >  2:312295/1‑150 (MQ=255)
                                                                                                             cATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaac                                          >  2:4889671/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  2:2521535/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  1:2029951/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  1:4670580/1‑150 (MQ=255)
                                                                                                              atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca                                         >  1:61594/1‑150 (MQ=255)
                                                                                                                atCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAAt                                       >  1:403307/1‑150 (MQ=255)
                                                                                                                atCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAAt                                       >  2:299542/1‑150 (MQ=255)
                                                                                                                 tCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa                                      >  1:5048202/1‑150 (MQ=255)
                                                                                                                     cGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATCATTTGAGCAACAATACCAg                                  >  2:2139424/1‑150 (MQ=255)
                                                                                                                            ttttATGGTAATACCTGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAAc                           >  1:3805769/1‑150 (MQ=255)
                                                                                                                            ttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAAc                           <  1:814649/150‑1 (MQ=255)
                                                                                                                            ttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAAc                           >  2:340707/1‑150 (MQ=255)
                                                                                                                               tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTc                        >  1:377680/1‑150 (MQ=255)
                                                                                                                                   gTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAAcc                    >  2:2427604/1‑150 (MQ=255)
                                                                                                                                   gTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAAcc                    >  1:2128844/1‑150 (MQ=255)
                                                                                                                                       tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACt                >  1:2416524/1‑150 (MQ=255)
                                                                                                                                        aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTc               >  2:4348640/1‑150 (MQ=255)
                                                                                                                                        aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTc               <  1:2253503/150‑1 (MQ=255)
                                                                                                                                           tGCTTCATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGc            >  1:3680002/1‑150 (MQ=255)
                                                                                                                                             cTTAATTTAAGATCCATGTCTCAATATTCTCGTGCCATTATCGTAACAGATCTAACGGACGACGGCGATAGAGTTTTTCGGAAAGAGCGCTGGGTTAGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCa          >  2:3459213/1‑150 (MQ=255)
                                                                                                                                                cATTTTAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAg       <  1:5100667/150‑1 (MQ=255)
                                                                                                                                                cATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAg       >  1:1299455/1‑150 (MQ=255)
                                                                                                                                                  tttAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTa     <  2:667917/150‑1 (MQ=255)
                                                                                                                                                     aaGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATACCAGTGGGAACTTCAACCCACTCGGCCAAAGTACCa  >  1:3344739/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: