Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,166 C→A 65.7% 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,1660CA65.7% 207.7 / 264.1 216D445Y (GAC→TAC) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base A (60/82);  total (97/119)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 220 total aligned reads displayed.

CGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACT  >  CP009273/2020017‑2020315
                                                                                                                                                     |                                                                                                                                                     
  cAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCAt                                                                                                                                                     >  2:3620120/1‑150 (MQ=255)
   aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCAtt                                                                                                                                                    >  1:1166460/1‑150 (MQ=255)
    gCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTa                                                                                                                                                   <  1:1565057/150‑1 (MQ=255)
      cTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAAt                                                                                                                                                 >  1:4700112/1‑150 (MQ=255)
       tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATc                                                                                                                                                <  1:624299/150‑1 (MQ=255)
       tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATc                                                                                                                                                <  1:325967/150‑1 (MQ=255)
         cGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGc                                                                                                                                              >  1:840647/1‑150 (MQ=255)
          gTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCt                                                                                                                                             <  2:4429826/150‑1 (MQ=255)
            cAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCttt                                                                                                                                           <  2:394887/150‑1 (MQ=255)
               aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaaa                                                                                                                                        >  2:4659958/1‑150 (MQ=255)
               aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaaa                                                                                                                                        >  2:4660046/1‑150 (MQ=255)
                cTCGCGCCTGGCAGAATCACACAACACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaaaa                                                                                                                                       <  2:3415311/150‑1 (MQ=255)
                cTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaaaa                                                                                                                                       <  2:3137938/150‑1 (MQ=255)
                 tCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAAt                                                                                                                                      >  2:4543899/1‑150 (MQ=255)
                 tCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAAt                                                                                                                                      <  1:514448/150‑1 (MQ=255)
                  cgcgCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATg                                                                                                                                     >  2:416790/1‑150 (MQ=255)
                  cgcgCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATg                                                                                                                                     <  1:961/150‑1 (MQ=255)
                   gcgcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATgg                                                                                                                                    <  2:2746356/150‑1 (MQ=255)
                      ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCa                                                                                                                                 >  2:4600689/1‑150 (MQ=255)
                       cTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCaa                                                                                                                                >  2:425158/1‑150 (MQ=255)
                       cTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCaa                                                                                                                                <  1:4554832/150‑1 (MQ=255)
                        tGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAg                                                                                                                               >  1:4020607/1‑150 (MQ=255)
                         ggCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGa                                                                                                                              <  2:4107213/150‑1 (MQ=255)
                          gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                             >  2:943474/1‑150 (MQ=255)
                            aGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCg                                                                                                                           <  2:3674294/150‑1 (MQ=255)
                             gAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGa                                                                                                                          <  1:3645383/150‑1 (MQ=255)
                              aaTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTGATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGAc                                                                                                                         >  2:286020/1‑150 (MQ=255)
                              aaTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGAc                                                                                                                         <  1:940073/150‑1 (MQ=255)
                              aaTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGAc                                                                                                                         >  2:4787861/1‑150 (MQ=255)
                                tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCt                                                                                                                       <  2:943966/150‑1 (MQ=255)
                                    acaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAAt                                                                                                                   <  2:1132314/150‑1 (MQ=255)
                                     caAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATg                                                                                                                  >  2:1098244/1‑150 (MQ=255)
                                     caAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATg                                                                                                                  >  2:3966833/1‑150 (MQ=255)
                                        aaCTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACa                                                                                                               >  2:3511717/1‑150 (MQ=255)
                                          ctcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGa                                                                                                             >  2:2401811/1‑150 (MQ=255)
                                             ctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaa                                                                                                          <  1:231726/150‑1 (MQ=255)
                                               caccacCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                        <  2:403206/150‑1 (MQ=255)
                                               caccacCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                        <  2:1227582/150‑1 (MQ=255)
                                                accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAgg                                                                                                       >  2:824334/1‑150 (MQ=255)
                                                  caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                                     <  2:626453/150‑1 (MQ=255)
                                                  caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACATGAATGACAGAAAAAGGAt                                                                                                     >  1:2820736/1‑150 (MQ=255)
                                                  caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                                     >  2:4109712/1‑150 (MQ=255)
                                                   accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                    <  1:3926420/150‑1 (MQ=255)
                                                   accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                    <  2:3230112/150‑1 (MQ=255)
                                                   accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                    <  1:2251873/150‑1 (MQ=255)
                                                     cGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGtt                                                                                                  <  2:3365987/150‑1 (MQ=255)
                                                     cGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGtt                                                                                                  >  1:3168159/1‑150 (MQ=255)
                                                      gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                                 <  1:952788/150‑1 (MQ=255)
                                                        cccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGt                                                                                               <  1:62064/150‑1 (MQ=255)
                                                          cgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAACAAGGATGTTGGTGc                                                                                             >  2:2653847/1‑150 (MQ=255)
                                                            cccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGt                                                                                           <  2:4339648/150‑1 (MQ=255)
                                                             ccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTc                                                                                          <  2:1835341/150‑1 (MQ=255)
                                                                 cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                      >  2:4288063/1‑150 (MQ=255)
                                                                 cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                      >  2:3081717/1‑150 (MQ=255)
                                                                 cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                      <  2:1789617/150‑1 (MQ=255)
                                                                  aaCGTCCTGCGCCCGCAAGGAACTGCTAATTAAGCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGaca                                                                                     <  1:2518425/150‑1 (MQ=255)
                                                                   aCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGgcac                                                                                    >  2:1196762/1‑150 (MQ=255)
                                                                            gcacgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcg                                                                           <  2:3276702/150‑1 (MQ=255)
                                                                             cacgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcga                                                                          <  1:1920626/150‑1 (MQ=255)
                                                                              acgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgag                                                                         >  1:383858/1‑150 (MQ=255)
                                                                                gcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagcg                                                                       >  2:2522912/1‑150 (MQ=255)
                                                                                gcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagcg                                                                       <  1:1639306/150‑1 (MQ=255)
                                                                                 caAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCgg                                                                      >  1:2406776/1‑150 (MQ=255)
                                                                                  aaGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACTCAATTTAGCGAGCGGa                                                                     <  2:788979/150‑1 (MQ=255)
                                                                                      aaCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACggg                                                                 >  2:2104980/1‑150 (MQ=255)
                                                                                        cTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAGGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCt                                                               >  2:1325443/1‑150 (MQ=255)
                                                                                        cTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCt                                                               <  2:4488529/150‑1 (MQ=255)
                                                                                        cTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCt                                                               >  2:2157412/1‑150 (MQ=255)
                                                                                           cTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtttt                                                            >  1:4218397/1‑150 (MQ=255)
                                                                                             aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                          >  2:597211/1‑150 (MQ=255)
                                                                                             aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                          <  2:498014/150‑1 (MQ=255)
                                                                                             aattaatCCGGCAGCATGAGAAAGAACACGGTCGCACGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTAAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                          <  1:1293718/150‑1 (MQ=255)
                                                                                                taatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCAGACACAATTTAGCGAGCGGACGGGCTTTTTCGaa                                                       >  1:3724010/1‑150 (MQ=255)
                                                                                                  atCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACa                                                     <  2:4645039/150‑1 (MQ=255)
                                                                                                   tCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAg                                                    <  2:3162305/150‑1 (MQ=255)
                                                                                                    ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                                   <  2:3163959/150‑1 (MQ=255)
                                                                                                    ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                                   <  1:3624529/150‑1 (MQ=255)
                                                                                                     cGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTg                                                  <  1:3069045/150‑1 (MQ=255)
                                                                                                      ggCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgc                                                 >  2:905586/1‑150 (MQ=255)
                                                                                                          gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCa                                             <  1:4034889/150‑1 (MQ=255)
                                                                                                          gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCa                                             >  2:757780/1‑150 (MQ=255)
                                                                                                           caTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCAc                                            >  2:3710921/1‑150 (MQ=255)
                                                                                                           caTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCAc                                            >  2:2968938/1‑150 (MQ=255)
                                                                                                               agaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGtt                                        <  1:3329034/150‑1 (MQ=255)
                                                                                                                gacagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta                                       <  2:3191725/147‑1 (MQ=255)
                                                                                                                gaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta                                       <  1:2137914/150‑1 (MQ=255)
                                                                                                                   agaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatat                                    <  1:2912918/150‑1 (MQ=255)
                                                                                                                       cacGGTCGCCCGCCTGATGACCAAAGAGGTAATTAATCTCTTTAAAATGGTCAAGATCGACCTGAATGACAGTAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGTTCTTTTTCTAACAGTGCGCCACGGTTTTATAAAc                                <  1:1269124/150‑1 (MQ=255)
                                                                                                                          ggTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcgcg                             >  1:2489789/1‑150 (MQ=255)
                                                                                                                           gTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGt                            <  1:4324868/150‑1 (MQ=255)
                                                                                                                                 cgccTGATGACCAAAGCGTTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgt                      <  2:916980/150‑1 (MQ=255)
                                                                                                                                 cgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgt                      <  1:1098244/150‑1 (MQ=255)
                                                                                                                                  gccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtg                     <  2:3945068/150‑1 (MQ=255)
                                                                                                                                    cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc                   >  1:1025782/1‑150 (MQ=255)
                                                                                                                                    cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc                   <  1:4116800/150‑1 (MQ=255)
                                                                                                                                      gatgaACAAAGCGTTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGt                 <  2:2715149/150‑1 (MQ=255)
                                                                                                                                         gaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc              <  2:4710635/150‑1 (MQ=255)
                                                                                                                                               aGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCt        <  2:1924137/150‑1 (MQ=255)
                                                                                                                                                 cGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGc      >  1:1845253/1‑150 (MQ=255)
                                                                                                                                                     cATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAATGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACt  <  2:3005999/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCACT  >  CP009273/2020017‑2020315

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: