Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 262,114 C→A 57.0% D51Y (GAT→TAT)  yafY ← lipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273262,1140CA57.0% 22.5 / inf 226D51Y (GAT→TAT) yafYlipoprotein, inner membrane; overproduction stimulates degP expression; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (44/53);  new base A (59/70);  total (103/123)
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.88e-01

Only 100 of 229 total aligned reads displayed.

TCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTC  >  CP009273/261967‑262262
                                                                                                                                                   |                                                                                                                                                    
tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatcat                                                                                                                                                    <  2:1402233/150‑1 (MQ=255)
 catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACataatc                                                                                                                                                   <  1:3616086/150‑1 (MQ=255)
   tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCt                                                                                                                                                 <  1:3581415/150‑1 (MQ=255)
     aGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCtgg                                                                                                                                               >  2:3781364/1‑150 (MQ=255)
        ggCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggtgg                                                                                                                                            >  2:4067666/1‑150 (MQ=255)
         gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGc                                                                                                                                           >  2:687135/1‑150 (MQ=255)
          cTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCt                                                                                                                                          <  1:932877/150‑1 (MQ=255)
          cTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCt                                                                                                                                          <  2:508194/150‑1 (MQ=255)
           tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTg                                                                                                                                         <  1:2324979/150‑1 (MQ=255)
           tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTg                                                                                                                                         <  1:3467648/150‑1 (MQ=255)
            ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGa                                                                                                                                        >  1:1787344/1‑150 (MQ=255)
                ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGTGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACcac                                                                                                                                    <  1:3934851/150‑1 (MQ=255)
                           gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAc                                                                                                                         >  2:1170996/1‑150 (MQ=255)
                              tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                                                                      <  1:2786177/150‑1 (MQ=255)
                              tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                                                                      <  1:4140057/150‑1 (MQ=255)
                              tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                                                                      <  1:3649466/150‑1 (MQ=255)
                                tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGt                                                                                                                    <  1:2523427/150‑1 (MQ=255)
                                  cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTgg                                                                                                                  >  1:2615297/1‑150 (MQ=255)
                                   aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                                                                                 <  1:2069384/150‑1 (MQ=255)
                                     gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                                                                               <  2:2614508/150‑1 (MQ=255)
                                      aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCTCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTg                                                                                                              <  1:4177201/150‑1 (MQ=255)
                                       ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                                                                             <  2:389098/150‑1 (MQ=255)
                                        gtagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAg                                                                                                            <  2:2576435/150‑1 (MQ=255)
                                          agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGt                                                                                                          >  1:2040112/1‑150 (MQ=255)
                                           gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTc                                                                                                         >  2:345351/1‑150 (MQ=255)
                                           gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTc                                                                                                         <  2:4029078/150‑1 (MQ=255)
                                            tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc                                                                                                        <  1:85496/150‑1 (MQ=255)
                                            tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc                                                                                                        <  2:3993329/150‑1 (MQ=255)
                                             agtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCt                                                                                                       <  1:705235/150‑1 (MQ=255)
                                              gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCtt                                                                                                      <  1:2859872/150‑1 (MQ=255)
                                              gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCtt                                                                                                      >  2:2517292/1‑150 (MQ=255)
                                               taACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTa                                                                                                     <  2:3279533/150‑1 (MQ=255)
                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                                                                    <  1:2908850/150‑1 (MQ=255)
                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                                                                    <  1:8445/150‑1 (MQ=255)
                                                 aCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaaa                                                                                                   <  1:1328702/150‑1 (MQ=255)
                                                 aCTGCCATCCTCGCAGCCATCGGGCGACAGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTACGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGTTGAGGTCCTTaaa                                                                                                   >  2:23587/1‑150 (MQ=255)
                                                   tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt                                                                                                 >  1:1500143/1‑150 (MQ=255)
                                                   tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt                                                                                                 <  1:1839172/150‑1 (MQ=255)
                                                   tGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACATGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGCTCCTTAAAtt                                                                                                 >  1:1086914/1‑150 (MQ=255)
                                                     ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTa                                                                                               <  2:3085088/150‑1 (MQ=255)
                                                       aTCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACt                                                                                             <  1:3397225/150‑1 (MQ=255)
                                                         ccTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGa                                                                                           <  2:102975/150‑1 (MQ=255)
                                                           tCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                                                                         <  1:2454314/150‑1 (MQ=255)
                                                               aGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                                                                     <  2:934775/150‑1 (MQ=255)
                                                                gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                                                                    >  2:1696372/1‑150 (MQ=255)
                                                                   aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAg                                                                                 <  2:3918577/150‑1 (MQ=255)
                                                                      gggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc                                                                              >  2:635315/1‑150 (MQ=255)
                                                                         cGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                           >  1:2360503/1‑150 (MQ=255)
                                                                         cGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                           <  2:1225055/150‑1 (MQ=255)
                                                                           aCTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCAc                                                                         >  2:2916394/1‑150 (MQ=255)
                                                                           aCTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCAc                                                                         <  2:3304224/150‑1 (MQ=255)
                                                                           aCTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCAc                                                                         <  2:1890387/150‑1 (MQ=255)
                                                                              ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTTGTGGCTGACCACATCACTGAAACGCGGTGGGGTGATGTCCTTAAATTTACTGATGGCCTTCAGGACATCACTTc                                                                      <  2:3538677/150‑1 (MQ=255)
                                                                              ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                      <  2:1334839/150‑1 (MQ=255)
                                                                              ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                      >  1:3848522/1‑150 (MQ=255)
                                                                                cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                                                                    >  1:3021421/1‑150 (MQ=255)
                                                                                   gTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                                                                 >  2:780754/1‑150 (MQ=255)
                                                                                    tACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAt                                                                >  1:1118834/1‑150 (MQ=255)
                                                                                     aCGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAACCGCGGTGGGGTGAGGTCCTTCAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc                                                               >  2:2195387/1‑150 (MQ=255)
                                                                                        tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                                                            >  1:668142/1‑150 (MQ=255)
                                                                                          aaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg                                                          <  2:3458392/150‑1 (MQ=255)
                                                                                           aaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGc                                                         <  2:1098385/150‑1 (MQ=255)
                                                                                           aaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGc                                                         >  1:2749727/1‑150 (MQ=255)
                                                                                               tGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg                                                     >  1:2564242/1‑150 (MQ=255)
                                                                                               tGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg                                                     <  2:3382916/150‑1 (MQ=255)
                                                                                                  aaGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgaga                                                  >  1:4051260/1‑150 (MQ=255)
                                                                                                     gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc                                               >  2:356/1‑150 (MQ=255)
                                                                                                     gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc                                               >  2:2199846/1‑150 (MQ=255)
                                                                                                     gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc                                               >  2:150/1‑150 (MQ=255)
                                                                                                     gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc                                               >  2:158/1‑150 (MQ=255)
                                                                                                         gACCGGATAAGTAGTCATCCTGTGACCAGTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                           <  2:244866/150‑1 (MQ=255)
                                                                                                         gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                           <  2:1764352/150‑1 (MQ=255)
                                                                                                           ccGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtg                                         <  2:627889/150‑1 (MQ=255)
                                                                                                           ccGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtg                                         <  2:2800174/150‑1 (MQ=255)
                                                                                                             ggATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                       >  1:1629275/1‑150 (MQ=255)
                                                                                                              gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                      <  2:1813947/150‑1 (MQ=255)
                                                                                                                 aaGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCa                                   <  2:2097364/150‑1 (MQ=255)
                                                                                                                    tagtCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAcc                                >  2:2429682/1‑150 (MQ=255)
                                                                                                                    tagtCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAcc                                >  1:2673821/1‑150 (MQ=255)
                                                                                                                    tagtCAACCTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACAGAAACGCGCTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCaaac                                >  1:3292818/1‑148 (MQ=255)
                                                                                                                     agtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCa                               >  2:2913236/1‑150 (MQ=255)
                                                                                                                     agtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCa                               <  2:2595533/150‑1 (MQ=255)
                                                                                                                           ccTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCca                         <  2:3609248/150‑1 (MQ=255)
                                                                                                                           ccTGTGACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCca                         >  1:325274/1‑150 (MQ=255)
                                                                                                                             tgtgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagc                       >  1:728809/1‑150 (MQ=255)
                                                                                                                              gtgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagca                      >  1:657440/1‑150 (MQ=255)
                                                                                                                               tgACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagcag                     <  1:727326/150‑1 (MQ=255)
                                                                                                                                gACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCaga                    <  1:1006908/150‑1 (MQ=255)
                                                                                                                                gACCATTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCaga                    >  1:2249659/1‑150 (MQ=255)
                                                                                                                                 aCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCagag                   <  2:3636716/150‑1 (MQ=255)
                                                                                                                                  ccATTCTTCGCTGACATAATCCTGGGGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAgg                  <  2:3979708/150‑1 (MQ=255)
                                                                                                                                    aTTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa                >  2:1968439/1‑150 (MQ=255)
                                                                                                                                    aTTCTTCGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa                <  1:1036834/150‑1 (MQ=255)
                                                                                                                                     ttcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg               <  1:2006107/150‑1 (MQ=255)
                                                                                                                                     ttcttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg               >  1:1455398/1‑150 (MQ=255)
                                                                                                                                       cttcGCTGACATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtg             <  1:3181506/150‑1 (MQ=255)
                                                                                                                                          cGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGttt          <  2:1911714/150‑1 (MQ=255)
                                                                                                                                             tGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc       >  2:2662871/1‑150 (MQ=255)
                                                                                                                                             tGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc       <  1:329791/150‑1 (MQ=255)
                                                                                                                                               aCATAATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCGGTCATCGCAGGCGATGAGAAACAGTGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACgg     >  2:313404/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTC  >  CP009273/261967‑262262

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: