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,4940GT98.5% 644.4 / 0.7 199C75* (TGC→TGAyafWCP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (96/100);  total (98/101)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 200 total aligned reads displayed.

GCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAGT  >  CP009273/259349‑259641
                                                                                                                                                 |                                                                                                                                                   
gcccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCa                                                                                                                                                   <  1:3355300/148‑1 (MQ=255)
      tGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGc                                                                                                                                             >  1:1880612/1‑148 (MQ=255)
       gAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCa                                                                                                                                            <  1:2659847/148‑1 (MQ=255)
       gAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCa                                                                                                                                            >  2:756245/1‑148 (MQ=255)
       gAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCa                                                                                                                                            >  1:505373/1‑148 (MQ=255)
          ttATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                         >  1:1128267/1‑148 (MQ=255)
           tATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggt                                                                                                                                        >  2:2306273/1‑148 (MQ=255)
            aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                       <  1:520774/148‑1 (MQ=255)
            aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                       >  2:3204796/1‑148 (MQ=255)
            aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                       >  2:3490846/1‑148 (MQ=255)
            aTTGCACGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtggtg                                                                                                                                       >  1:975965/1‑148 (MQ=255)
              tGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                     >  2:1199999/1‑148 (MQ=255)
                cAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTACCTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggggg                                                                                                                                   <  2:1370466/148‑1 (MQ=255)
                    tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                               <  1:1301465/148‑1 (MQ=255)
                    tAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGa                                                                                                                               <  2:30038/148‑1 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                              >  1:689085/1‑148 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                              <  2:896286/148‑1 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCAAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTCCAGGGTCACTCAGTGGGCCTGGTGGGGGCTGAg                                                                                                                              >  2:3229480/1‑148 (MQ=25)
                      aaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGt                                                                                                                             <  1:2729908/148‑1 (MQ=255)
                       aGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGtt                                                                                                                            <  1:1197839/148‑1 (MQ=255)
                        gTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTc                                                                                                                           >  1:1006903/1‑148 (MQ=255)
                         ttttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                          <  2:3435301/148‑1 (MQ=255)
                         ttttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                          <  1:3280777/148‑1 (MQ=255)
                         tggtCATGGTTGCTCTCGTGGAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                          <  2:2889300/145‑1 (MQ=255)
                            tCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCg                                                                                                                       >  2:3319104/1‑148 (MQ=255)
                              aTGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggt                                                                                                                     <  1:29464/148‑1 (MQ=255)
                                 tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                                  >  1:3412088/1‑148 (MQ=255)
                                 tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                                  <  2:2696948/148‑1 (MQ=255)
                                 tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                                  >  1:2778594/1‑148 (MQ=255)
                                      tctcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCa                                                                                                             <  2:2650907/148‑1 (MQ=255)
                                      tctcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCa                                                                                                             <  2:1717302/148‑1 (MQ=255)
                                       ctcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCAt                                                                                                            <  2:50174/148‑1 (MQ=255)
                                       ctcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCAt                                                                                                            >  2:2156539/1‑148 (MQ=255)
                                        tcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATc                                                                                                           <  2:3372764/148‑1 (MQ=255)
                                           tgtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTACGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGAt                                                                                                        >  2:2339482/1‑148 (MQ=255)
                                              aGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTcc                                                                                                     <  2:1175332/148‑1 (MQ=255)
                                              aGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTcc                                                                                                     >  1:3502160/1‑148 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                                >  1:3117802/1‑148 (MQ=255)
                                                       aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                            <  2:484155/148‑1 (MQ=255)
                                                        cGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagga                                                                                           >  1:3423029/1‑148 (MQ=255)
                                                        cGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagga                                                                                           >  1:371123/1‑148 (MQ=255)
                                                          cTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATa                                                                                         <  2:310107/148‑1 (MQ=255)
                                                            gAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaaa                                                                                       <  2:2524486/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                      <  2:2855136/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                      <  1:3040221/148‑1 (MQ=255)
                                                              gTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGt                                                                                     <  2:1499947/148‑1 (MQ=255)
                                                              gTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGt                                                                                     <  1:2637167/148‑1 (MQ=255)
                                                                   gTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCggg                                                                                >  2:210796/1‑148 (MQ=255)
                                                                    tAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGa                                                                               <  2:1291298/148‑1 (MQ=255)
                                                                    tAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGa                                                                               <  1:1115608/148‑1 (MQ=255)
                                                                      aaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAAt                                                                             >  2:2554699/1‑148 (MQ=255)
                                                                       aTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATg                                                                            <  1:1702411/148‑1 (MQ=255)
                                                                           cGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTc                                                                        <  2:3224484/148‑1 (MQ=255)
                                                                           cGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTc                                                                        >  1:295911/1‑148 (MQ=255)
                                                                             aTGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCAt                                                                      <  1:1134115/148‑1 (MQ=255)
                                                                                  taCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                                 <  1:770046/148‑1 (MQ=255)
                                                                                     gTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTc                                                              <  2:297123/148‑1 (MQ=255)
                                                                                       aGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCg                                                            <  2:2756548/148‑1 (MQ=255)
                                                                                       aGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCg                                                            >  2:416263/1‑148 (MQ=255)
                                                                                       aGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCg                                                            >  2:3266551/1‑148 (MQ=255)
                                                                                         ccGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGc                                                          >  1:176656/1‑148 (MQ=255)
                                                                                          cGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCa                                                         <  1:1453149/148‑1 (MQ=255)
                                                                                              aCTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCg                                                     <  1:1568314/148‑1 (MQ=255)
                                                                                               cTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGa                                                    <  2:2028444/148‑1 (MQ=255)
                                                                                                  ccAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGt                                                 <  2:708193/148‑1 (MQ=255)
                                                                                                  ccAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGt                                                 >  2:2473843/1‑148 (MQ=255)
                                                                                                   cAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTc                                                >  1:3265900/1‑148 (MQ=255)
                                                                                                     aGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCAc                                              <  1:427687/148‑1 (MQ=255)
                                                                                                     aGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCAc                                              >  1:333478/1‑148 (MQ=255)
                                                                                                     aGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCAc                                              <  1:2196069/148‑1 (MQ=255)
                                                                                                      gggTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACt                                             >  1:2514137/1‑148 (MQ=255)
                                                                                                          gtCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAAc                                         <  1:2146411/148‑1 (MQ=255)
                                                                                                           tCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACt                                        <  1:297149/148‑1 (MQ=255)
                                                                                                           tCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACt                                        <  1:297154/148‑1 (MQ=255)
                                                                                                           tCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACt                                        <  1:1332301/148‑1 (MQ=255)
                                                                                                            cGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACtt                                       >  1:2867226/1‑148 (MQ=255)
                                                                                                             ggCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTg                                      >  2:648505/1‑148 (MQ=255)
                                                                                                             ggCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTg                                      >  2:647851/1‑148 (MQ=255)
                                                                                                             ggCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTg                                      >  1:3152717/1‑148 (MQ=255)
                                                                                                              gCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGc                                     <  2:2614434/148‑1 (MQ=255)
                                                                                                                ttCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCg                                   >  1:2167492/1‑148 (MQ=255)
                                                                                                                   gCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTg                                >  1:3295344/1‑148 (MQ=255)
                                                                                                                    cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                               >  2:277764/1‑148 (MQ=255)
                                                                                                                    cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                               >  1:2849532/1‑148 (MQ=255)
                                                                                                                      ggTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATg                             >  2:960938/1‑148 (MQ=255)
                                                                                                                      ggTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATg                             <  2:1128258/148‑1 (MQ=255)
                                                                                                                        ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                           >  1:233455/1‑148 (MQ=255)
                                                                                                                        ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                           >  1:1538898/1‑148 (MQ=255)
                                                                                                                        ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                           >  2:130139/1‑148 (MQ=255)
                                                                                                                              cGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCg                     >  1:80487/1‑148 (MQ=255)
                                                                                                                                 tgtgGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTc                  <  2:3163482/148‑1 (MQ=255)
                                                                                                                                     gTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGcc              <  2:1040878/148‑1 (MQ=255)
                                                                                                                                        cAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGa           >  1:1260635/1‑148 (MQ=255)
                                                                                                                                          gggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATa         <  1:3197580/148‑1 (MQ=255)
                                                                                                                                           ggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATaa        <  2:3139465/148‑1 (MQ=255)
                                                                                                                                           ggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATaa        >  2:2153450/1‑148 (MQ=255)
                                                                                                                                           ggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATaa        >  2:1897590/1‑148 (MQ=255)
                                                                                                                                            gTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAAt       <  1:1489689/148‑1 (MQ=255)
                                                                                                                                             tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATg      <  1:98225/148‑1 (MQ=255)
                                                                                                                                              aaCTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGc     <  1:656402/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAGT  >  CP009273/259349‑259641

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: