Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,166 C→A 93.3% D445Y (GAC→TAC)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,020,1660CA93.3% 822.4 / 45.3 272D445Y (GAC→TAC) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (7/12);  new base A (109/144);  total (116/156)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 274 total aligned reads displayed.

AGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCT  >  CP009273/2020020‑2020309
                                                                                                                                                  |                                                                                                                                               
aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAtt                                                                                                                                              >  2:3476730/1‑150 (MQ=255)
aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAtt                                                                                                                                              <  1:330458/150‑1 (MQ=255)
  ccTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGACCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTaa                                                                                                                                            >  2:4116555/1‑150 (MQ=255)
   cTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGtaattaat                                                                                                                                           >  1:1107123/1‑150 (MQ=255)
       gTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCt                                                                                                                                       >  1:2556587/1‑150 (MQ=255)
          aGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTa                                                                                                                                    <  2:797696/150‑1 (MQ=255)
             cTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTaaaa                                                                                                                                 <  2:305628/150‑1 (MQ=255)
                gcgcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATgg                                                                                                                              >  2:3905445/1‑150 (MQ=255)
                 cgcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGt                                                                                                                             <  1:281603/150‑1 (MQ=255)
                  gcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGACAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTc                                                                                                                            >  2:1133392/1‑150 (MQ=255)
                   ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCa                                                                                                                           >  1:2368944/1‑150 (MQ=255)
                    cTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCaa                                                                                                                          >  1:2491427/1‑150 (MQ=255)
                      ggCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGa                                                                                                                        <  1:1463686/150‑1 (MQ=255)
                       gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                       <  2:4232306/150‑1 (MQ=255)
                       gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                       <  1:2797024/150‑1 (MQ=255)
                        cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATc                                                                                                                      <  2:1510230/150‑1 (MQ=255)
                        cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATc                                                                                                                      >  1:3934574/1‑150 (MQ=255)
                        cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATc                                                                                                                      <  2:5023124/150‑1 (MQ=255)
                        cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATc                                                                                                                      <  2:1177055/150‑1 (MQ=255)
                                cacaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGaa                                                                                                              <  1:1178599/150‑1 (MQ=255)
                                cacaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGaa                                                                                                              <  1:4750376/150‑1 (MQ=255)
                                         cctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaa                                                                                                     >  2:1590930/1‑150 (MQ=255)
                                         cctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTACTTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaa                                                                                                     >  2:4072149/1‑150 (MQ=255)
                                          ctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTACTTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaa                                                                                                    >  1:4622750/1‑150 (MQ=255)
                                           tcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                                                   >  1:2447545/1‑150 (MQ=255)
                                           tcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                                                   <  2:4634438/150‑1 (MQ=255)
                                           tcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                                                   <  2:2079375/150‑1 (MQ=255)
                                           tcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                                                   >  1:1233223/1‑150 (MQ=255)
                                               caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                               <  1:2181178/150‑1 (MQ=255)
                                               caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                               <  1:2181096/150‑1 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                              >  1:1588698/1‑150 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                              >  1:4261034/1‑150 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                              <  1:1001307/150‑1 (MQ=255)
                                                   gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGTCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                           >  2:540890/1‑150 (MQ=255)
                                                     cccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGt                                                                                         <  1:659975/150‑1 (MQ=255)
                                                      ccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTg                                                                                        <  1:601667/150‑1 (MQ=255)
                                                       cgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGc                                                                                       <  1:2911367/150‑1 (MQ=255)
                                                        gcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGATTGACAGAAAAAGGATGTTGGTGCg                                                                                      <  2:3381150/150‑1 (MQ=255)
                                                         cccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGt                                                                                     >  1:1706777/1‑150 (MQ=255)
                                                             gCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGa                                                                                 >  2:4176071/1‑150 (MQ=255)
                                                              cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                <  2:4362062/150‑1 (MQ=255)
                                                              cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                >  2:1101868/1‑150 (MQ=255)
                                                              cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                >  2:1087565/1‑150 (MQ=255)
                                                               aaCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGaca                                                                               <  1:1598467/150‑1 (MQ=255)
                                                               aaCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGaca                                                                               >  1:3110816/1‑150 (MQ=255)
                                                                 cGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGAcaca                                                                             >  1:4548441/1‑150 (MQ=255)
                                                                   tCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAt                                                                           <  1:4898677/150‑1 (MQ=255)
                                                                     cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAttt                                                                         <  2:1503990/150‑1 (MQ=255)
                                                                         gcacgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcg                                                                     <  2:500138/150‑1 (MQ=255)
                                                                                aGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGAc                                                              <  1:3853712/150‑1 (MQ=255)
                                                                                aGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGAc                                                              <  2:4864565/150‑1 (MQ=255)
                                                                                 ggAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACg                                                             >  1:4463086/1‑150 (MQ=255)
                                                                                  gAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACgg                                                            >  1:787326/1‑150 (MQ=255)
                                                                                   aaCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACggg                                                           <  1:209328/150‑1 (MQ=255)
                                                                                      tGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtt                                                        <  1:1871904/150‑1 (MQ=255)
                                                                                         taattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttttt                                                     <  1:4427948/150‑1 (MQ=255)
                                                                                         taattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttttt                                                     <  2:4437695/150‑1 (MQ=255)
                                                                                          aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCACAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                    <  2:2429337/150‑1 (MQ=255)
                                                                                          aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                    <  2:5113626/150‑1 (MQ=255)
                                                                                          aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                    <  1:5065562/150‑1 (MQ=255)
                                                                                          aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                    <  2:2192211/150‑1 (MQ=255)
                                                                                            ttaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGa                                                  <  2:49390/150‑1 (MQ=255)
                                                                                            ttaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGa                                                  <  2:2253404/150‑1 (MQ=255)
                                                                                                 ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                             <  1:981312/150‑1 (MQ=255)
                                                                                                 ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                             >  2:3475242/1‑150 (MQ=255)
                                                                                                    gcagcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgcg                                          <  1:3744396/150‑1 (MQ=255)
                                                                                                     cagcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgcgc                                         <  2:140337/150‑1 (MQ=255)
                                                                                                      agcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGcc                                        >  2:1702748/1‑150 (MQ=255)
                                                                                                         aTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACg                                     >  2:2655/1‑150 (MQ=255)
                                                                                                          tGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACgg                                    <  1:4100359/150‑1 (MQ=255)
                                                                                                           gagaAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                   >  2:4314236/1‑150 (MQ=255)
                                                                                                           gagaAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                   <  1:3591491/150‑1 (MQ=255)
                                                                                                              aaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtat                                <  2:1094231/150‑1 (MQ=255)
                                                                                                              aaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtat                                >  1:1315820/1‑150 (MQ=255)
                                                                                                                agaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatat                              >  2:759752/1‑150 (MQ=255)
                                                                                                                    cacGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc                          <  2:4456045/150‑1 (MQ=255)
                                                                                                                     acGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcg                         >  2:1433853/1‑150 (MQ=255)
                                                                                                                     acGGTCGCCCGCCTGATGACCAAAGCGGGAAGTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGCACAGTGCGCCACGGGtctatccacg                         >  2:1281312/1‑145 (MQ=255)
                                                                                                                       ggTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcgcg                       >  2:3192483/1‑150 (MQ=255)
                                                                                                                           gcccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTaa                   <  1:4151113/150‑1 (MQ=255)
                                                                                                                             ccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTCTATAAACGCGTGAAgg                 >  2:297109/1‑150 (MQ=255)
                                                                                                                              cgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgt                <  1:3190998/150‑1 (MQ=255)
                                                                                                                               gccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtg               >  2:392788/1‑150 (MQ=255)
                                                                                                                                ccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtgt              <  2:54870/150‑1 (MQ=255)
                                                                                                                                 cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGACAAAGGAGGTTGGTGCGTCTGACACACTTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTGCTATCAACGCGTGAAGGAGTc             >  2:3406061/1‑150 (MQ=255)
                                                                                                                                 cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc             >  1:4359709/1‑150 (MQ=255)
                                                                                                                                 cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc             >  2:3268919/1‑150 (MQ=255)
                                                                                                                                   gatgaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGt           >  2:2811373/1‑150 (MQ=255)
                                                                                                                                   gatgaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGt           <  2:1930636/150‑1 (MQ=255)
                                                                                                                                   gatgaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGt           <  2:344502/150‑1 (MQ=255)
                                                                                                                                    atgaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTg          <  2:3222975/150‑1 (MQ=255)
                                                                                                                                     tgaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGc         <  2:1059758/150‑1 (MQ=255)
                                                                                                                                      gaCCAACGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc        <  1:3114786/150‑1 (MQ=255)
                                                                                                                                      gaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc        <  1:2268892/150‑1 (MQ=255)
                                                                                                                                      gaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc        <  2:3168582/150‑1 (MQ=255)
                                                                                                                                       aCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCa       >  1:3025643/1‑150 (MQ=255)
                                                                                                                                        ccAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCAc      <  2:257662/150‑1 (MQ=255)
                                                                                                                                        ccAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCAc      <  1:4824244/150‑1 (MQ=255)
                                                                                                                                         cAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACg     >  2:3793159/1‑150 (MQ=255)
                                                                                                                                            aGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCt  <  2:4752492/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCT  >  CP009273/2020020‑2020309

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: