Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,274,926 A→C 100% I72M (ATT→ATG yjcH ← DUF485 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,274,9260AC97.7% 1092.7 / 4.2 276I72M (ATT→ATGyjcHDUF485 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (130/139);  total (133/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 279 total aligned reads displayed.

CGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCG  >  CP009273/4274779‑4275072
                                                                                                                                                   |                                                                                                                                                  
cgTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCAt                                                                                                                                                  <  2:1571463/150‑1 (MQ=255)
  ttAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCa                                                                                                                                                <  2:87611/150‑1 (MQ=255)
           aaaGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAacac                                                                                                                                       >  1:385348/1‑150 (MQ=255)
            aaGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACAcc                                                                                                                                      >  2:2220016/1‑150 (MQ=255)
                gAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAAtt                                                                                                                                  >  1:1106183/1‑150 (MQ=255)
                 aGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTg                                                                                                                                 <  1:3504163/150‑1 (MQ=255)
                 aGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTg                                                                                                                                 >  2:4151532/1‑150 (MQ=255)
                   tgtgGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGa                                                                                                                               <  2:3353908/150‑1 (MQ=255)
                       gcggcAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATAc                                                                                                                           <  2:1467653/150‑1 (MQ=255)
                         ggcAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACcg                                                                                                                         >  1:1998273/1‑150 (MQ=255)
                         ggcAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACcg                                                                                                                         <  2:1661381/150‑1 (MQ=255)
                           cAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACcgcg                                                                                                                       <  1:165965/150‑1 (MQ=255)
                                gcCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGa                                                                                                                  <  2:3637105/150‑1 (MQ=255)
                                 ccGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGAc                                                                                                                 >  1:2216203/1‑150 (MQ=255)
                                   gTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGc                                                                                                               >  1:233600/1‑150 (MQ=255)
                                            ctTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCg                                                                                                      <  1:3632086/150‑1 (MQ=255)
                                              tttCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGa                                                                                                    <  2:3802101/150‑1 (MQ=255)
                                              tttCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGa                                                                                                    >  2:3220873/1‑150 (MQ=255)
                                                tCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGAtt                                                                                                  >  2:2428446/1‑150 (MQ=255)
                                                 cATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTc                                                                                                 >  2:1333108/1‑150 (MQ=255)
                                                 cATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTc                                                                                                 <  1:2391575/150‑1 (MQ=255)
                                                     atgCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAgcg                                                                                             >  2:2012499/1‑150 (MQ=255)
                                                     atgCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAgcg                                                                                             <  1:1228549/150‑1 (MQ=255)
                                                       gCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAgcggc                                                                                           >  2:4469194/1‑150 (MQ=255)
                                                                  aTGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGAGGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagcca                                                                                <  2:2835419/150‑1 (MQ=255)
                                                                   tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                               >  1:2558078/1‑150 (MQ=255)
                                                                   tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                               >  1:3169590/1‑150 (MQ=255)
                                                                   tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                               <  2:4279941/150‑1 (MQ=255)
                                                                   tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                               <  2:454222/150‑1 (MQ=255)
                                                                   tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                               >  2:805047/1‑150 (MQ=255)
                                                                     cAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCagccagcc                                                                             >  2:3711136/1‑150 (MQ=255)
                                                                     cAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCagccagcc                                                                             <  1:4308079/150‑1 (MQ=255)
                                                                      aGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCg                                                                            >  1:1174659/1‑150 (MQ=255)
                                                                      aGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCg                                                                            >  2:718692/1‑150 (MQ=255)
                                                                      aGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCg                                                                            >  1:645831/1‑150 (MQ=255)
                                                                       ggACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCgg                                                                           >  1:3205236/1‑150 (MQ=255)
                                                                         aCTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCACTTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGgc                                                                         <  2:1065030/150‑1 (MQ=255)
                                                                             cATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCACTTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGa                                                                     <  1:919855/150‑1 (MQ=255)
                                                                              attattAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGaa                                                                    <  2:940456/150‑1 (MQ=255)
                                                                              attattAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGaa                                                                    >  1:53174/1‑150 (MQ=255)
                                                                               ttattaAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATTCCGGTGATAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCAGGGCGCGAAg                                                                   <  1:4203861/150‑1 (MQ=255)
                                                                                 attaAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGc                                                                 >  2:3593552/1‑150 (MQ=255)
                                                                                  ttaAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCg                                                                <  1:3982383/150‑1 (MQ=255)
                                                                                   taAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGa                                                               <  1:4426973/150‑1 (MQ=255)
                                                                                   taAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGa                                                               >  1:3749948/1‑150 (MQ=255)
                                                                                    aaGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGAt                                                              <  1:1298362/150‑1 (MQ=255)
                                                                                     aGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATc                                                             >  1:3441921/1‑150 (MQ=255)
                                                                                       aCGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAg                                                           <  1:4637247/150‑1 (MQ=255)
                                                                                        cGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGt                                                          <  2:3031219/150‑1 (MQ=255)
                                                                                        cGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGt                                                          <  2:4629652/150‑1 (MQ=255)
                                                                                         ggTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTa                                                         <  1:232394/150‑1 (MQ=255)
                                                                                         ggTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTa                                                         >  2:4529503/1‑150 (MQ=255)
                                                                                         ggTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTa                                                         >  2:3927598/1‑150 (MQ=255)
                                                                                          gTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTaa                                                        <  1:2257500/150‑1 (MQ=255)
                                                                                           tCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTaaa                                                       <  1:839249/150‑1 (MQ=255)
                                                                                               aTTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATaa                                                   <  1:2063685/150‑1 (MQ=255)
                                                                                                   gCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGcc                                               >  2:1136591/1‑150 (MQ=255)
                                                                                                    ccGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCg                                              >  2:1776615/1‑150 (MQ=255)
                                                                                                     cGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGa                                             >  2:409317/1‑150 (MQ=255)
                                                                                                      gTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGat                                            <  1:978540/150‑1 (MQ=255)
                                                                                                      gTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGat                                            >  2:551020/1‑150 (MQ=255)
                                                                                                      gTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGat                                            <  1:15499/150‑1 (MQ=255)
                                                                                                        tCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGatat                                          <  2:2520357/150‑1 (MQ=255)
                                                                                                          gcccgccAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATaa                                        <  1:2118137/150‑1 (MQ=255)
                                                                                                            ccgccAGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAAc                                      >  2:3430020/1‑150 (MQ=255)
                                                                                                                cAGATGTAGATACCGGTGAGAACAAAGTAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGcc                                  >  2:3320040/1‑150 (MQ=255)
                                                                                                                 aGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGACGGCGATCAGTAAAATAAAGCCGATATAAACTGCCa                                 <  1:1389715/150‑1 (MQ=255)
                                                                                                                 aGATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCa                                 >  1:317638/1‑150 (MQ=255)
                                                                                                                  gATGTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGGGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAg                                >  2:3777889/1‑150 (MQ=255)
                                                                                                                     gTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCat                             <  2:4254656/150‑1 (MQ=255)
                                                                                                                     gTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCat                             <  2:2252685/150‑1 (MQ=255)
                                                                                                                     gTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCat                             >  2:3123256/1‑150 (MQ=255)
                                                                                                                     gTAGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCac                             <  2:2125958/150‑2 (MQ=255)
                                                                                                                       aGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATCTAAACTGCCAGCataa                           >  1:4596989/1‑150 (MQ=255)
                                                                                                                       aGATACCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCataa                           >  2:853984/1‑150 (MQ=255)
                                                                                                                           aCCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATaataat                       <  1:4151532/150‑1 (MQ=255)
                                                                                                                           aCCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATaataat                       <  2:2886569/150‑1 (MQ=255)
                                                                                                                           aCCGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATaataat                       <  1:866519/150‑1 (MQ=255)
                                                                                                                            ccGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAATGAAAGCCGATATAAACTGCCAGCATAATAATc                      <  2:3670509/150‑1 (MQ=255)
                                                                                                                            ccGGTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATc                      <  2:2353816/150‑1 (MQ=255)
                                                                                                                               gTGAGAACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGAc                   <  1:4076837/150‑1 (MQ=255)
                                                                                                                                agagaACCAAGGAGAACAACATCACTCCAACACCAAATGGAATACCGCGTGTGACGCTGGTGTTCTGATACAGCGTCTTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATAAAAACAGCCAGCAAAAAAATCTACa                  >  2:2228744/2‑150 (MQ=255)
                                                                                                                                  agaACAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAgg                >  1:3716169/1‑150 (MQ=255)
                                                                                                                                    aaCAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAt              <  2:4028118/150‑1 (MQ=255)
                                                                                                                                     aCAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtg             >  1:998965/1‑150 (MQ=255)
                                                                                                                                     aCAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtg             <  2:338421/150‑1 (MQ=255)
                                                                                                                                     aCAAAGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtg             <  2:4003130/150‑1 (MQ=255)
                                                                                                                                     aCAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtg             <  1:3897690/150‑1 (MQ=255)
                                                                                                                                       aaaGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggt           <  1:722055/150‑1 (MQ=255)
                                                                                                                                       aaaGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggt           >  2:2496170/1‑150 (MQ=255)
                                                                                                                                        aaGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggtg          <  1:2996206/150‑1 (MQ=255)
                                                                                                                                         aGGAGATCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggtgg         >  2:3521424/1‑150 (MQ=255)
                                                                                                                                           gagaTCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCa       >  1:1647889/1‑150 (MQ=255)
                                                                                                                                            agaTCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGGGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCaa      <  2:355515/150‑1 (MQ=255)
                                                                                                                                              aTCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAAc    <  1:1843116/150‑1 (MQ=255)
                                                                                                                                               tCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAAcc   <  2:1666109/150‑1 (MQ=255)
                                                                                                                                               tCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAAcc   >  2:4176895/1‑150 (MQ=255)
                                                                                                                                               tCACCATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAAcc   <  1:3292518/150‑1 (MQ=255)
                                                                                                                                                caccATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCg  <  1:3345552/150‑1 (MQ=255)
                                                                                                                                                caccATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCg  >  2:3516923/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCG  >  CP009273/4274779‑4275072

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: