Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,133,918 C→A 12.4% E81* (GAA→TAA)  hybE ← hydrogenase 2‑specific chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,133,9180CA12.4% 778.9 / 56.1 224E81* (GAA→TAA) hybEhydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base C (98/98);  new base A (11/17);  total (109/115)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

GCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGAAAAGAGA  >  CP009273/3133769‑3134067
                                                                                                                                                     |                                                                                                                                                     
 cAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTtca                                                                                                                                                      >  2:474144/1‑150 (MQ=255)
 cAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTtca                                                                                                                                                      >  1:1379080/1‑150 (MQ=255)
       tCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTcactcact                                                                                                                                                >  1:4126212/1‑150 (MQ=255)
        ccGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACtt                                                                                                                                               <  2:4266890/150‑1 (MQ=255)
         cGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACttt                                                                                                                                              <  2:2864213/150‑1 (MQ=255)
         cGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACttt                                                                                                                                              <  2:2864186/150‑1 (MQ=255)
           tCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTgc                                                                                                                                            >  1:2300187/1‑150 (MQ=255)
                cGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGc                                                                                                                                       <  1:2454181/150‑1 (MQ=255)
                   tGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGc                                                                                                                                    <  2:1146121/150‑1 (MQ=255)
                   tGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGc                                                                                                                                    <  1:314977/150‑1 (MQ=255)
                      cccTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGAGCAGCGGCCag                                                                                                                                 >  2:1754971/1‑150 (MQ=255)
                       ccTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCaga                                                                                                                                >  1:1277632/1‑150 (MQ=255)
                        ctctTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCagag                                                                                                                               >  2:1550974/1‑150 (MQ=255)
                        ctctTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCagag                                                                                                                               >  1:949529/1‑150 (MQ=255)
                           ttCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTg                                                                                                                            <  1:1209604/150‑1 (MQ=255)
                           ttCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTg                                                                                                                            <  1:2337765/150‑1 (MQ=255)
                           ttCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTg                                                                                                                            <  2:1520538/150‑1 (MQ=255)
                            tCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGa                                                                                                                           >  2:4412329/1‑150 (MQ=255)
                              aaTCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATc                                                                                                                         >  2:1322375/1‑150 (MQ=255)
                               aTCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATAGTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATcc                                                                                                                        <  1:4492268/150‑1 (MQ=255)
                                   aCATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATccgggc                                                                                                                    <  1:625347/150‑1 (MQ=255)
                                   aCATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATccgggc                                                                                                                    >  2:4219255/1‑150 (MQ=255)
                                    cATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAACGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATccgggcc                                                                                                                   >  1:4078371/1‑150 (MQ=255)
                                     aTGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATccgggccg                                                                                                                  >  2:3063881/1‑150 (MQ=255)
                                     aTGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATccgggccg                                                                                                                  <  2:925254/150‑1 (MQ=255)
                                      tGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATccgggccgg                                                                                                                 <  1:2836717/150‑1 (MQ=255)
                                         tgtgCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGa                                                                                                              >  2:1435993/1‑150 (MQ=255)
                                          gtgCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGaa                                                                                                             >  2:1387037/1‑150 (MQ=255)
                                           tgCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAg                                                                                                            >  1:1973899/1‑150 (MQ=255)
                                              gAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATa                                                                                                         >  2:2952667/1‑150 (MQ=255)
                                               aaaGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATaa                                                                                                        >  1:1303545/1‑150 (MQ=255)
                                                aaGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAAc                                                                                                       >  1:3876169/1‑150 (MQ=255)
                                                aaGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAAc                                                                                                       <  1:3977787/150‑1 (MQ=255)
                                                 aGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACt                                                                                                      <  1:1754876/150‑1 (MQ=255)
                                                  gcggcgACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTg                                                                                                     <  1:445406/150‑1 (MQ=255)
                                                    ggcgACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCAGAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCa                                                                                                   >  2:3174115/1‑150 (MQ=255)
                                                     gcgACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCAc                                                                                                  <  1:785306/150‑1 (MQ=255)
                                                     gcgACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCAc                                                                                                  <  1:3107964/150‑1 (MQ=255)
                                                         cATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAg                                                                                              >  2:4322845/1‑150 (MQ=255)
                                                          aTCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGc                                                                                             >  1:4141962/1‑150 (MQ=255)
                                                           tCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCa                                                                                            <  1:3520043/150‑1 (MQ=255)
                                                             aGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATc                                                                                          <  2:2644099/150‑1 (MQ=255)
                                                             aGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTACACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATc                                                                                          >  2:1114716/1‑150 (MQ=255)
                                                               cGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCa                                                                                        <  1:3301382/150‑1 (MQ=255)
                                                                 aaCAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACg                                                                                      <  2:4005842/150‑1 (MQ=255)
                                                                    aGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgggg                                                                                   <  2:3659581/150‑1 (MQ=255)
                                                                    aGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgggg                                                                                   >  1:3202097/1‑150 (MQ=255)
                                                                    aGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgggg                                                                                   >  1:1502099/1‑150 (MQ=255)
                                                                    aGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgggg                                                                                   >  2:605768/1‑150 (MQ=255)
                                                                       agagaTATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGa                                                                                >  1:2850006/1‑150 (MQ=255)
                                                                        gagaTATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGAt                                                                               <  2:4357697/150‑1 (MQ=255)
                                                                           atatTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcac                                                                            >  2:2516318/1‑150 (MQ=255)
                                                                                 cGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAccc                                                                      >  2:3997546/1‑150 (MQ=255)
                                                                                  gAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcccc                                                                     >  1:589595/1‑150 (MQ=255)
                                                                                  gAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcccc                                                                     >  1:3065115/1‑150 (MQ=255)
                                                                                    aaCACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGt                                                                   >  1:4531010/1‑150 (MQ=255)
                                                                                       acCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCa                                                                <  2:2737718/150‑1 (MQ=255)
                                                                                       acCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCa                                                                >  2:1281483/1‑150 (MQ=255)
                                                                                         cGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACt                                                              <  2:923994/150‑1 (MQ=255)
                                                                                         cGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACt                                                              <  1:4399937/150‑1 (MQ=255)
                                                                                          gTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTg                                                             >  1:2990461/1‑150 (MQ=255)
                                                                                          gTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTg                                                             >  2:2140768/1‑150 (MQ=255)
                                                                                          gTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTg                                                             >  2:2692089/1‑150 (MQ=255)
                                                                                           tCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGa                                                            <  1:1088567/150‑1 (MQ=255)
                                                                                             cAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGAcc                                                          <  2:2552267/150‑1 (MQ=255)
                                                                                               gTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCtt                                                        <  2:2423855/150‑1 (MQ=255)
                                                                                                ttCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                                                       <  2:1610145/150‑1 (MQ=255)
                                                                                                ttCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                                                       <  2:3368036/150‑1 (MQ=255)
                                                                                                ttCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                                                       <  1:2041121/150‑1 (MQ=255)
                                                                                                    ccAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAAc                                                   <  2:2495452/150‑1 (MQ=255)
                                                                                                     cAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACa                                                  <  1:2946580/150‑1 (MQ=255)
                                                                                                      aaCGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAg                                                 <  1:1358165/150‑1 (MQ=255)
                                                                                                        cGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCg                                               >  2:2267986/1‑150 (MQ=255)
                                                                                                            aaaGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaa                                           >  1:2115301/1‑150 (MQ=255)
                                                                                                            aaaGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaa                                           >  1:4187178/1‑150 (MQ=255)
                                                                                                              aGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAAt                                         <  2:147087/150‑1 (MQ=255)
                                                                                                              aGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAAt                                         <  2:3515456/150‑1 (MQ=255)
                                                                                                                gTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCa                                       >  2:4138484/1‑150 (MQ=255)
                                                                                                                 tCATAGTGCCATAAGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCTCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAg                                      <  2:1061153/150‑1 (MQ=255)
                                                                                                                    tAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGaaa                                   >  1:4183704/1‑150 (MQ=255)
                                                                                                                    tAGTGCCATACGGCAGTTGCAGACCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGaaa                                   <  2:4487094/150‑1 (MQ=255)
                                                                                                                     aGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAAc                                  <  1:1825337/150‑1 (MQ=255)
                                                                                                                      gTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACa                                 >  1:1027450/1‑150 (MQ=255)
                                                                                                                       tGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACAt                                <  2:4077599/150‑1 (MQ=255)
                                                                                                                          cATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATaca                             >  2:622345/1‑150 (MQ=255)
                                                                                                                            tACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACAcc                           <  1:1436218/150‑1 (MQ=255)
                                                                                                                                   gTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTg                    >  2:3887037/1‑150 (MQ=255)
                                                                                                                                    ttGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGa                   <  2:721640/150‑1 (MQ=255)
                                                                                                                                     tGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGaa                  >  1:1075589/1‑150 (MQ=255)
                                                                                                                                         gACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACTCACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGAt              <  2:751901/150‑1 (MQ=255)
                                                                                                                                         gACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGAt              <  1:4356325/150‑1 (MQ=255)
                                                                                                                                         gACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGAt              <  2:1687115/150‑1 (MQ=255)
                                                                                                                                         gACCAATTTTTTCACTCACTGTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGAt              <  2:15168/150‑1 (MQ=255)
                                                                                                                                          aCCAATTTTTTAACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATg             <  2:1062173/150‑1 (MQ=255)
                                                                                                                                               ttttttCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGaa        >  1:748924/1‑150 (MQ=255)
                                                                                                                                               ttttttAACTCACGTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACTGGGTGATCACACACCCCGTCCCCTGACATTCTAACAGCGTAAAATCAGAAACATACACCGGCATTGAATGATGCAGaa        <  2:2122253/150‑1 (MQ=255)
                                                                                                                                                 ttttCACTGAATTTGCGCAGCGGCCAGTGTTGATCCGGGTCGGGGAAGATTACTGCACTCAGCATCCCCGGGGTGATCACACACCCCTTCCACTGACGTTCGAACAGCGTATATTCAGAAACATTCACCGGCATTGAAGGATGCAGaaaa      <  2:2318778/150‑1 (MQ=255)
                                                                                                                                                  tttCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGAAAag     <  2:2782845/150‑1 (MQ=255)
                                                                                                                                                    tcactcacTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGAAAagag   <  2:1202340/150‑1 (MQ=255)
                                                                                                                                                     cactcactTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGAAAagaga  <  2:2116394/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACCGGCATTGAAGGATGCAGAAAAGAGA  >  CP009273/3133769‑3134067

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: