Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 263,112 G→T 7.8% S202* (TCG→TAG)  yafZ ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273263,1120GT 7.8% 600.5 / 17.3 173S202* (TCG→TAG) yafZCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (90/68);  new base T (8/7);  total (98/75)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

ATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATA  >  CP009273/262965‑263257
                                                                                                                                                   |                                                                                                                                                 
     gcgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCa                                                                                                                                              <  1:3308211/148‑1 (MQ=255)
      cgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCTGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAg                                                                                                                                             <  1:886796/148‑1 (MQ=255)
      cgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGCCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAg                                                                                                                                             <  1:1423229/148‑1 (MQ=255)
       gAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGt                                                                                                                                            <  1:2552973/148‑1 (MQ=255)
       gAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGt                                                                                                                                            <  1:1861279/148‑1 (MQ=255)
       gAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGt                                                                                                                                            <  2:3311889/148‑1 (MQ=255)
            gCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCg                                                                                                                                       >  2:452745/1‑148 (MQ=255)
             cACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                                                                                                                      <  1:1087542/148‑1 (MQ=255)
             cACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                                                                                                                      <  1:2148944/148‑1 (MQ=255)
             cACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                                                                                                                      <  1:2147674/148‑1 (MQ=255)
              aCGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGc                                                                                                                                     <  2:2281417/148‑1 (MQ=255)
                 ggTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtgg                                                                                                                                  <  2:639551/148‑1 (MQ=255)
                   tATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtggtg                                                                                                                                <  1:3182122/148‑1 (MQ=255)
                     tgtgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTggtggt                                                                                                                              >  2:2078681/1‑148 (MQ=255)
                       tgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCt                                                                                                                            >  1:3180379/1‑148 (MQ=255)
                         aCCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTc                                                                                                                          >  1:1939032/1‑148 (MQ=255)
                                tcctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAg                                                                                                                   <  2:738754/148‑1 (MQ=255)
                                tcctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAg                                                                                                                   <  1:3054179/148‑1 (MQ=255)
                                   tttGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCgg                                                                                                                <  1:1746046/148‑1 (MQ=255)
                                   tttGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCgg                                                                                                                >  1:1630142/1‑148 (MQ=255)
                                      ggCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAt                                                                                                             >  1:2438865/1‑148 (MQ=255)
                                       gCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtg                                                                                                            <  2:2081962/148‑1 (MQ=255)
                                          ttACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGa                                                                                                         >  1:1687311/1‑148 (MQ=255)
                                          ttACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCACTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGa                                                                                                         <  1:1683602/148‑1 (MQ=255)
                                           tACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAg                                                                                                        >  2:1594158/1‑148 (MQ=255)
                                           tACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAg                                                                                                        >  1:1085496/1‑148 (MQ=255)
                                           tACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAg                                                                                                        <  2:2833975/148‑1 (MQ=255)
                                               gCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCGTCACCAAAGCGGTATGTGAGggcg                                                                                                    >  1:103315/1‑148 (MQ=255)
                                                ccACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg                                                                                                   <  2:2686413/148‑1 (MQ=255)
                                                ccACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg                                                                                                   <  2:3002509/148‑1 (MQ=255)
                                                   cTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCtt                                                                                                <  2:2386375/148‑1 (MQ=255)
                                                    tGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCttt                                                                                               <  2:2039722/148‑1 (MQ=255)
                                                       gcccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGc                                                                                            >  2:25369/1‑148 (MQ=255)
                                                        cccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGcc                                                                                           >  2:2586982/1‑148 (MQ=255)
                                                        cccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGcc                                                                                           >  2:873775/1‑148 (MQ=255)
                                                        cccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGAGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGcc                                                                                           >  2:1848883/1‑148 (MQ=255)
                                                           gcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGt                                                                                        >  1:1366963/1‑148 (MQ=255)
                                                            cccTTAATCAGGTTCTCCTGAATACTCTGGTACGTGGTCCACAGGTCATTGCTCTCCTCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTTCCAGTg                                                                                       <  2:1399906/148‑1 (MQ=255)
                                                            cccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTg                                                                                       <  2:1399909/148‑1 (MQ=255)
                                                              cTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGcc                                                                                     >  1:2799537/1‑148 (MQ=255)
                                                                   tCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCctgctg                                                                                >  1:1876119/1‑148 (MQ=255)
                                                                    cAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCtgt                                                                               >  2:1713148/1‑148 (MQ=255)
                                                                        ttCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCg                                                                           >  1:855948/1‑148 (MQ=255)
                                                                           tcCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCgggg                                                                        >  2:1042683/1‑148 (MQ=255)
                                                                            ccTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCggggg                                                                       >  1:2677560/1‑148 (MQ=255)
                                                                            ccTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCggggg                                                                       >  1:1975156/1‑148 (MQ=255)
                                                                              tGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTg                                                                     >  2:2559886/1‑148 (MQ=255)
                                                                               gAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTgg                                                                    >  2:2576817/1‑148 (MQ=255)
                                                                                aaTACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGc                                                                   <  1:629700/148‑1 (MQ=255)
                                                                                  tACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCaa                                                                 >  1:556661/1‑148 (MQ=255)
                                                                                   aCGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAAc                                                                <  1:1307860/148‑1 (MQ=255)
                                                                                    cGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACa                                                               >  2:377842/1‑148 (MQ=255)
                                                                                    cGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACa                                                               >  2:377894/1‑148 (MQ=255)
                                                                                         gTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAc                                                          >  1:2083788/1‑148 (MQ=255)
                                                                                         gTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAAc                                                          >  2:1836894/1‑148 (MQ=255)
                                                                                           aCGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGa                                                        >  1:2563990/1‑148 (MQ=255)
                                                                                            cGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGAc                                                       >  1:865971/1‑148 (MQ=255)
                                                                                                  ccACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATg                                                 >  1:1988544/1‑148 (MQ=255)
                                                                                                   cacaGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATgg                                                >  1:3102260/1‑148 (MQ=255)
                                                                                                    acaGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGc                                               >  1:2966367/1‑148 (MQ=255)
                                                                                                     caGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCa                                              >  1:482076/1‑148 (MQ=255)
                                                                                                     caGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCa                                              >  1:2463913/1‑148 (MQ=255)
                                                                                                      aGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCAt                                             >  2:1684614/1‑148 (MQ=255)
                                                                                                           aTTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGt                                        >  1:3074996/1‑148 (MQ=255)
                                                                                                                tctcATCGTGCCAGCGGCGTGGGTAGAGGATCTGCGATTCAGTCACGGGGTGGTGGTCTTAACCAAAGGGGTTTGTGTGGGCGGCTTTTGCCTGTGCCTGCTGTGCCGGGGGTGCCAACAGCAACGACTGCATGGCATCCCGTTtctc                                   <  1:2675073/148‑1 (MQ=21)
                                                                                                                tctcATCCTGCCAGCGGCGTGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTtctc                                   <  2:2137974/148‑1 (MQ=255)
                                                                                                                 ctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcc                                  >  2:3304951/1‑148 (MQ=255)
                                                                                                                 ctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcc                                  >  1:3008143/1‑148 (MQ=255)
                                                                                                                 ctcATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcc                                  >  1:3008140/1‑148 (MQ=255)
                                                                                                                   caaccTGCCCGCGGCGTGGGGAGAGTATCTGCTATTCCGCCCCCGGGTGGGGGTCTACACCAACGGGGTTTGTGAGGGGGGCTTTTGCCAGGGCGTGCTGGGCGGGGGGTGCAAACAGCAACGACTGCATGGCATCCCGTTTctcctc                                <  1:452745/145‑1 (MQ=16)
                                                                                                                       cTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAACGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCGGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCAcc                            <  2:1720615/148‑1 (MQ=255)
                                                                                                                                gcgTGGGGCGAGTATGTGCGATTCAGACCCGGGGTGGTGTTTTTCACCAACGCGGTTTGTGTGGGCGGCTTTTGCCTGTGCCTGCTGGGCGGGGGGTGGAAACAGCAACGACTGCATAGCATCCCGTTTCTCCTCCACCCGGTCaaaa                   <  2:285771/148‑1 (MQ=255)
                                                                                                                                gcgAGGGGAGAGTATCTGCGATTAATTCACGGGCTGGTGGTCTTCACCAAAGCTTTTTGTGTGGGCGGCTTTTGCCTGTGCGTGCTGTGCGGGGGGTGGAAACAGCAACGACTGCATGGCAACCCGTTTCTCCTCCACCCGGTCaaaa                   <  2:566338/148‑1 (MQ=255)
                                                                                                                                 cgAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCaaaaa                  >  1:2983673/1‑148 (MQ=255)
                                                                                                                                  gAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAAt                 <  2:2130514/148‑1 (MQ=255)
                                                                                                                                  gAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAAt                 <  2:2438865/148‑1 (MQ=255)
                                                                                                                                   aGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATc                >  2:1223029/1‑148 (MQ=255)
                                                                                                                                    ggggAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAACGCGGTATGTGCGGGCGGCTTTTGCCAGGGCGTGCGGTGCCGGGGGTGGAAACCGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATcc               <  1:3242378/148‑1 (MQ=25)
                                                                                                                                     gggTGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACAAAAGGGGTTTGTGAGGGCGGCTTTTGCCTGTGCCTGCTGTGCGGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATccc              <  1:1858603/148‑1 (MQ=25)
                                                                                                                                     gggAGAGGATCTGCGATTCAGTCACGGGGTGGTGGTCTTCACCAACGCGGTTTGTGTGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGAAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATccc              <  1:2696940/148‑1 (MQ=25)
                                                                                                                                      ggAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATcccc             >  1:3223887/1‑148 (MQ=255)
                                                                                                                                      ggAGAGGATCTGCTATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATcccc             >  1:3226282/1‑148 (MQ=255)
                                                                                                                                        agagGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAg           >  1:2741769/1‑148 (MQ=255)
                                                                                                                                        agagGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAg           >  1:2741756/1‑148 (MQ=255)
                                                                                                                                        agagGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAg           >  1:2782427/1‑148 (MQ=255)
                                                                                                                                        agagGATCTGCGATTCAGTCACCGGCTGGCGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAg           >  1:192119/1‑148 (MQ=255)
                                                                                                                                          agGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCa         <  2:1908313/148‑1 (MQ=255)
                                                                                                                                           ggATCTGCGATTCAGTCCCCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCGGGGGGTGCCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCAc        <  1:186768/148‑1 (MQ=255)
                                                                                                                                           ggATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCAc        >  1:3235157/1‑148 (MQ=255)
                                                                                                                                            gATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCAcc       <  2:1801614/148‑1 (MQ=255)
                                                                                                                                            gATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCAcc       >  2:187884/1‑148 (MQ=255)
                                                                                                                                              tCTGCGATTCCGCCCCGGGGTGGTGTTCTACACCAACGCGTTTTGTGCGGGGGGCTTTTGCCAGTGCCGGGTGTGCGGGGGGGGGCAACACCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTc     <  1:1504984/148‑1 (MQ=255)
                                                                                                                                              tCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTc     >  2:3240955/1‑148 (MQ=255)
                                                                                                                                              tCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTc     <  2:3189465/148‑1 (MQ=255)
                                                                                                                                               gtGCGATTCAGTCACGGGCTGGTGGTCTTCACCAAAGCGGTTTGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa    <  2:727762/147‑1 (MQ=255)
                                                                                                                                               cTGCTATTCAGTCACGGGGTGGTGTTCTTCACCAACGGGTTATGTGTGGGCGGTTTTTGCCTGTGCCTGGTGGGCGGGGGGTGGAAACACCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa    <  1:3088142/148‑1 (MQ=255)
                                                                                                                                               cTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa    >  2:3260461/1‑148 (MQ=255)
                                                                                                                                               cTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa    >  2:3085732/1‑148 (MQ=255)
                                                                                                                                               cTGCGATTCAGCCCCGGGGTGGTGTTCTTCACAAAAGGGTTTTGTGTGGGGGGTTTTTGCCTGTGCCTGCTGGGCGGGGGGTGGAAACACCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCa    <  1:373623/148‑1 (MQ=255)
                                                                                                                                                 gCGATTCAGTCACGGGCTGGTGGTCTTCACAAAAGGGGTATGTGTGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGAAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATa  <  2:784911/148‑1 (MQ=25)
                                                                                                                                                   |                                                                                                                                                 
ATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGGTCAAAAATCCCCAGCACCTCATA  >  CP009273/262965‑263257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: