Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,334 C→A 20.9% I141I (ATC→ATA yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,3340CA20.9% 775.7 / 146.9 285I141I (ATC→ATAyagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (109/116);  new base A (27/33);  total (136/149)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

TCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTGCC  >  CP009273/279186‑279483
                                                                                                                                                    |                                                                                                                                                     
tCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCg                                                                                                                                                      <  2:4602846/150‑1 (MQ=255)
   gCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAgtgt                                                                                                                                                   <  2:562464/150‑1 (MQ=255)
       gcagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTcg                                                                                                                                               <  2:3853810/150‑1 (MQ=255)
            tgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCgccgcc                                                                                                                                          <  1:3941397/150‑1 (MQ=255)
            tgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCgccgcc                                                                                                                                          <  1:3178972/150‑1 (MQ=255)
             gcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCt                                                                                                                                         <  2:1100963/150‑1 (MQ=255)
              cagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTg                                                                                                                                        <  1:5054537/150‑1 (MQ=255)
               agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGa                                                                                                                                       <  2:2161405/150‑1 (MQ=255)
               agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGa                                                                                                                                       <  2:5116534/150‑1 (MQ=255)
               agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGa                                                                                                                                       <  1:2708506/150‑1 (MQ=255)
                 gcggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATc                                                                                                                                     <  2:1054866/150‑1 (MQ=255)
                  cggcggcgAAGGAGAGCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATcc                                                                                                                                    <  2:3228725/150‑1 (MQ=255)
                     cggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCt                                                                                                                                 >  2:4371749/1‑150 (MQ=255)
                       gcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTcc                                                                                                                               >  1:2338297/1‑150 (MQ=255)
                        cgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCATTCGTCAGCGATCCGTGCGACGGGAGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTccc                                                                                                                              >  1:4263851/1‑150 (MQ=255)
                         gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                             >  2:2053254/1‑150 (MQ=255)
                           aGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGTTTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGc                                                                                                                           >  2:3640030/1‑150 (MQ=255)
                            ggAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGcc                                                                                                                          >  2:4037723/1‑150 (MQ=255)
                             gagaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCg                                                                                                                         <  1:1582036/150‑1 (MQ=255)
                              agaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGa                                                                                                                        >  2:1127492/1‑150 (MQ=255)
                              agaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGa                                                                                                                        >  2:3314752/1‑150 (MQ=255)
                               gaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAc                                                                                                                       <  2:58237/150‑1 (MQ=255)
                                aTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAcg                                                                                                                      <  1:2227806/150‑1 (MQ=255)
                                   aCCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcgg                                                                                                                   >  2:1830658/1‑150 (MQ=255)
                                     ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggcg                                                                                                                 >  2:4915670/1‑150 (MQ=255)
                                        cAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGCTCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCggcgggc                                                                                                              <  1:2636101/150‑1 (MQ=255)
                                        cAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTACGCCGCCTGATCCGCTCCCTGCCGACGCggcgggc                                                                                                              <  2:1241765/150‑1 (MQ=255)
                                         aaTGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGTACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCggcgggcg                                                                                                             >  1:956098/1‑150 (MQ=255)
                                         aaTGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCggcgggcg                                                                                                             >  1:1166272/1‑150 (MQ=255)
                                           tGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGt                                                                                                           <  2:4760883/150‑1 (MQ=255)
                                              cGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGAt                                                                                                        <  1:3861131/150‑1 (MQ=255)
                                               gggATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATc                                                                                                       <  2:1597804/150‑1 (MQ=255)
                                                   tCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCg                                                                                                   >  1:1828190/1‑150 (MQ=255)
                                                     ccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTa                                                                                                 <  1:3327039/150‑1 (MQ=255)
                                                      cGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAg                                                                                                >  2:2501596/1‑150 (MQ=255)
                                                      cGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAg                                                                                                <  2:4269564/150‑1 (MQ=255)
                                                       gTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGc                                                                                               <  1:2092034/150‑1 (MQ=255)
                                                             ggCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTg                                                                                         <  2:1405730/150‑1 (MQ=255)
                                                             ggCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTg                                                                                         >  1:2458518/1‑150 (MQ=255)
                                                               ccTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCg                                                                                       <  2:3730027/150‑1 (MQ=255)
                                                                   cGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaa                                                                                   <  1:698867/150‑1 (MQ=255)
                                                                    gTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaaa                                                                                  <  1:1437310/150‑1 (MQ=255)
                                                                     tCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAg                                                                                 <  1:5138709/150‑1 (MQ=255)
                                                                         cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCt                                                                             <  2:2296404/150‑1 (MQ=255)
                                                                          gATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTg                                                                            >  2:4237575/1‑150 (MQ=255)
                                                                          gATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTg                                                                            >  1:3981766/1‑150 (MQ=255)
                                                                           aTCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGc                                                                           <  2:2095403/150‑1 (MQ=255)
                                                                             ccGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGccc                                                                         <  2:1155486/150‑1 (MQ=255)
                                                                              cGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCg                                                                        <  2:2759812/150‑1 (MQ=255)
                                                                              cGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCg                                                                        <  2:632767/150‑1 (MQ=255)
                                                                              cGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCg                                                                        <  2:3737089/150‑1 (MQ=255)
                                                                                   gACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccacc                                                                   >  2:339932/1‑150 (MQ=255)
                                                                                      gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                                <  1:4952757/150‑1 (MQ=255)
                                                                                          gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTg                                                            <  2:3961806/150‑1 (MQ=255)
                                                                                           cTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgc                                                           >  1:3020700/1‑150 (MQ=255)
                                                                                             cGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcgc                                                         >  1:3225381/1‑150 (MQ=255)
                                                                                               cAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTg                                                       <  1:899321/150‑1 (MQ=255)
                                                                                               cAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTg                                                       <  1:4452066/150‑1 (MQ=255)
                                                                                                  ggCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGcc                                                    <  1:1360933/150‑1 (MQ=255)
                                                                                                   gcacgcacGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcg                                                   <  2:1935997/150‑1 (MQ=255)
                                                                                                    cacgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                                  >  1:2350504/1‑150 (MQ=255)
                                                                                                     acgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgcg                                                 <  1:2146917/150‑1 (MQ=255)
                                                                                                       gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                               <  1:734597/150‑1 (MQ=255)
                                                                                                       gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                               <  1:1188950/150‑1 (MQ=255)
                                                                                                       gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                               >  1:419115/1‑150 (MQ=255)
                                                                                                        cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                              >  2:904100/1‑150 (MQ=255)
                                                                                                        cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                              >  1:1545553/1‑150 (MQ=255)
                                                                                                        cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                              >  2:4015357/1‑150 (MQ=255)
                                                                                                         acgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGc                                             <  2:3054789/150‑1 (MQ=255)
                                                                                                         acgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGc                                             <  2:4570872/150‑1 (MQ=255)
                                                                                                         acgGTATGTTCGATTCCCTGCCGGACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGc                                             <  2:2593612/150‑1 (MQ=255)
                                                                                                           ggTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCAc                                           >  2:2842866/1‑150 (MQ=255)
                                                                                                           ggTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCAc                                           <  2:1891919/150‑1 (MQ=255)
                                                                                                               tGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGAcc                                       <  1:4352610/150‑1 (MQ=255)
                                                                                                                  tCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGc                                    >  1:1964508/1‑150 (MQ=255)
                                                                                                                   cGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGcc                                   >  1:4823205/1‑150 (MQ=255)
                                                                                                                   cGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGcc                                   <  2:4389641/150‑1 (MQ=255)
                                                                                                                     aTTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCTCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGa                                 >  2:4785638/1‑150 (MQ=255)
                                                                                                                      ttCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGAc                                >  2:2658720/1‑150 (MQ=255)
                                                                                                                       tCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACt                               <  2:1695154/150‑1 (MQ=255)
                                                                                                                        cccTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTa                              >  2:5013432/1‑150 (MQ=255)
                                                                                                                          cTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAtt                            >  1:200494/1‑150 (MQ=255)
                                                                                                                           tGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTc                           <  1:2257502/150‑1 (MQ=255)
                                                                                                                            gCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCt                          <  2:3358732/150‑1 (MQ=255)
                                                                                                                             ccGTACCGCAACGACGCGGCGATCGTGTTACGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTg                         >  2:1955138/1‑150 (MQ=255)
                                                                                                                              cGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTgg                        >  1:1072793/1‑150 (MQ=255)
                                                                                                                                tACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTg                      <  1:4478299/150‑1 (MQ=255)
                                                                                                                                  ccGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGcc                    <  1:2973022/150‑1 (MQ=255)
                                                                                                                                  ccGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGcc                    <  2:371538/150‑1 (MQ=255)
                                                                                                                                   cGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcg                   >  1:4044478/1‑150 (MQ=255)
                                                                                                                                         gacgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggcggg             >  1:1744085/1‑150 (MQ=255)
                                                                                                                                         gacgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggcggg             >  1:3282446/1‑150 (MQ=255)
                                                                                                                                          acgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcgggg            <  1:2482105/150‑1 (MQ=255)
                                                                                                                                           cgcgGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcggggc           >  1:3657074/1‑150 (MQ=255)
                                                                                                                                             cggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTCTTCTGGTGCCGGGCGGGGCGa         >  2:3545286/1‑150 (MQ=255)
                                                                                                                                             cggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGa         <  2:54165/150‑1 (MQ=255)
                                                                                                                                             cggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGa         <  2:54430/150‑1 (MQ=255)
                                                                                                                                                 gATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCt     >  1:4637854/1‑150 (MQ=255)
                                                                                                                                                 gATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCt     >  2:410236/1‑150 (MQ=255)
                                                                                                                                                    cGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTgcc  >  1:1684866/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTGCC  >  CP009273/279186‑279483

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: