Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,247 A→T 93.4% V57E (GTG→GAG)  yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,2470AT93.4% 830.2 / 37.7 242V57E (GTG→GAG) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (5/11);  new base T (120/106);  total (125/117)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 244 total aligned reads displayed.

TTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTA  >  CP009273/3994098‑3994393
                                                                                                                                                     |                                                                                                                                                  
    gATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCg                                                                                                                                                <  1:802306/150‑1 (MQ=255)
          cAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGt                                                                                                                                          >  1:1612461/1‑150 (MQ=255)
           aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                         <  1:1858405/150‑1 (MQ=255)
            aaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGttt                                                                                                                                        >  2:1193881/1‑150 (MQ=255)
             aCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTg                                                                                                                                       >  1:2691396/1‑150 (MQ=255)
                  gcggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAt                                                                                                                                  <  1:1133761/150‑1 (MQ=255)
                  gcggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAt                                                                                                                                  <  1:231087/150‑1 (MQ=255)
                         aTAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGt                                                                                                                           <  2:3111672/150‑1 (MQ=255)
                         aTAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGGTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGt                                                                                                                           <  2:1372228/150‑1 (MQ=255)
                          tAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGta                                                                                                                          >  2:2870520/1‑150 (MQ=255)
                           aGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtat                                                                                                                         <  1:16432/150‑1 (MQ=255)
                           aGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtat                                                                                                                         >  2:395028/1‑150 (MQ=255)
                           aGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtat                                                                                                                         <  2:3524276/150‑1 (MQ=255)
                            gCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtata                                                                                                                        >  2:2663474/1‑150 (MQ=255)
                                  atcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTg                                                                                                                  <  2:1039453/150‑1 (MQ=255)
                                    catAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                                <  2:1454322/150‑1 (MQ=255)
                                      tAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACg                                                                                                              <  1:1176930/150‑1 (MQ=255)
                                      tAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACg                                                                                                              <  1:1126571/150‑1 (MQ=255)
                                          tGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCa                                                                                                          <  2:3399440/150‑1 (MQ=255)
                                               ccgTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTc                                                                                                     <  1:1012795/150‑1 (MQ=255)
                                                 gTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCgg                                                                                                   >  2:53580/1‑150 (MQ=255)
                                                   tCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGAt                                                                                                 >  1:4048150/1‑150 (MQ=255)
                                                    cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATc                                                                                                <  2:2226630/150‑1 (MQ=255)
                                                       cGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTga                                                                                             >  1:616945/1‑150 (MQ=255)
                                                         ccATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTgaga                                                                                           <  1:265808/150‑1 (MQ=255)
                                                           aTCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGt                                                                                         >  1:3710571/1‑150 (MQ=255)
                                                              tGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcg                                                                                      >  2:626399/1‑150 (MQ=255)
                                                                   ctgaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatc                                                                                 >  2:116756/1‑150 (MQ=255)
                                                                   ctgaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatc                                                                                 <  1:1320516/150‑1 (MQ=255)
                                                                    tgaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatca                                                                                >  2:2232672/1‑150 (MQ=255)
                                                                     gaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcat                                                                               <  1:3533616/150‑1 (MQ=255)
                                                                      accctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcatc                                                                              >  1:1653092/1‑150 (MQ=255)
                                                                         ctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGc                                                                           >  2:3424557/1‑150 (MQ=255)
                                                                          tgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCg                                                                          <  2:4092528/150‑1 (MQ=255)
                                                                          tgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCGTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCg                                                                          <  1:2227792/150‑1 (MQ=255)
                                                                           gacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCgg                                                                         >  1:1561006/1‑150 (MQ=255)
                                                                             cTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCa                                                                       <  1:1504007/150‑1 (MQ=255)
                                                                                    tatCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCaa                                                                >  1:26326/1‑150 (MQ=255)
                                                                                         gCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAg                                                           >  2:1662106/1‑150 (MQ=255)
                                                                                          ccAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAgg                                                          >  1:3860816/1‑150 (MQ=255)
                                                                                            agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        >  1:3247464/1‑150 (MQ=255)
                                                                                            agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        >  2:3866746/1‑150 (MQ=255)
                                                                                              cagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAg                                                      >  2:1101939/1‑150 (MQ=255)
                                                                                              cagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAg                                                      >  2:3200643/1‑150 (MQ=255)
                                                                                                 agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTaa                                                   <  2:4209244/150‑1 (MQ=255)
                                                                                                 agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTaa                                                   <  2:4135677/150‑1 (MQ=255)
                                                                                                 agCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTaa                                                   >  1:2524151/1‑150 (MQ=255)
                                                                                                    aTGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGt                                                >  1:1026607/1‑150 (MQ=255)
                                                                                                    aTGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGt                                                >  1:346149/1‑150 (MQ=255)
                                                                                                     tGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTg                                               >  2:4306868/1‑150 (MQ=255)
                                                                                                     tGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTg                                               <  1:3209519/150‑1 (MQ=255)
                                                                                                          cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGc                                          <  2:4165122/150‑1 (MQ=255)
                                                                                                          cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGc                                          <  1:821985/150‑1 (MQ=255)
                                                                                                          cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGc                                          <  2:3031773/150‑1 (MQ=255)
                                                                                                           gTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                         <  2:3831368/150‑1 (MQ=255)
                                                                                                           gTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                         >  1:1064804/1‑150 (MQ=255)
                                                                                                            ttCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCga                                        >  2:472706/1‑150 (MQ=255)
                                                                                                            ttCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCga                                        >  1:1197976/1‑150 (MQ=255)
                                                                                                                ccACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaa                                    <  2:3120875/150‑1 (MQ=255)
                                                                                                                  aCGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaaaa                                  >  1:2982190/1‑150 (MQ=255)
                                                                                                                  aCGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaaaa                                  >  1:1756537/1‑150 (MQ=255)
                                                                                                                   cGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAAt                                 >  2:2448989/1‑150 (MQ=255)
                                                                                                                    gCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATc                                <  1:3894027/150‑1 (MQ=255)
                                                                                                                     cTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATcc                               >  1:274668/1‑150 (MQ=255)
                                                                                                                      ttCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                              >  1:4241349/1‑150 (MQ=255)
                                                                                                                      ttCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                              <  1:217290/150‑1 (MQ=255)
                                                                                                                      ttCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                              <  1:2507687/150‑1 (MQ=255)
                                                                                                                        cGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGtt                            <  2:1564284/150‑1 (MQ=255)
                                                                                                                        cGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGtt                            <  1:1935087/150‑1 (MQ=255)
                                                                                                                         gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                           <  1:378947/150‑1 (MQ=255)
                                                                                                                         gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                           <  1:212181/150‑1 (MQ=255)
                                                                                                                           cATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTga                         >  1:1249201/1‑150 (MQ=255)
                                                                                                                             ttGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTgaga                       <  1:3055418/150‑1 (MQ=255)
                                                                                                                               ggCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAtt                     >  2:1870099/1‑150 (MQ=255)
                                                                                                                               ggCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAtt                     <  2:1582184/150‑1 (MQ=255)
                                                                                                                                gCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTc                    >  1:632120/1‑150 (MQ=255)
                                                                                                                                gCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTc                    <  1:4041268/150‑1 (MQ=255)
                                                                                                                                 cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                   <  1:1598779/150‑1 (MQ=255)
                                                                                                                                 cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                   >  1:64481/1‑150 (MQ=255)
                                                                                                                                  tttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCaa                  >  1:2380962/1‑150 (MQ=255)
                                                                                                                                    tttGCCAGTACATTTTCTCCCGCTCTGTTTTAGGATTACGGGTATACGCCTGTAGGTGCCGGGTCGTATCGGATCGTGCGCAATCTTCAGCGGCAACAGAAAACACGGTATGTATGGTGTGGGCGAGAAAAATCCGTTTGAGATTCAAgg                <  2:279313/150‑1 (MQ=255)
                                                                                                                                     ttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGt               <  2:2374032/150‑1 (MQ=255)
                                                                                                                                      tGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg              <  2:768529/150‑1 (MQ=255)
                                                                                                                                      tGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg              <  2:3167693/150‑1 (MQ=255)
                                                                                                                                      tGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg              <  1:1795330/150‑1 (MQ=255)
                                                                                                                                       gCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGa             <  1:4111170/150‑1 (MQ=255)
                                                                                                                                          aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa          <  1:190052/150‑1 (MQ=255)
                                                                                                                                          aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa          >  2:1975930/1‑150 (MQ=255)
                                                                                                                                            tACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAtt        >  1:1268578/1‑150 (MQ=255)
                                                                                                                                            tACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAtt        <  2:800519/150‑1 (MQ=255)
                                                                                                                                            tACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAtt        <  2:2974411/150‑1 (MQ=255)
                                                                                                                                            tACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAtt        >  1:3826621/1‑150 (MQ=255)
                                                                                                                                             aCATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTc       <  2:2152566/150‑1 (MQ=255)
                                                                                                                                              cATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc      >  2:3918356/1‑150 (MQ=255)
                                                                                                                                              cATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc      <  2:1796318/150‑1 (MQ=255)
                                                                                                                                                ttttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGTTTCAAGGTGATCATTCCtt    >  2:666319/1‑150 (MQ=255)
                                                                                                                                                 tttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCttt   <  2:1201480/150‑1 (MQ=255)
                                                                                                                                                 tttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCttt   >  1:483001/1‑150 (MQ=255)
                                                                                                                                                  ttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTa  >  1:3022345/1‑150 (MQ=255)
                                                                                                                                                  ttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTa  <  1:3948998/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTA  >  CP009273/3994098‑3994393

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: ‑