Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 300,770 C→A 6.1% V49V (GTG→GTT ecpE ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273300,7700CA 6.1% 1129.9 / 31.4 295V49V (GTG→GTTecpEECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (155/122);  new base A (8/10);  total (163/132)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 299 total aligned reads displayed.

AAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGA  >  CP009273/300621‑300918
                                                                                                                                                     |                                                                                                                                                    
aaaTACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTc                                                                                                                                                      >  1:1311091/1‑150 (MQ=255)
 aaTACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCa                                                                                                                                                     >  1:3827565/1‑150 (MQ=255)
  aTACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCAc                                                                                                                                                    >  1:4154811/1‑150 (MQ=255)
  aTACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAAc                                                                                                                                                    >  1:4532574/1‑150 (MQ=255)
   tACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACt                                                                                                                                                   <  1:3788821/150‑1 (MQ=255)
     cTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTac                                                                                                                                                 >  2:476236/1‑150 (MQ=255)
         ctctCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGt                                                                                                                                             <  2:1362093/150‑1 (MQ=255)
          tctcACCATCCTGCAACGCCAGCTGGCGGGTGGCGAAAATCAGTTCACCATCCACCGGTCGGGTCCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTATTGTTATTCACTACACGtt                                                                                                                                            >  1:3692589/1‑150 (MQ=255)
               ccAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTg                                                                                                                                       <  1:180413/150‑1 (MQ=255)
                cAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGa                                                                                                                                      >  2:1694063/1‑150 (MQ=255)
                cAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGa                                                                                                                                      >  1:4027612/1‑150 (MQ=255)
                cAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGa                                                                                                                                      >  1:2379626/1‑150 (MQ=255)
                cAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGAGCGCTTTTGTTGTTATTCACTACACGTTTGCTGa                                                                                                                                      >  2:2532409/1‑150 (MQ=255)
                 aGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGAc                                                                                                                                     >  1:590540/1‑150 (MQ=255)
                  gCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACa                                                                                                                                    >  2:3388197/1‑150 (MQ=255)
                    cTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACaaa                                                                                                                                  <  2:938243/150‑1 (MQ=255)
                     tGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAag                                                                                                                                 <  1:2885052/150‑1 (MQ=255)
                       cAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTGATTCACTACACGTTTGCTGACAAagtc                                                                                                                               <  2:2321741/150‑1 (MQ=255)
                        aaCGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAagtca                                                                                                                              <  1:1966183/150‑1 (MQ=255)
                          cGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAagtcagt                                                                                                                            >  1:1390038/1‑150 (MQ=255)
                          cGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAagtcagt                                                                                                                            >  2:1098826/1‑150 (MQ=255)
                              aGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCg                                                                                                                        >  1:2755583/1‑150 (MQ=255)
                              aGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCg                                                                                                                        >  2:3200582/1‑150 (MQ=255)
                               gCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCgg                                                                                                                       <  2:2012907/150‑1 (MQ=255)
                                 tGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGAc                                                                                                                     >  2:3618081/1‑150 (MQ=255)
                                  ggCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACg                                                                                                                    >  2:412311/1‑150 (MQ=255)
                                        gggCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCg                                                                                                              >  2:2467468/1‑150 (MQ=255)
                                        gggCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCg                                                                                                              >  2:1811308/1‑150 (MQ=255)
                                            gAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGaaaa                                                                                                          >  1:1475806/1‑150 (MQ=255)
                                                 gCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCa                                                                                                     <  1:713749/150‑1 (MQ=255)
                                                    gTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGAt                                                                                                  <  1:2865230/150‑1 (MQ=255)
                                                         ccATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAATATGTCAGATTGccg                                                                                             <  1:1560162/150‑1 (MQ=255)
                                                              cACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGccgaccga                                                                                        <  1:3567816/150‑1 (MQ=255)
                                                              cACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGccgaccga                                                                                        >  2:1759583/1‑150 (MQ=255)
                                                              cACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGccgaccga                                                                                        <  2:1828627/150‑1 (MQ=255)
                                                               aCCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaa                                                                                       >  2:580159/1‑150 (MQ=255)
                                                               aCCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaa                                                                                       >  2:4457177/1‑150 (MQ=255)
                                                                  ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                    <  1:3471102/150‑1 (MQ=255)
                                                                  ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                    <  1:259669/150‑1 (MQ=255)
                                                                  ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                    >  2:3501427/1‑150 (MQ=255)
                                                                        gTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCtgt                                                                              >  2:2244707/1‑150 (MQ=255)
                                                                           cgcAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTc                                                                           <  1:3232158/150‑1 (MQ=255)
                                                                            gcAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCa                                                                          >  2:4104121/1‑150 (MQ=255)
                                                                            gcAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCa                                                                          <  1:3936765/150‑1 (MQ=255)
                                                                              aaTTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCaaa                                                                        >  1:3325836/1‑150 (MQ=255)
                                                                               aTTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAg                                                                       >  1:4069899/1‑150 (MQ=255)
                                                                               aTTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTAAACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAg                                                                       <  2:376487/150‑1 (MQ=255)
                                                                                   aCTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgc                                                                   >  1:1631589/1‑150 (MQ=255)
                                                                                    ctgctgCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTCCACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgcc                                                                  <  1:4453860/150‑1 (MQ=255)
                                                                                    ctgctgCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgcc                                                                  <  2:2897168/150‑1 (MQ=255)
                                                                                    ctgctgCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgcc                                                                  <  1:903690/150‑1 (MQ=255)
                                                                                     tgctgcCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgccg                                                                 >  1:2294881/1‑150 (MQ=255)
                                                                                     tgctgcCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgccg                                                                 <  2:2654037/150‑1 (MQ=255)
                                                                                       ctgcCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTgccgcc                                                               <  2:2951602/150‑1 (MQ=255)
                                                                                        tgcCCGGCATATCAATAGCACTGATGGCAATACGTTCTAACCGAGCGCTTGTGGTTTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCCGACGGCAGCGAAAATCACAGATAGCCGACCGAAATCGCCTGTGACAAAGGTGCCGCCa                                                              >  2:2414389/1‑150 (MQ=255)
                                                                                            cGGGCTATCACTCGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCAt                                                          <  2:2915208/150‑1 (MQ=255)
                                                                                            cGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCAt                                                          >  1:2943738/1‑150 (MQ=255)
                                                                                            cGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCAt                                                          >  2:1892842/1‑150 (MQ=255)
                                                                                             gggCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATa                                                         >  1:1004501/1‑150 (MQ=255)
                                                                                              ggCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATAc                                                        <  1:1719750/150‑1 (MQ=255)
                                                                                               gCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACa                                                       <  1:12736/150‑1 (MQ=255)
                                                                                                 tATCAATAGCACTGATGGCAAAACGCTATATCCGTGCGCTTTTGTTGTTATTCACTACACTTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCAGACAGAAATCGCCTGAGTCAAAGGTGCCGCCAGCATACAGa                                                     >  1:4344430/1‑150 (MQ=255)
                                                                                                  aTCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGAc                                                    <  1:989207/150‑1 (MQ=255)
                                                                                                   tCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGAcc                                                   <  1:3825589/150‑1 (MQ=255)
                                                                                                            cTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTcagc                                          >  1:1082989/1‑150 (MQ=255)
                                                                                                             tGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTcagca                                         >  1:910205/1‑150 (MQ=255)
                                                                                                                tgtcAATACGGGATATCCGTGCGCTTTTGTTGTTAGTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAggg                                      <  2:1655524/147‑1 (MQ=255)
                                                                                                                tGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAggg                                      <  2:3126024/150‑1 (MQ=255)
                                                                                                                  gCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGcc                                    >  2:1964004/1‑150 (MQ=255)
                                                                                                                     aTACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTa                                 <  1:1827676/150‑1 (MQ=255)
                                                                                                                     aTACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTa                                 >  2:748366/1‑150 (MQ=255)
                                                                                                                      tACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAg                                <  1:2106785/150‑1 (MQ=255)
                                                                                                                      tACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAg                                >  1:3602537/1‑150 (MQ=255)
                                                                                                                       aCGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGt                               <  2:618846/150‑1 (MQ=255)
                                                                                                                         ggTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTa                             >  2:637295/1‑150 (MQ=255)
                                                                                                                         ggTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTa                             >  2:1488832/1‑150 (MQ=255)
                                                                                                                          gTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTaa                            >  1:2535837/1‑150 (MQ=255)
                                                                                                                           tataTCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAAt                           >  2:2152703/1‑150 (MQ=255)
                                                                                                                             tatCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGt                         <  1:39069/150‑1 (MQ=255)
                                                                                                                              atCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                        <  2:626821/150‑1 (MQ=255)
                                                                                                                              atCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                        >  1:2500360/1‑150 (MQ=255)
                                                                                                                              atCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                        >  1:4362139/1‑150 (MQ=255)
                                                                                                                              atCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                        <  2:5252/150‑1 (MQ=255)
                                                                                                                                 cGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACg                     >  2:1527312/1‑150 (MQ=255)
                                                                                                                                 cGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACg                     >  2:3848397/1‑150 (MQ=255)
                                                                                                                                  gTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGc                    >  2:4466321/1‑150 (MQ=255)
                                                                                                                                      gcTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACg                <  2:2782514/150‑1 (MQ=255)
                                                                                                                                       cTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGt               <  1:1008838/150‑1 (MQ=255)
                                                                                                                                       cTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGt               >  1:459374/1‑150 (MQ=255)
                                                                                                                                        ttttGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTc              >  2:4315128/1‑150 (MQ=255)
                                                                                                                                        ttttGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTc              <  1:1173402/150‑1 (MQ=255)
                                                                                                                                         tttGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCg             >  2:4548631/1‑150 (MQ=255)
                                                                                                                                         tttGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCg             >  2:3794245/1‑150 (MQ=255)
                                                                                                                                             ttgttATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTg         <  2:1238144/150‑1 (MQ=255)
                                                                                                                                              tgttATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGa        >  2:3261649/1‑150 (MQ=255)
                                                                                                                                               gttATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGaa       >  2:819508/1‑150 (MQ=255)
                                                                                                                                                  aTTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACAt    >  2:3534351/1‑150 (MQ=255)
                                                                                                                                                  aTTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACAt    <  1:1512698/150‑1 (MQ=255)
                                                                                                                                                   ttCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATg   >  1:2469101/1‑150 (MQ=255)
                                                                                                                                                   gtCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATg   <  1:3889807/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAATACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGA  >  CP009273/300621‑300918

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: