Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,285,476 C→A 10.5% I259I (ATC→ATA yraK → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,285,4760CA10.5% 848.1 / 51.6 238I259I (ATC→ATAyraKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (94/119);  new base A (12/13);  total (106/132)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 244 total aligned reads displayed.

TTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCT  >  CP009273/3285327‑3285625
                                                                                                                                                     |                                                                                                                                                     
ttATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATc                                                                                                                                                       >  2:2785071/1‑150 (MQ=255)
   tCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACt                                                                                                                                                    <  1:2315539/150‑1 (MQ=255)
    cGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACtt                                                                                                                                                   <  1:4028692/150‑1 (MQ=255)
    cGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACtt                                                                                                                                                   <  1:1517871/150‑1 (MQ=255)
    cGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACtt                                                                                                                                                   >  2:3401735/1‑150 (MQ=255)
     gATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTa                                                                                                                                                  <  2:472473/150‑1 (MQ=255)
     gATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTa                                                                                                                                                  >  1:981035/1‑150 (MQ=255)
       ttCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATg                                                                                                                                                <  2:1380409/150‑1 (MQ=255)
        tCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTAGTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGa                                                                                                                                               >  1:2881719/1‑150 (MQ=255)
         cAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGAc                                                                                                                                              >  2:3865961/1‑150 (MQ=255)
         cAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGAc                                                                                                                                              >  1:88828/1‑150 (MQ=255)
          aGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACt                                                                                                                                             <  1:3779529/150‑1 (MQ=255)
          aGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACt                                                                                                                                             <  1:1680554/150‑1 (MQ=255)
            ggAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACtgt                                                                                                                                           <  1:1975405/150‑1 (MQ=255)
             gAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACtgtg                                                                                                                                          <  2:3449522/150‑1 (MQ=255)
              aaaTCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTgtgg                                                                                                                                         >  2:174278/1‑150 (MQ=255)
               aaTCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTgtggt                                                                                                                                        <  1:862605/150‑1 (MQ=255)
                  ctcACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTAACACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGAt                                                                                                                                     >  1:85440/1‑150 (MQ=255)
                   tcACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATa                                                                                                                                    >  1:3447118/1‑150 (MQ=255)
                     aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACt                                                                                                                                  >  1:2113166/1‑150 (MQ=255)
                     aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACt                                                                                                                                  <  1:1838555/150‑1 (MQ=255)
                     aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACt                                                                                                                                  >  2:1062310/1‑150 (MQ=255)
                     aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACt                                                                                                                                  >  2:4091969/1‑150 (MQ=255)
                       cGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTc                                                                                                                                <  1:2944622/150‑1 (MQ=255)
                       cGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTc                                                                                                                                >  2:3172210/1‑150 (MQ=255)
                               aGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTa                                                                                                                        >  2:847794/1‑150 (MQ=255)
                                gTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGATGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaa                                                                                                                       <  2:3557033/150‑1 (MQ=255)
                                  cGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaaaa                                                                                                                     <  2:2691024/150‑1 (MQ=255)
                                   gTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaaaaa                                                                                                                    >  1:2979844/1‑150 (MQ=255)
                                   gTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaaaaa                                                                                                                    >  1:1743182/1‑150 (MQ=255)
                                    tGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAAc                                                                                                                   <  2:2649706/150‑1 (MQ=255)
                                        aaaTTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGt                                                                                                               <  2:740308/150‑1 (MQ=255)
                                        aaaTTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGt                                                                                                               >  1:2716869/1‑150 (MQ=255)
                                          aTTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTg                                                                                                             <  2:1401466/150‑1 (MQ=255)
                                          aTTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTg                                                                                                             >  2:1512699/1‑150 (MQ=255)
                                           ttAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCAAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAAATCGTTGc                                                                                                            >  2:2108813/1‑150 (MQ=255)
                                           ttAATCAGGGAGATTTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGc                                                                                                            <  2:1390522/150‑1 (MQ=255)
                                             taTCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCaa                                                                                                          <  2:200095/149‑1 (MQ=255)
                                             aaTCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCaa                                                                                                          <  2:1762983/150‑1 (MQ=255)
                                               tCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAAt                                                                                                        >  2:3522652/1‑150 (MQ=255)
                                               tCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAAt                                                                                                        >  2:1869151/1‑150 (MQ=255)
                                                cAGGGCGATGTTATTTAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCTCCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATg                                                                                                       <  1:3532553/150‑1 (MQ=255)
                                                 aGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATgc                                                                                                      <  2:3121024/150‑1 (MQ=255)
                                                  gggCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATgcg                                                                                                     <  2:4141064/150‑1 (MQ=255)
                                                   ggCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATgcgc                                                                                                    >  2:2159986/1‑150 (MQ=255)
                                                    gCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCa                                                                                                   <  2:72888/150‑1 (MQ=255)
                                                               aaGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACa                                                                                        <  2:3189656/150‑1 (MQ=255)
                                                               aaGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACa                                                                                        >  2:3284625/1‑150 (MQ=255)
                                                                 ggTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAAc                                                                                      <  2:2525898/150‑1 (MQ=255)
                                                                  gTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACt                                                                                     <  2:3548917/150‑1 (MQ=255)
                                                                         ttGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGgtag                                                                              >  1:1276958/1‑150 (MQ=255)
                                                                          tGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGgtagt                                                                             >  1:650903/1‑150 (MQ=255)
                                                                                catcaATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCa                                                                       >  2:2377201/1‑150 (MQ=255)
                                                                                    aaTGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTAc                                                                   >  1:1190392/1‑150 (MQ=255)
                                                                                      tGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACaa                                                                 <  1:3042209/150‑1 (MQ=255)
                                                                                        gTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAAcc                                                               <  1:544674/150‑1 (MQ=255)
                                                                                          cAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTg                                                             <  1:1403652/150‑1 (MQ=255)
                                                                                           aGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTgg                                                            >  2:2397733/1‑150 (MQ=255)
                                                                                           aGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTgg                                                            >  2:3958452/1‑150 (MQ=255)
                                                                                            gAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGt                                                           >  1:3065455/1‑150 (MQ=255)
                                                                                             aaGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTc                                                          <  2:942282/150‑1 (MQ=255)
                                                                                             aaGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTc                                                          >  1:3304307/1‑150 (MQ=255)
                                                                                                tttACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGcc                                                       <  2:1273957/150‑1 (MQ=255)
                                                                                                 ttACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCa                                                      <  1:1811017/150‑1 (MQ=255)
                                                                                                  tACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAg                                                     <  1:4157818/150‑1 (MQ=255)
                                                                                                  tACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAg                                                     <  1:2231075/150‑1 (MQ=255)
                                                                                                    ccaccCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGc                                                   >  1:3961711/1‑150 (MQ=255)
                                                                                                     caccCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCg                                                  >  2:308619/1‑150 (MQ=255)
                                                                                                     caccCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCg                                                  >  2:1151699/1‑150 (MQ=255)
                                                                                                     caccCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCg                                                  <  2:2625349/150‑1 (MQ=255)
                                                                                                      accCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCga                                                 <  2:1369790/150‑1 (MQ=255)
                                                                                                       cccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCgaa                                                >  2:1379393/1‑150 (MQ=255)
                                                                                                        ccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCgaag                                               <  2:1461419/150‑1 (MQ=255)
                                                                                                        ccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCgaag                                               <  1:218948/150‑1 (MQ=255)
                                                                                                        ccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCgaag                                               <  1:912805/150‑1 (MQ=255)
                                                                                                          gCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCgaagaa                                             >  2:3048425/1‑150 (MQ=255)
                                                                                                           cAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTTATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGaagaag                                            >  1:2393027/1‑150 (MQ=255)
                                                                                                           cAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGaagaag                                            <  1:454292/150‑1 (MQ=255)
                                                                                                            aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGt                                           >  1:2616446/1‑150 (MQ=255)
                                                                                                            aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGt                                           <  2:3652730/150‑1 (MQ=255)
                                                                                                            aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCACATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGt                                           <  1:3868480/150‑1 (MQ=255)
                                                                                                              tgccatgccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTc                                         <  1:2370381/150‑1 (MQ=255)
                                                                                                              tgccatgccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTc                                         <  1:3765701/150‑1 (MQ=255)
                                                                                                                      caGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtg                                 >  1:1900556/1‑150 (MQ=255)
                                                                                                                       aGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAGGAAGTTCAGACAAtgt                                >  2:605865/1‑150 (MQ=255)
                                                                                                                       aGACGGTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtgt                                >  2:2457136/1‑150 (MQ=255)
                                                                                                                           ggTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGccc                            <  2:1213068/150‑1 (MQ=255)
                                                                                                                           ggTTTTACTCCAGTAGACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGccc                            >  2:1704634/1‑150 (MQ=255)
                                                                                                                             ttttACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGa                          <  2:254970/150‑1 (MQ=255)
                                                                                                                                 aCTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTc                      <  2:4016059/150‑1 (MQ=255)
                                                                                                                                     cAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTa                  >  1:3280772/1‑150 (MQ=255)
                                                                                                                                      aGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTAt                 <  2:390872/150‑1 (MQ=255)
                                                                                                                                       gTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTAtt                <  2:591960/150‑1 (MQ=255)
                                                                                                                                        tAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGGAGTTCAGACAATGTGCCCGATGTCGGTATTc               >  2:2740169/1‑150 (MQ=255)
                                                                                                                                          gACTTTGATATAACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGt             >  1:1162493/1‑150 (MQ=255)
                                                                                                                                            cTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTAt           >  2:2852340/1‑150 (MQ=255)
                                                                                                                                               tGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGa        <  1:3142604/150‑1 (MQ=255)
                                                                                                                                                gATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaa       >  1:1000039/1‑150 (MQ=255)
                                                                                                                                                  tatCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaaa     >  1:3208856/1‑150 (MQ=255)
                                                                                                                                                     cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtct  >  1:1434313/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCT  >  CP009273/3285327‑3285625

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: