Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,643,336 C→A 7.0% D771Y (GAC→TAC)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,643,3360CA 7.0% 1025.0 / 33.6 269D771Y (GAC→TAC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (117/133);  new base A (9/10);  total (126/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 271 total aligned reads displayed.

ACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGAT  >  CP009273/2643187‑2643485
                                                                                                                                                     |                                                                                                                                                     
 cGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCg                                                                                                                                                      >  1:3821356/1‑150 (MQ=255)
   aTTCCAGGTTTTATCGACCGGATTCGTCAGATCCAGCCCTTGAGCCCGAACCTCAATGTCTTGCCACCTCAGCGGCCATTCCCTGGTCTAGACCATCGCATTCCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGcc                                                                                                                                                    <  2:3046907/150‑1 (MQ=255)
     tCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAg                                                                                                                                                  <  2:1558871/150‑1 (MQ=255)
      ccAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTACTGGTGTCGCCAgg                                                                                                                                                 <  2:1599784/150‑1 (MQ=255)
      ccAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAgg                                                                                                                                                 <  2:4176856/150‑1 (MQ=255)
       cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAggg                                                                                                                                                >  1:490232/1‑150 (MQ=255)
        aGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGc                                                                                                                                               >  2:2199479/1‑150 (MQ=255)
         ggTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCg                                                                                                                                              >  2:3157839/1‑150 (MQ=255)
           ttttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATTGTGTCGCCAGGGCGAt                                                                                                                                            <  1:3010477/150‑1 (MQ=255)
           ttttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt                                                                                                                                            <  1:3520962/150‑1 (MQ=255)
           ttttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt                                                                                                                                            >  2:1895625/1‑150 (MQ=255)
            tttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATa                                                                                                                                           <  2:4421155/150‑1 (MQ=255)
             ttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATaa                                                                                                                                          <  1:4057753/150‑1 (MQ=255)
                 cGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTgg                                                                                                                                      >  2:3630968/1‑150 (MQ=255)
                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGc                                                                                                                                     <  2:590696/150‑1 (MQ=255)
                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGc                                                                                                                                     >  2:3344208/1‑150 (MQ=255)
                   aCCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCt                                                                                                                                    >  1:4304651/1‑150 (MQ=255)
                   aCCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCt                                                                                                                                    <  2:3531870/150‑1 (MQ=255)
                    ccGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCtt                                                                                                                                   <  2:3719021/150‑1 (MQ=255)
                    ccGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCtt                                                                                                                                   <  1:559426/150‑1 (MQ=255)
                     cGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGGTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCttt                                                                                                                                  >  1:806186/1‑150 (MQ=255)
                       gAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTAt                                                                                                                                >  1:225315/1‑150 (MQ=255)
                       gAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTAt                                                                                                                                <  2:3296254/150‑1 (MQ=255)
                         aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATcc                                                                                                                              <  2:2492004/150‑1 (MQ=255)
                          tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCa                                                                                                                             >  1:2204386/1‑150 (MQ=255)
                          tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCa                                                                                                                             >  1:374235/1‑150 (MQ=255)
                          tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCa                                                                                                                             <  1:4572211/150‑1 (MQ=255)
                           cGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGAATAACCTTTACCCTACGTTGTCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAg                                                                                                                            >  2:3092459/1‑150 (MQ=255)
                            gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGt                                                                                                                           >  2:669393/1‑150 (MQ=255)
                             tCAGATCCAGCCCTTGAGCCCGAACATCACTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGtt                                                                                                                          <  2:2119738/150‑1 (MQ=255)
                                gATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCa                                                                                                                       <  2:850882/150‑1 (MQ=255)
                                   ccAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCg                                                                                                                    <  1:323869/150‑1 (MQ=255)
                                     aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTg                                                                                                                  >  1:3406373/1‑150 (MQ=255)
                                     aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTg                                                                                                                  <  2:1046309/150‑1 (MQ=255)
                                        ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGaca                                                                                                               >  2:451575/1‑150 (MQ=255)
                                         cTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGacac                                                                                                              >  2:403618/1‑150 (MQ=255)
                                         cTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGacac                                                                                                              <  1:503264/150‑1 (MQ=255)
                                          ttGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACg                                                                                                             >  1:1026967/1‑150 (MQ=255)
                                            gAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggt                                                                                                           >  2:4254535/1‑150 (MQ=255)
                                             aGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtc                                                                                                          <  1:1261636/150‑1 (MQ=255)
                                              gCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcg                                                                                                         <  1:3954840/150‑1 (MQ=255)
                                                ccGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcggg                                                                                                       >  1:2629771/1‑150 (MQ=255)
                                                 cGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGACGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggt                                                                                                      >  1:296656/1‑150 (MQ=255)
                                                  gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                                                                     >  1:336236/1‑150 (MQ=255)
                                                  gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                                                                     <  1:2298022/150‑1 (MQ=255)
                                                   aaCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtcg                                                                                                    <  1:1898923/150‑1 (MQ=255)
                                                                tGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTa                                                                                       >  2:3682786/1‑150 (MQ=255)
                                                                 gCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAc                                                                                      >  2:2418597/1‑150 (MQ=255)
                                                                 gCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAAACTTTACCCGCCGTAGGCGCGGCGATAAGCAACTTAATGGTGTCGCCAGCGCTATAACACGCTTTATCCAGATACACCGTGACACGGTCGGGTCGCACTGCCCCGCTAc                                                                                      >  1:1642909/1‑150 (MQ=255)
                                                                  ccaccaCAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAcc                                                                                     >  1:531039/1‑150 (MQ=255)
                                                                     ccacAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTc                                                                                  <  2:180937/150‑1 (MQ=255)
                                                                         aGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                              >  2:3873119/1‑150 (MQ=255)
                                                                         aGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                              <  1:1171848/150‑1 (MQ=255)
                                                                         aGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                              >  1:2987218/1‑150 (MQ=255)
                                                                         aGCGGCCCTTCACTGGACTCGACCATCGCATAACCTATACCCGCCGTTGGCGCGGCGATATGCAACTTAAGGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                              >  2:1516816/1‑150 (MQ=255)
                                                                          gCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgt                                                                             <  1:1042072/150‑1 (MQ=255)
                                                                               ccTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTcc                                                                        >  1:2513306/1‑150 (MQ=255)
                                                                                cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCt                                                                       <  1:1355412/150‑1 (MQ=255)
                                                                                cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCt                                                                       <  1:1355403/150‑1 (MQ=255)
                                                                                  tCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGc                                                                     >  1:1140717/1‑150 (MQ=255)
                                                                                       ggACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCt                                                                <  2:252273/150‑1 (MQ=255)
                                                                                       ggACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCt                                                                <  1:3695445/150‑1 (MQ=255)
                                                                                        gACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCta                                                               <  1:2011487/150‑1 (MQ=255)
                                                                                        gACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCta                                                               <  1:2218428/150‑1 (MQ=255)
                                                                                                cATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagccagc                                                       >  1:181884/1‑150 (MQ=255)
                                                                                                  tCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGccc                                                     <  2:1167358/150‑1 (MQ=255)
                                                                                                  tCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGccc                                                     <  2:160457/150‑1 (MQ=255)
                                                                                                   cGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCa                                                    <  1:4474111/150‑1 (MQ=255)
                                                                                                      aTAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGaa                                                 >  1:1770575/1‑150 (MQ=255)
                                                                                                      aTAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGaa                                                 <  2:3110042/150‑1 (MQ=255)
                                                                                                        aaCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAAc                                               <  2:4060596/150‑1 (MQ=255)
                                                                                                        aaCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAAc                                               <  1:3971508/150‑1 (MQ=255)
                                                                                                         cccTTTACCCGCCGTTGGCGCGGCGAAATGCAACATAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                              <  2:2895036/149‑1 (MQ=255)
                                                                                                         aCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                              <  2:2894801/150‑1 (MQ=255)
                                                                                                         aCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                              >  1:4025763/1‑150 (MQ=255)
                                                                                                          ccTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCAAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGTTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGa                                             <  1:1766710/150‑1 (MQ=255)
                                                                                                           cTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGaa                                            <  1:721721/150‑1 (MQ=255)
                                                                                                            tttACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAac                                           >  1:2372431/1‑150 (MQ=255)
                                                                                                               aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACact                                        <  2:634719/150‑1 (MQ=255)
                                                                                                               aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACact                                        >  2:1950965/1‑150 (MQ=255)
                                                                                                               aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACact                                        <  1:4134506/150‑1 (MQ=255)
                                                                                                                cccGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACacta                                       >  2:2321609/1‑150 (MQ=255)
                                                                                                                 ccgccgTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                      <  2:1271689/150‑1 (MQ=255)
                                                                                                                  cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactact                                     <  2:1293413/150‑1 (MQ=255)
                                                                                                                    ccgTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGa                                   >  2:3255637/1‑150 (MQ=255)
                                                                                                                    ccgTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGa                                   >  2:1123502/1‑150 (MQ=255)
                                                                                                                    ccgTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGa                                   <  1:4060678/150‑1 (MQ=255)
                                                                                                                             cggcgATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCCGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCAt                          >  1:1183264/1‑150 (MQ=255)
                                                                                                                             cggcgATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCAt                          >  1:2682072/1‑150 (MQ=255)
                                                                                                                               gcgATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTc                        <  1:308263/150‑1 (MQ=255)
                                                                                                                                 gATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCgg                      >  2:9645/1‑150 (MQ=255)
                                                                                                                                 gATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCgg                      >  2:157786/1‑150 (MQ=255)
                                                                                                                                       cAACTTACTGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCttt                <  2:2774982/150‑1 (MQ=255)
                                                                                                                                         aCTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGa              >  1:1418191/1‑150 (MQ=255)
                                                                                                                                           ttAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACt            <  1:2479754/150‑1 (MQ=255)
                                                                                                                                           ttAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACt            <  1:3633261/150‑1 (MQ=255)
                                                                                                                                              aTGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTcc         >  2:1451067/1‑150 (MQ=255)
                                                                                                                                                   gtCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATATCCAGCCCAGAAAACAACACAACTGACCGCTTAATTCGGCGCTTTTACTTCCAGACg    >  1:936472/1‑150 (MQ=255)
                                                                                                                                                   gtCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACg    <  1:3091899/150‑1 (MQ=255)
                                                                                                                                                   gtCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACg    <  1:3104755/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGAT  >  CP009273/2643187‑2643485

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: