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.0% 970.8 / 5.4 300C75* (TGC→TGAyafWCP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base T (156/138);  total (161/139)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 300 total aligned reads displayed.

CCGCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAG  >  CP009273/259347‑259640
                                                                                                                                                   |                                                                                                                                                  
ccgcccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCa                                                                                                                                                  >  1:1522428/1‑150 (MQ=255)
  gcccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGt                                                                                                                                                >  1:4999373/1‑150 (MQ=255)
     cgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAACTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTggg                                                                                                                                             <  2:2238305/150‑1 (MQ=255)
     cgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTggg                                                                                                                                             <  1:4414378/150‑1 (MQ=255)
        tGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAt                                                                                                                                          >  1:3130686/1‑150 (MQ=255)
        tGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAt                                                                                                                                          <  1:4617290/150‑1 (MQ=255)
         gAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAta                                                                                                                                         <  2:4684876/150‑2 (MQ=255)
         gAGTTATAGCAGGTAAAGTTTTCATGTGTGCTCTCGTGGAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtg                                                                                                                                         <  2:3894068/150‑1 (MQ=255)
          tgTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                        >  1:2779125/2‑150 (MQ=255)
          tgTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGAAGCCGCAACTGCCAAGGGTGTCTGCTACGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                        >  1:1361152/2‑150 (MQ=255)
          aGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                        >  2:3339302/1‑150 (MQ=255)
          aGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                        >  2:3391750/1‑150 (MQ=255)
              atggCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                    <  2:1712821/147‑1 (MQ=255)
              aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGTCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                    >  2:258877/1‑150 (MQ=255)
              aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGTGCATGGTggg                                                                                                                                    >  1:3027128/1‑150 (MQ=255)
              aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                    <  1:3530556/150‑1 (MQ=255)
              aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                    <  1:2217905/150‑1 (MQ=255)
              aTTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                    >  2:2172956/1‑150 (MQ=255)
               ttGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTgggg                                                                                                                                   <  2:4022225/150‑1 (MQ=255)
                  cAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCt                                                                                                                                <  1:2962478/150‑1 (MQ=255)
                   aGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTg                                                                                                                               >  2:4225795/1‑150 (MQ=255)
                     gTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                             >  2:1848732/1‑150 (MQ=255)
                       aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGtt                                                                                                                           >  2:3310037/1‑150 (MQ=255)
                        aaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTc                                                                                                                          <  2:1818765/150‑1 (MQ=255)
                        aaGTTTTCATGGTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTc                                                                                                                          <  2:2629558/150‑1 (MQ=255)
                         aGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCa                                                                                                                         <  1:789092/150‑1 (MQ=255)
                           ttttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCAcc                                                                                                                       <  1:4054085/150‑1 (MQ=255)
                             ttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAATTCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCgg                                                                                                                     <  1:4302989/150‑1 (MQ=255)
                             ttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCgg                                                                                                                     >  1:5078882/1‑150 (MQ=255)
                               cATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                   >  1:1797510/1‑150 (MQ=255)
                               cATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                   <  2:4577313/150‑1 (MQ=255)
                                 tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                                 >  1:1441055/1‑150 (MQ=255)
                                 tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                                 >  1:2488714/1‑150 (MQ=255)
                                   tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCTCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTTGGCATGGTGGGGTCTGAGATCAGCGGTGGTCa                                                                                                               >  2:4398286/1‑150 (MQ=25)
                                   tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCa                                                                                                               <  1:1325982/150‑1 (MQ=255)
                                    ttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAg                                                                                                              >  2:3688623/1‑150 (MQ=255)
                                     tGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGc                                                                                                             >  1:3450558/1‑150 (MQ=255)
                                     tGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGc                                                                                                             >  1:2468909/1‑150 (MQ=255)
                                     tGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGc                                                                                                             <  1:4808634/150‑1 (MQ=255)
                                         ctcGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCg                                                                                                         >  2:3839056/1‑150 (MQ=255)
                                              gtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTcc                                                                                                    <  1:1209233/150‑1 (MQ=255)
                                               tAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                   <  2:2459918/150‑1 (MQ=255)
                                                   tAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                               <  1:3941181/150‑1 (MQ=255)
                                                   tAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                               <  2:5147115/150‑1 (MQ=255)
                                                    aaTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                              <  1:3657125/150‑1 (MQ=255)
                                                    aaTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                              >  2:2892413/1‑150 (MQ=255)
                                                      ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggag                                                                                            >  1:2923092/1‑150 (MQ=255)
                                                       tAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                           >  2:1092126/1‑150 (MQ=255)
                                                       tAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                           <  1:2688058/150‑1 (MQ=255)
                                                       tAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCAGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                           >  2:3898852/1‑150 (MQ=255)
                                                       tAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                           >  1:4948259/1‑150 (MQ=255)
                                                         aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGAt                                                                                         >  2:4960009/1‑150 (MQ=255)
                                                         aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGAt                                                                                         <  2:547223/150‑1 (MQ=255)
                                                          cGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATa                                                                                        <  2:848080/150‑1 (MQ=255)
                                                           gatGAGTGGGTTATATTGCGATGTATACGTAGCCGCATCTGCCAAGGGTGTCGGCTTCGCAGTTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaa                                                                                       <  1:1174663/148‑1 (MQ=255)
                                                            cTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaaa                                                                                      <  1:1867283/150‑1 (MQ=255)
                                                                gTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGc                                                                                  <  1:603175/150‑1 (MQ=255)
                                                                     gTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGaa                                                                             <  1:204512/150‑1 (MQ=255)
                                                                       aaaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATg                                                                           <  1:977147/150‑1 (MQ=255)
                                                                        aaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGt                                                                          >  1:2586985/1‑150 (MQ=255)
                                                                          tGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGttt                                                                        >  2:1970620/1‑150 (MQ=255)
                                                                            gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCa                                                                      <  2:2197086/150‑1 (MQ=255)
                                                                            gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCa                                                                      <  2:2957410/150‑1 (MQ=255)
                                                                                  tataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                                >  2:3412945/1‑150 (MQ=255)
                                                                                  tataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                                >  2:1867162/1‑150 (MQ=255)
                                                                                   ataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCt                                                               <  2:4863296/150‑1 (MQ=255)
                                                                                   ataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCt                                                               >  1:1084951/1‑150 (MQ=255)
                                                                                   ataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCt                                                               >  2:1884276/1‑150 (MQ=255)
                                                                                   ataCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCt                                                               >  2:760990/1‑150 (MQ=255)
                                                                                    taCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCtt                                                              <  1:3688623/150‑1 (MQ=255)
                                                                                       gTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCg                                                           <  1:3326353/150‑1 (MQ=255)
                                                                                           ccGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCAt                                                       >  2:563652/1‑150 (MQ=255)
                                                                                            cGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCAtt                                                      <  1:1004767/150‑1 (MQ=255)
                                                                                            cGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCAtt                                                      <  1:758416/150‑1 (MQ=255)
                                                                                              cAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCg                                                    >  2:2356381/1‑150 (MQ=255)
                                                                                                     cAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCAc                                             <  2:1775172/150‑1 (MQ=255)
                                                                                                        gggTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGa                                          <  2:2153098/150‑1 (MQ=255)
                                                                                                         ggTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGaa                                         >  2:4421353/1‑150 (MQ=255)
                                                                                                            gtCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACtt                                      >  2:3038213/1‑150 (MQ=255)
                                                                                                               ggCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGcc                                   >  1:2867498/1‑150 (MQ=255)
                                                                                                                 cTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCCTGGTGGGGGCTGAGTTCACCGGTGGTCAGCCTCGATTCCATCTGGAGGATACAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCgg                                 >  2:823413/1‑150 (MQ=255)
                                                                                                                  ttCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGt                                <  2:2530594/150‑1 (MQ=255)
                                                                                                                   tCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTg                               >  2:3525464/1‑150 (MQ=255)
                                                                                                                      cAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATg                            <  2:5073411/150‑1 (MQ=255)
                                                                                                                        ggTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCt                          >  1:474840/1‑150 (MQ=255)
                                                                                                                         gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTg                         <  2:3195507/150‑1 (MQ=255)
                                                                                                                          ttAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTgg                        >  2:3086099/1‑150 (MQ=255)
                                                                                                                           tAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGc                       <  2:2800152/150‑1 (MQ=255)
                                                                                                                           tAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGc                       >  1:3883323/1‑150 (MQ=255)
                                                                                                                             aaTCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGccc                     <  2:1084951/150‑1 (MQ=255)
                                                                                                                              aTCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCg                    <  2:2779125/150‑1 (MQ=255)
                                                                                                                              aTCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCg                    <  2:1361152/150‑1 (MQ=255)
                                                                                                                                cGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGt                  >  2:2041828/1‑150 (MQ=255)
                                                                                                                                 gTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTc                 >  2:4100289/1‑150 (MQ=255)
                                                                                                                                  ttGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCa                >  2:2895501/1‑150 (MQ=255)
                                                                                                                                       gTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAg           <  1:4114671/150‑1 (MQ=255)
                                                                                                                                            gggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAAt      <  2:2848544/150‑1 (MQ=255)
                                                                                                                                              gTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGc    <  2:1441055/150‑1 (MQ=255)
                                                                                                                                               tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCa   <  2:4137145/150‑1 (MQ=255)
                                                                                                                                               tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCa   <  2:2404505/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAG  >  CP009273/259347‑259640

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: