Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,985 C→T 7.8% E172K (GAA→AAA)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,020,9850CT 7.8% 1044.7 / 55.8 293E172K (GAA→AAA) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (120/150);  new base T (13/10);  total (133/160)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

TACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGC  >  CP009273/2020836‑2021134
                                                                                                                                                     |                                                                                                                                                     
    gcgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAAcc                                                                                                                                                   <  1:4449056/150‑1 (MQ=255)
    gcgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAAcc                                                                                                                                                   <  2:1276606/150‑1 (MQ=255)
      gcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCAt                                                                                                                                                 <  1:3198587/150‑1 (MQ=255)
         gTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATcga                                                                                                                                              >  2:3607941/1‑150 (MQ=255)
         gTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcta                                                                                                                                              >  2:5049777/1‑148 (MQ=255)
         gTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcga                                                                                                                                              >  2:3323678/1‑150 (MQ=255)
         gTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcga                                                                                                                                              <  2:2047140/150‑1 (MQ=255)
            aCGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacga                                                                                                                                           <  2:2306773/150‑1 (MQ=255)
              gATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgagg                                                                                                                                         >  1:130031/1‑150 (MQ=255)
              gATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgagg                                                                                                                                         >  1:2443147/1‑150 (MQ=255)
               aTTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgagga                                                                                                                                        <  2:3465387/150‑1 (MQ=255)
                ttttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACgaggag                                                                                                                                       >  2:1763610/1‑150 (MQ=255)
                ttttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgaggag                                                                                                                                       >  1:4785740/1‑150 (MQ=255)
                  ttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGtt                                                                                                                                     <  2:2678292/150‑1 (MQ=255)
                  ttCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGtt                                                                                                                                     >  2:4048728/1‑150 (MQ=255)
                   tctcGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAAACATCGACGAGGAGTtg                                                                                                                                    >  2:1114664/1‑150 (MQ=255)
                    ctcGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTtgt                                                                                                                                   <  1:2061955/150‑1 (MQ=255)
                      cGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCTGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCTCCTGTTCAACCATCGACGAGGAGTTGTgc                                                                                                                                 >  2:3169225/1‑150 (MQ=255)
                      cGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTgc                                                                                                                                 <  1:4387521/150‑1 (MQ=255)
                          gtgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTa                                                                                                                             >  1:882927/1‑150 (MQ=255)
                           tgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAg                                                                                                                            <  1:1031559/150‑1 (MQ=255)
                           tgAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAg                                                                                                                            <  1:1031608/150‑1 (MQ=255)
                              aTGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCg                                                                                                                         <  2:428627/150‑1 (MQ=255)
                                  ccAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGCACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCgcagc                                                                                                                     >  1:4835934/1‑150 (MQ=255)
                                    aaTAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCaa                                                                                                                   <  2:2183653/150‑1 (MQ=255)
                                     aTAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAg                                                                                                                  <  1:4139109/150‑1 (MQ=255)
                                        aaCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTaa                                                                                                               <  2:2641370/150‑1 (MQ=255)
                                        aaCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTaa                                                                                                               <  2:386252/150‑1 (MQ=255)
                                        aaCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTaa                                                                                                               >  2:4522862/1‑150 (MQ=255)
                                         aCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAAc                                                                                                              <  2:3116226/150‑1 (MQ=255)
                                           cAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCa                                                                                                            >  2:544298/1‑150 (MQ=255)
                                           cAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCa                                                                                                            >  2:3519895/1‑150 (MQ=255)
                                           cAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCa                                                                                                            >  2:3912675/1‑150 (MQ=255)
                                            aGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCaa                                                                                                           <  1:2560878/150‑1 (MQ=255)
                                            aGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCaa                                                                                                           >  2:1605979/1‑150 (MQ=255)
                                             gggTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAAc                                                                                                          <  2:3787529/150‑1 (MQ=255)
                                              ggTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACt                                                                                                         <  2:1383823/150‑1 (MQ=255)
                                               gTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGTTCGGCTGCGTCGAAAAGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACtt                                                                                                        <  1:1847692/150‑1 (MQ=255)
                                                ttGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTc                                                                                                       <  1:2834790/150‑1 (MQ=255)
                                                 tGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTcc                                                                                                      <  1:3558985/150‑1 (MQ=255)
                                                 tGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTcc                                                                                                      <  1:1145576/150‑1 (MQ=255)
                                                  gggTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCa                                                                                                     <  1:3120934/150‑1 (MQ=255)
                                                    gTGACATAGCCGTAATAACGCGATGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGt                                                                                                   >  1:5120919/1‑150 (MQ=255)
                                                         aTAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGc                                                                                              >  1:1600514/1‑150 (MQ=255)
                                                          tAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCg                                                                                             >  1:2775373/1‑150 (MQ=255)
                                                          tAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCg                                                                                             <  1:593745/150‑1 (MQ=255)
                                                             ccGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTa                                                                                          <  1:2408144/150‑1 (MQ=255)
                                                              cGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTaa                                                                                         >  1:4272085/1‑150 (MQ=255)
                                                              cGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTaa                                                                                         <  2:4771383/150‑1 (MQ=255)
                                                               gTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAt                                                                                        <  1:4088454/150‑1 (MQ=255)
                                                                taataaCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAtt                                                                                       >  1:3056536/1‑150 (MQ=255)
                                                                taataaCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAtt                                                                                       >  1:4856702/1‑150 (MQ=255)
                                                                taataaCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAtt                                                                                       >  2:882532/1‑150 (MQ=255)
                                                                 aataaCGCGTTGGAACAGTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAAttt                                                                                      <  1:197261/150‑1 (MQ=255)
                                                                  ataaCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTc                                                                                     <  1:1649079/150‑1 (MQ=255)
                                                                     aCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCACGTAACCAACTTCCAGTGCAGCGGtcatttcatt                                                                                  >  1:4233448/1‑150 (MQ=255)
                                                                     aCGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCAtt                                                                                  >  1:2899835/1‑150 (MQ=255)
                                                                        cgTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTc                                                                               <  2:3431789/150‑1 (MQ=255)
                                                                        cgTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTc                                                                               >  2:4385259/1‑150 (MQ=255)
                                                                        cgTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCACCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTc                                                                               >  2:5115973/1‑150 (MQ=255)
                                                                           tGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTTAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAg                                                                            >  1:3526706/1‑150 (MQ=255)
                                                                            ggAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAgg                                                                           <  1:3684221/150‑1 (MQ=255)
                                                                            ggAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAgg                                                                           <  1:3091441/150‑1 (MQ=255)
                                                                            ggAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAgg                                                                           >  1:2961246/1‑150 (MQ=255)
                                                                             gAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGc                                                                          <  2:907863/150‑1 (MQ=255)
                                                                             gAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGc                                                                          <  1:2322805/150‑1 (MQ=255)
                                                                              aaccgtGCGCGTAAACCAGGTCGGCTGCGTCGAACCGTACCATCCGGCACGCGAGACATCCATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCACGTAACCACCGTCCAGTGCAGCGGTAATTTCATTGTCGAGGCt                                                                         <  2:3547681/145‑1 (MQ=255)
                                                                              aaCATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCt                                                                         <  2:3565000/150‑1 (MQ=255)
                                                                                   tGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCAt                                                                    <  1:3271360/150‑1 (MQ=255)
                                                                                   tGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCAt                                                                    <  1:810573/150‑1 (MQ=255)
                                                                                     cgcgTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCAttt                                                                  <  2:2442534/150‑1 (MQ=255)
                                                                                     cgcgTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCAttt                                                                  >  2:4697868/1‑150 (MQ=255)
                                                                                                 gTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAgg                                                      <  1:4690359/150‑1 (MQ=255)
                                                                                                  tCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCACCTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGga                                                     <  1:1479973/150‑1 (MQ=255)
                                                                                                  tCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGga                                                     >  2:1605662/1‑150 (MQ=255)
                                                                                                      cTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGAtt                                                 <  1:3786978/150‑1 (MQ=255)
                                                                                                        gCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTcc                                               <  1:1236015/150‑1 (MQ=255)
                                                                                                        gCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTcc                                               <  2:120097/150‑1 (MQ=255)
                                                                                                         cGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTccc                                              <  1:1887892/150‑1 (MQ=255)
                                                                                                          gTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCCGCGGTAATTGCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTcccc                                             >  2:2747008/1‑150 (MQ=255)
                                                                                                          gTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTcccc                                             <  1:4274894/150‑1 (MQ=255)
                                                                                                            cGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTc                                           >  2:3221986/1‑150 (MQ=255)
                                                                                                            cGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTc                                           <  1:4013684/150‑1 (MQ=255)
                                                                                                             gAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCg                                          >  2:4671186/1‑150 (MQ=255)
                                                                                                                aCGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTg                                       <  1:1891920/150‑1 (MQ=255)
                                                                                                                aCGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTg                                       <  1:2653875/150‑1 (MQ=255)
                                                                                                                  gTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGAc                                     >  1:3931196/1‑150 (MQ=255)
                                                                                                                     aaaTCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTa                                   >  1:736194/1‑149 (MQ=255)
                                                                                                                      aaTCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAAt                                 >  2:2996530/1‑150 (MQ=255)
                                                                                                                                         acatCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTaa              >  1:2573005/1‑150 (MQ=255)
                                                                                                                                           atCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACt            <  2:1073649/150‑1 (MQ=255)
                                                                                                                                           atCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACt            <  2:669087/150‑1 (MQ=255)
                                                                                                                                           atCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACt            >  1:451397/1‑150 (MQ=255)
                                                                                                                                            tCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTg           >  1:2278620/1‑150 (MQ=255)
                                                                                                                                            tCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTg           >  1:2351980/1‑150 (MQ=255)
                                                                                                                                             cGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTgg          <  1:1421665/150‑1 (MQ=255)
                                                                                                                                              gCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGc         >  2:2472245/1‑150 (MQ=255)
                                                                                                                                                   ttCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGt    <  2:2842212/150‑1 (MQ=255)
                                                                                                                                                    tCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTg   <  2:2906185/150‑1 (MQ=255)
                                                                                                                                                     cAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGc  >  1:2048641/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGTGCAGCGGTAATTTCATTGTCGAGGCTTTCATTTTCGCGAGACAGGAGATTCCCCTCGCTGACTAATGCATCCGACACACCATTAACTGGCAGAGTGC  >  CP009273/2020836‑2021134

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: