Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 259,494 G→T 100% C75* (TGC→TGA yafW ← CP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273259,4940GT97.6% 660.0 / 6.0 208C75* (TGC→TGAyafWCP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (98/104);  total (103/105)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 208 total aligned reads displayed.

CCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGC  >  CP009273/259350‑259638
                                                                                                                                                |                                                                                                                                                
cccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAg                                                                                                                                               >  1:1901736/1‑148 (MQ=255)
     tGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGc                                                                                                                                          <  1:3736081/148‑1 (MQ=255)
         ttATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                      >  2:3010859/1‑148 (MQ=255)
           aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                    <  2:3585330/148‑1 (MQ=255)
           aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                    <  1:436799/148‑1 (MQ=255)
              gCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTgggg                                                                                                                                 <  2:2486623/148‑1 (MQ=255)
               cAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggggg                                                                                                                                <  2:847890/148‑1 (MQ=255)
               cAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggggg                                                                                                                                >  1:3535155/1‑148 (MQ=255)
                 ggTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCt                                                                                                                              >  2:1310023/1‑148 (MQ=255)
                  gTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTg                                                                                                                             <  1:536618/148‑1 (MQ=255)
                   tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCGGGTTAAATGGTTGTGGGACAGGGTAACTCCGTGGGCATGGTgggggggga                                                                                                                            >  2:2100201/1‑144 (MQ=25)
                   tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                            >  1:1099881/1‑148 (MQ=255)
                   tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                            >  2:1113755/1‑148 (MQ=255)
                   tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                            >  1:2358990/1‑148 (MQ=255)
                   tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                            >  1:1670870/1‑148 (MQ=255)
                    aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                           <  2:2635083/148‑1 (MQ=255)
                    aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAGGGTGGGGGCTGAg                                                                                                                           >  1:2351684/1‑148 (MQ=255)
                        ttttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                       >  1:1005940/1‑148 (MQ=255)
                        ttttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAACGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACACAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                       >  1:158128/1‑148 (MQ=255)
                         tttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCAc                                                                                                                      <  2:1099541/148‑1 (MQ=255)
                              tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                 <  2:1582778/148‑1 (MQ=255)
                               gTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtgg                                                                                                                >  1:2112755/1‑148 (MQ=255)
                               gTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtgg                                                                                                                >  2:140697/1‑148 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               <  2:2668429/148‑1 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               <  1:513351/148‑1 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               <  1:2597616/148‑1 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               >  1:565729/1‑148 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               >  2:3410046/1‑148 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               >  1:311224/1‑148 (MQ=255)
                                tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                               <  1:1645724/148‑1 (MQ=255)
                                       tcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATc                                                                                                        >  2:403788/1‑148 (MQ=255)
                                         gtgtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGa                                                                                                      <  1:3650590/148‑1 (MQ=255)
                                         gtgtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGa                                                                                                      <  1:1195555/148‑1 (MQ=255)
                                          tgtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGAt                                                                                                     <  2:1157963/148‑1 (MQ=255)
                                            tAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTc                                                                                                   <  2:3289896/148‑1 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 >  2:1679782/1‑148 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 >  1:1413410/1‑148 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 >  2:2235181/1‑148 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 >  2:308127/1‑148 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 >  2:561180/1‑148 (MQ=255)
                                              ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                 <  2:1351181/148‑1 (MQ=255)
                                               gTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCAt                                                                                                <  2:2964224/148‑1 (MQ=255)
                                               gTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCAt                                                                                                <  1:3532122/148‑1 (MQ=255)
                                               gTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCTATTCCAt                                                                                                >  2:1092684/1‑148 (MQ=255)
                                                 aaTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCt                                                                                              <  2:1606258/148‑1 (MQ=255)
                                                   ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                            <  2:1901751/148‑1 (MQ=255)
                                                   ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                            >  1:2208449/1‑148 (MQ=255)
                                                   ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                            >  1:2285668/1‑148 (MQ=255)
                                                   ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGATTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGATTTCACCGGTGGTAAGCATCGATTCCATCTgg                                                                                            <  2:3460341/148‑1 (MQ=255)
                                                      aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                         <  1:1490820/148‑1 (MQ=255)
                                                       cGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagga                                                                                        <  1:2781896/148‑1 (MQ=255)
                                                           gAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaaa                                                                                    <  1:2471558/148‑1 (MQ=255)
                                                              tGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTg                                                                                 >  1:2274839/1‑148 (MQ=255)
                                                              tGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTg                                                                                 <  1:2065207/148‑1 (MQ=255)
                                                              tGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTg                                                                                 <  1:2333428/148‑1 (MQ=255)
                                                               ggggTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGc                                                                                <  1:2869246/148‑1 (MQ=255)
                                                                gggTAACTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCg                                                                               <  2:2564256/148‑1 (MQ=255)
                                                                    aaaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGaa                                                                           >  2:290190/1‑148 (MQ=255)
                                                                     aaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGact                                                                          >  2:3736696/1‑146 (MQ=255)
                                                                      aTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATg                                                                         <  1:2404107/148‑1 (MQ=255)
                                                                         gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGttt                                                                      >  2:180115/1‑148 (MQ=255)
                                                                            aTGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCAt                                                                   >  1:1410657/1‑148 (MQ=255)
                                                                              gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                 <  2:2985500/148‑1 (MQ=255)
                                                                              gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                 >  1:2640682/1‑148 (MQ=255)
                                                                              gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                 >  1:2641068/1‑148 (MQ=255)
                                                                              gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAGGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                 >  2:2020198/1‑148 (MQ=255)
                                                                                 taCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                              >  2:508605/1‑148 (MQ=255)
                                                                                 taCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                              >  2:1513305/1‑148 (MQ=255)
                                                                                    gTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTc                                                           <  2:180721/148‑1 (MQ=255)
                                                                                    gTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTc                                                           <  2:2843455/148‑1 (MQ=255)
                                                                                       gCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCgg                                                        >  1:3374258/1‑148 (MQ=255)
                                                                                        ccGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGc                                                       <  1:3544404/148‑1 (MQ=255)
                                                                                         cGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCa                                                      <  1:2741798/148‑1 (MQ=255)
                                                                                             aCTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCCTGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCg                                                  >  2:749109/1‑148 (MQ=255)
                                                                                              cTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGa                                                 <  2:2837916/148‑1 (MQ=255)
                                                                                                 ccAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGt                                              <  1:3190992/148‑1 (MQ=255)
                                                                                                   aaGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCa                                            >  2:570245/1‑148 (MQ=255)
                                                                                                   aaGGGTGTCGGCTTCGAAGGTTAAATCGTTGTGGTACAGGGAAACTCAGTGGGCAGGGGGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCa                                            >  2:198126/1‑148 (MQ=255)
                                                                                                       gtgtCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGa                                        <  2:192587/148‑1 (MQ=255)
                                                                                                              cTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGcc                                 >  1:811982/1‑148 (MQ=255)
                                                                                                                  gCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTg                             <  1:1261751/148‑1 (MQ=255)
                                                                                                                   cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                            <  1:3106303/148‑1 (MQ=255)
                                                                                                                   cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                            <  1:756474/148‑1 (MQ=255)
                                                                                                                   cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                            <  1:872016/148‑1 (MQ=255)
                                                                                                                     ggTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATg                          <  1:3702217/148‑1 (MQ=255)
                                                                                                                      gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                         >  2:102950/1‑148 (MQ=255)
                                                                                                                      gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                         >  1:2554490/1‑148 (MQ=255)
                                                                                                                      gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                         <  1:180115/148‑1 (MQ=255)
                                                                                                                      gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                         >  2:2447928/1‑148 (MQ=255)
                                                                                                                       ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                        <  1:290190/148‑1 (MQ=255)
                                                                                                                       ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                        <  2:792091/148‑1 (MQ=255)
                                                                                                                        tAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTg                       >  1:1512120/1‑148 (MQ=255)
                                                                                                                        tAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTg                       >  1:3420855/1‑148 (MQ=255)
                                                                                                                        tAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGAGTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTg                       <  2:3675446/148‑1 (MQ=255)
                                                                                                                            tCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGccc                   >  1:2261084/1‑148 (MQ=255)
                                                                                                                             cGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCg                  >  2:327929/1‑148 (MQ=255)
                                                                                                                                     tACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCa          >  2:3593574/1‑148 (MQ=255)
                                                                                                                                       cAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGa        <  1:3410046/148‑1 (MQ=255)
                                                                                                                                         gggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATa      <  2:1876626/148‑1 (MQ=255)
                                                                                                                                             aaCTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGc  <  2:2275114/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGC  >  CP009273/259350‑259638

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: