Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,036,052 C→A 47.4% S152S (TCG→TCT yqfA ← hemolysin III family HylIII inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,036,0520CA47.4% 141.9 / inf 234S152S (TCG→TCTyqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (63/60);  new base A (43/68);  total (106/128)
Fisher's exact test for biased strand distribution p-value = 6.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 239 total aligned reads displayed.

AGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAA  >  CP009273/3035903‑3036201
                                                                                                                                                     |                                                                                                                                                     
    cGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAgagaaa                                                                                                                                                   <  1:2713240/150‑3 (MQ=255)
     gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcag                                                                                                                                                  >  2:1079580/1‑150 (MQ=255)
     gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGAcag                                                                                                                                                  <  2:555912/150‑1 (MQ=255)
         ccGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGAcagccag                                                                                                                                              <  2:1495512/150‑1 (MQ=255)
            tGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                                                                                                                           <  2:2333748/150‑1 (MQ=255)
             gCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCa                                                                                                                                          >  1:4155943/1‑150 (MQ=255)
                aGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCg                                                                                                                                       >  1:280246/1‑150 (MQ=255)
                aGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCg                                                                                                                                       <  1:1035017/150‑1 (MQ=255)
                  aTGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc                                                                                                                                     <  1:956059/150‑1 (MQ=255)
                   tGGCATGGTTGTATGGGATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCa                                                                                                                                    >  2:128077/1‑150 (MQ=255)
                    ggCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCACGAGCTAATATACCACGCCGCCTACCGCCAGTAAGGTAAGGCTGCCCGCCGCGTGCTTAACTGACATTTCATTAATTACCACCAGAGAAATCCAGCCCATCGCCat                                                                                                                                   <  2:2349523/150‑2 (MQ=255)
                      cATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                                                                                                                 <  2:4313352/150‑1 (MQ=255)
                       aTGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATa                                                                                                                                >  2:4523930/1‑150 (MQ=255)
                        tGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAg                                                                                                                               <  2:359124/150‑1 (MQ=255)
                         ggTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg                                                                                                                              <  2:4437246/150‑1 (MQ=255)
                         ggTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg                                                                                                                              <  1:2364214/150‑1 (MQ=255)
                          gTTGTATGGAATGCGTTTGCAGACGTAGACAATCACCCCGAGCGATTAAAGCACGCCGACTACAGCCAGTAAGGTAACGCTGCCCGCCGCGATCTTAACTGCCTTTTCATAATTTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGt                                                                                                                             <  2:3356201/150‑1 (MQ=255)
                               aTGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCa                                                                                                                        <  2:3675163/150‑1 (MQ=255)
                               aTGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCa                                                                                                                        >  1:671724/1‑150 (MQ=255)
                                 ggAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCaga                                                                                                                      >  1:844309/1‑150 (MQ=255)
                                    aTGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGAt                                                                                                                   <  2:1033515/150‑1 (MQ=255)
                                        gTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaaaa                                                                                                               >  1:2294392/1‑150 (MQ=255)
                                         tttGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGCTAAAAt                                                                                                              >  2:1348240/1‑150 (MQ=255)
                                           tGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAttt                                                                                                            >  2:3477322/1‑150 (MQ=255)
                                                aCGTAGAAAATCACCCCGACCGAATACACCACGCCGCCTACAGCAAGAAAGGAAACTCTTCCCTCCGCAAGCTTAACTGCCATATCATAAATTACAACCAGCGACAGCCATCCCAACTCAAGATATTACACCAGAGATAAAATATAGaaa                                                                                                       >  1:1282404/1‑149 (MQ=255)
                                                 cGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATc                                                                                                      >  1:1107527/1‑150 (MQ=255)
                                                 cGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATc                                                                                                      <  1:1453379/150‑1 (MQ=255)
                                                  gTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCg                                                                                                     >  2:2004646/1‑150 (MQ=255)
                                                  gTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTCCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCGAGATAGGTCACCAGAGATAAAATTTTGAATCg                                                                                                     >  2:88414/1‑150 (MQ=255)
                                                     gAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATACAATTTTGAATCtgtg                                                                                                  >  2:399493/1‑150 (MQ=255)
                                                      aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc                                                                                                 <  1:1512388/150‑1 (MQ=255)
                                                      aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc                                                                                                 <  2:4067478/150‑1 (MQ=255)
                                                      aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc                                                                                                 >  2:1760340/1‑150 (MQ=255)
                                                      aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc                                                                                                 <  1:430667/150‑1 (MQ=255)
                                                      aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc                                                                                                 <  1:1596060/150‑1 (MQ=255)
                                                       aaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg                                                                                                >  2:860411/1‑150 (MQ=255)
                                                        aaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCCGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCCGAGATAAAATTTTGAATCGGTgcgc                                                                                               >  2:3538493/1‑150 (MQ=255)
                                                        aaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcgc                                                                                               >  2:772799/1‑150 (MQ=255)
                                                            aCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATg                                                                                           >  2:4237373/1‑150 (MQ=255)
                                                              cccGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGt                                                                                         >  2:414341/1‑150 (MQ=255)
                                                               ccGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTc                                                                                        >  1:1404225/1‑150 (MQ=255)
                                                               ccGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTc                                                                                        <  1:4481610/150‑1 (MQ=255)
                                                                cGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATACAATTTTGAATCGGTGCGCGATGGTCa                                                                                       <  1:2468237/150‑1 (MQ=255)
                                                                 gAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAg                                                                                      >  2:2015064/1‑150 (MQ=255)
                                                                  aGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGt                                                                                     >  2:1200006/1‑150 (MQ=255)
                                                                   gCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGtt                                                                                    <  1:1435590/150‑1 (MQ=255)
                                                                    cGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGttt                                                                                   <  1:3560756/150‑1 (MQ=255)
                                                                      aaTAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTaa                                                                                 <  2:225612/150‑1 (MQ=255)
                                                                      aaTAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCAAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTaa                                                                                 >  2:4614796/1‑150 (MQ=255)
                                                                        tAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAAc                                                                               <  2:2733421/150‑1 (MQ=255)
                                                                         aaaCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACa                                                                              >  2:3253121/1‑150 (MQ=255)
                                                                             cACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAAt                                                                          >  1:3661739/1‑150 (MQ=255)
                                                                               cgccgcCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATAc                                                                        <  1:1932290/150‑1 (MQ=255)
                                                                                gccgccTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATAcc                                                                       <  1:486203/150‑1 (MQ=255)
                                                                                  cgccTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCa                                                                     <  1:2946801/150‑1 (MQ=255)
                                                                                  cgccTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCa                                                                     >  1:1771582/1‑150 (MQ=255)
                                                                                   gccTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATAccatg                                                                    <  1:303026/150‑4 (MQ=255)
                                                                                       aCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAAt                                                                >  1:3567867/1‑150 (MQ=255)
                                                                                       aCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAAt                                                                <  1:3803857/150‑1 (MQ=255)
                                                                                          gCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGcc                                                             <  1:2671275/150‑1 (MQ=255)
                                                                                          gCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATGTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGcc                                                             >  2:2875513/1‑150 (MQ=255)
                                                                                            cAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAg                                                           <  2:3651909/150‑1 (MQ=255)
                                                                                            cAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAg                                                           >  1:1776079/1‑150 (MQ=255)
                                                                                             aGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAgg                                                          <  1:342584/150‑1 (MQ=255)
                                                                                                aaGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTc                                                       >  1:2489564/1‑150 (MQ=255)
                                                                                                 aGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGcgat                                                      <  2:4525800/150‑4 (MQ=255)
                                                                                                 aGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTcc                                                      >  1:463609/1‑150 (MQ=255)
                                                                                                  ggTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCa                                                     >  1:2067302/1‑150 (MQ=255)
                                                                                                  ggTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCa                                                     <  2:2573995/150‑1 (MQ=255)
                                                                                                  ggTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCa                                                     <  1:4016864/150‑1 (MQ=255)
                                                                                                   gTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAg                                                    <  2:3356675/150‑1 (MQ=255)
                                                                                                     aaCGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGAt                                                  <  2:618829/150‑1 (MQ=255)
                                                                                                     aaCGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGAt                                                  >  1:915084/1‑150 (MQ=255)
                                                                                                         cTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACa                                              <  1:1385766/150‑1 (MQ=255)
                                                                                                         cTGCACGCCGCGAGCTTAACTGCCATTTCATAAATTACCAACAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGCGATAAAATTTTGAATCGGTGCTCGATGGTCAGTTTAAACAGAATACCCAGCTATGCCTGGCTCCAGATAACa                                              <  2:509966/150‑1 (MQ=255)
                                                                                                          tGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACaa                                             >  1:4153808/1‑150 (MQ=255)
                                                                                                              cgccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAatcat                                         <  2:841836/150‑1 (MQ=255)
                                                                                                                ccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAAtcatca                                       >  1:1983424/1‑150 (MQ=255)
                                                                                                                ccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAAtcatca                                       <  2:2794220/150‑1 (MQ=255)
                                                                                                                ccgcGAGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAAtcatca                                       <  2:2793652/150‑1 (MQ=255)
                                                                                                                 cgcgAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCaa                                      <  1:669479/150‑1 (MQ=255)
                                                                                                                     aGCTTAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCcg                                  <  2:3936795/150‑1 (MQ=255)
                                                                                                                         tAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCcgcgcg                              <  1:3367271/150‑1 (MQ=255)
                                                                                                                         tAACTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCcgcgcg                              <  1:1131113/150‑1 (MQ=255)
                                                                                                                           aCTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCt                            <  1:32619/150‑1 (MQ=255)
                                                                                                                            cTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTa                           <  2:1508200/150‑1 (MQ=255)
                                                                                                                            cTGCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTa                           <  2:2538242/150‑1 (MQ=255)
                                                                                                                              gCCATTTCATAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAAc                         <  1:553678/150‑1 (MQ=255)
                                                                                                                                cATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACgg                       >  1:3543774/1‑150 (MQ=255)
                                                                                                                                cATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACgg                       >  2:725790/1‑150 (MQ=255)
                                                                                                                                   ttCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGgaga                    <  2:1880851/150‑1 (MQ=255)
                                                                                                                                       tAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATcc                >  2:4570288/1‑150 (MQ=255)
                                                                                                                                       tAAATTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATcc                <  2:1341657/150‑1 (MQ=255)
                                                                                                                                        aaaTTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCa               >  1:1586327/1‑150 (MQ=255)
                                                                                                                                         aaTTACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAAGCCAg              <  2:2240533/150‑1 (MQ=255)
                                                                                                                                           ttACCACCAGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGcc            <  2:3883306/150‑1 (MQ=255)
                                                                                                                                              ccaccaGAGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCcca         >  1:683341/1‑150 (MQ=255)
                                                                                                                                               caccaGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCccac        >  1:4317130/1‑150 (MQ=255)
                                                                                                                                                    gCGACAGCCAGCCGATCGCCAGATAGGTCACCAGAGATAATATTTTGTATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCa   <  1:815769/150‑1 (MQ=255)
                                                                                                                                                     cGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCaa  <  2:3047784/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATGGTCAGTTTAAACAGAATACCCAGCAATGCCAGGCTCCAGATAACAATCATCAACCCGCGCGCTAACGGAGAATCCAGCCCCACCAGCAA  >  CP009273/3035903‑3036201

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: