Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,421,791 G→T 32.0% R4S (CGC→AGC)  ytfE ← iron‑sulfur cluster repair protein RIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,421,7910GT32.0% 548.4 / 271.5 300R4S (CGC→AGC) ytfEiron‑sulfur cluster repair protein RIC
Reads supporting (aligned to +/- strand):  ref base G (78/126);  new base T (38/58);  total (116/184)
Fisher's exact test for biased strand distribution p-value = 8.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

AGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGTTT  >  CP009273/4421642‑4421940
                                                                                                                                                     |                                                                                                                                                     
  cTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGAt                                                                                                                                                     <  2:2476839/150‑1 (MQ=255)
   ttCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATa                                                                                                                                                    <  1:2975384/150‑1 (MQ=255)
     cAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAg                                                                                                                                                  <  2:2898240/150‑1 (MQ=255)
     cAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAg                                                                                                                                                  <  2:3026957/150‑1 (MQ=255)
      aaTGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGc                                                                                                                                                 >  1:3136270/1‑150 (MQ=255)
       aTGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGcc                                                                                                                                                >  1:958071/1‑150 (MQ=255)
         gACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCAt                                                                                                                                              <  1:4722841/150‑1 (MQ=255)
          aCCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATa                                                                                                                                             <  1:1478416/150‑1 (MQ=255)
           ccTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAg                                                                                                                                            <  1:2274281/150‑1 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCt                                                                                                                                          >  1:2940801/1‑150 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCt                                                                                                                                          <  1:61148/150‑1 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCt                                                                                                                                          <  2:424732/150‑1 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCt                                                                                                                                          <  1:1012589/150‑1 (MQ=255)
               aaCATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGa                                                                                                                                        <  1:3707629/150‑1 (MQ=255)
                aCATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGAt                                                                                                                                       <  2:4753156/150‑1 (MQ=255)
                aCATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGAt                                                                                                                                       <  1:1674631/150‑1 (MQ=255)
                 taTCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATa                                                                                                                                      >  2:4140664/2‑150 (MQ=255)
                 cATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATa                                                                                                                                      <  1:2621951/150‑1 (MQ=255)
                  aTCCAGTACTTTACGAGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATAc                                                                                                                                     >  2:1372490/1‑150 (MQ=255)
                    ccAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTGACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCt                                                                                                                                   <  2:4713574/150‑1 (MQ=255)
                      aGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCa                                                                                                                                 >  2:1839368/1‑150 (MQ=255)
                       gTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCAt                                                                                                                                <  1:4008504/150‑1 (MQ=255)
                       gTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                                                                                                                >  2:812983/1‑150 (MQ=255)
                       gTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                                                                                                                >  2:2728177/1‑150 (MQ=255)
                       gTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                                                                                                                >  2:70463/1‑150 (MQ=255)
                        ttCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAtt                                                                                                                               <  1:1185025/150‑1 (MQ=255)
                          cTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATATAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCt                                                                                                                             <  1:3731359/150‑1 (MQ=255)
                          cTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCt                                                                                                                             >  2:2620830/1‑150 (MQ=255)
                           tttACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCtt                                                                                                                            <  2:3375122/150‑1 (MQ=255)
                            ttACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTta                                                                                                                           <  2:375084/150‑1 (MQ=255)
                             tACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaa                                                                                                                          <  2:3937460/150‑1 (MQ=255)
                              aCGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaat                                                                                                                         <  2:3408721/150‑1 (MQ=255)
                               cGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaata                                                                                                                        <  1:2442272/150‑1 (MQ=255)
                                gTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTtaataa                                                                                                                       <  2:4003714/150‑1 (MQ=255)
                                gTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaataa                                                                                                                       <  1:180204/150‑1 (MQ=255)
                                 tGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGACTAGAGAGCGCCAGTTCACCTAAAGGATGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAataag                                                                                                                      >  2:1773844/1‑150 (MQ=255)
                                 tGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAataag                                                                                                                      >  1:4187026/1‑150 (MQ=255)
                                 tGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataag                                                                                                                      <  2:4500298/150‑1 (MQ=255)
                                   ccgccgCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAataagat                                                                                                                    <  2:1949169/150‑1 (MQ=255)
                                   ccgccgCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataatat                                                                                                                    <  1:2311529/150‑1 (MQ=255)
                                   ccgccgCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagat                                                                                                                    <  2:62027/150‑1 (MQ=255)
                                    cgccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagata                                                                                                                   <  2:1216415/150‑1 (MQ=255)
                                     gccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAataagataa                                                                                                                  <  2:3856943/150‑1 (MQ=255)
                                     gccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagataa                                                                                                                  <  1:1911779/150‑1 (MQ=255)
                                      ccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagataag                                                                                                                 <  2:2009078/150‑1 (MQ=255)
                                           cgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTa                                                                                                            <  1:3784642/150‑1 (MQ=255)
                                             cgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                                                                                                          <  2:3697297/150‑1 (MQ=255)
                                             cgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                                                                                                          <  2:556852/150‑1 (MQ=255)
                                                cAGCGTCCGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTa                                                                                                       <  1:2567106/150‑1 (MQ=255)
                                                 aGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaa                                                                                                      >  2:1296666/1‑150 (MQ=255)
                                                 aGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaa                                                                                                      <  1:4536610/150‑1 (MQ=255)
                                                  gCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa                                                                                                     <  1:2919560/150‑1 (MQ=255)
                                                     tCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGc                                                                                                  >  2:1835211/1‑150 (MQ=255)
                                                      cTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCa                                                                                                 >  2:3597117/1‑150 (MQ=255)
                                                          ttACCGCCACAGCAGTAATCCATATCATATTAACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAAAAAGATAAGATGTATATTAAATGCATCta                                                                                             >  1:4297487/1‑149 (MQ=255)
                                                            aCCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTa                                                                                           >  2:882843/1‑150 (MQ=255)
                                                            aCCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTa                                                                                           <  1:714668/150‑1 (MQ=255)
                                                               gCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAgg                                                                                        <  1:1250490/150‑1 (MQ=255)
                                                                 cacaGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCa                                                                                      <  2:4071270/150‑1 (MQ=255)
                                                                        gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCta                                                                               <  1:3976855/150‑1 (MQ=255)
                                                                           aTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATaa                                                                            <  1:928065/150‑1 (MQ=255)
                                                                             ccATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAAcc                                                                          >  1:4755062/1‑150 (MQ=255)
                                                                                  tcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTAc                                                                     >  1:3647360/1‑150 (MQ=255)
                                                                                     tatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTa                                                                  <  2:2241074/150‑1 (MQ=255)
                                                                                      atTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAg                                                                 >  1:3514208/1‑150 (MQ=255)
                                                                                      atTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAg                                                                 >  1:2229041/1‑150 (MQ=255)
                                                                                        ttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCa                                                               >  1:2686459/1‑150 (MQ=255)
                                                                                        ttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCa                                                               <  2:2294617/150‑1 (MQ=255)
                                                                                            gAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAgg                                                           >  2:4376494/1‑150 (MQ=255)
                                                                                            gAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAgg                                                           >  2:235619/1‑150 (MQ=255)
                                                                                             aaaCAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGc                                                          <  1:1146666/150‑1 (MQ=255)
                                                                                                cAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTa                                                       <  1:3074923/150‑1 (MQ=255)
                                                                                                cAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTa                                                       >  2:370663/1‑150 (MQ=255)
                                                                                                  gagCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTAtt                                                     <  1:3444478/150‑1 (MQ=255)
                                                                                                       gAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACt                                                >  1:3555952/1‑150 (MQ=255)
                                                                                                       gAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACt                                                <  1:3271092/150‑1 (MQ=255)
                                                                                                          gcgcGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCAc                                             <  1:1645955/150‑1 (MQ=255)
                                                                                                            gcgAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTg                                           <  1:3850908/150‑1 (MQ=255)
                                                                                                             cgAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGc                                          <  2:1571054/150‑1 (MQ=255)
                                                                                                              gAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGcc                                         <  2:2571724/150‑1 (MQ=255)
                                                                                                               aGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCa                                        <  1:1372490/150‑1 (MQ=255)
                                                                                                                  aaTAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCg                                     >  1:1500768/1‑150 (MQ=255)
                                                                                                                      gagagCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTc                                 <  1:1382051/150‑1 (MQ=255)
                                                                                                                       agagCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCa                                <  1:181358/150‑1 (MQ=255)
                                                                                                                        gagCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCaa                               >  2:3091255/1‑150 (MQ=255)
                                                                                                                                gTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTga                       <  1:245080/150‑1 (MQ=255)
                                                                                                                                gTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTga                       <  1:2637988/150‑1 (MQ=255)
                                                                                                                                 ttCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTGTAAGGCAAAAGCTATAACCCTTACGTAGCAGCGGCTTATTTAACTCACTGCCCGCGTGTCAACTTTCTgag                      >  2:3658123/1‑150 (MQ=255)
                                                                                                                                    aCCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACa                   <  1:2011238/150‑1 (MQ=255)
                                                                                                                                    aCCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACa                   <  1:3199893/150‑1 (MQ=255)
                                                                                                                                    aCCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACa                   <  1:106901/150‑1 (MQ=255)
                                                                                                                                      cTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGa                 >  2:460942/1‑150 (MQ=255)
                                                                                                                                       tAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGAt                >  1:615937/1‑150 (MQ=255)
                                                                                                                                        aaaGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGAtt               >  2:2781500/1‑150 (MQ=255)
                                                                                                                                         aaGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTc              >  2:4328333/1‑150 (MQ=255)
                                                                                                                                          aGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCa             >  2:1129386/1‑150 (MQ=255)
                                                                                                                                          aGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTGTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCa             <  1:1283802/150‑1 (MQ=255)
                                                                                                                                           ggTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAg            >  1:987283/1‑150 (MQ=255)
                                                                                                                                                  cgcgATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAgg     >  1:1697075/1‑150 (MQ=255)
                                                                                                                                                   gcgATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTACGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGt    >  2:4463020/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGTTT  >  CP009273/4421642‑4421940

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: