Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,604,731 C→T 8.5% R223W (CGG→TGG)  hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,604,7310CT 8.5% 1031.6 / 62.3 295R223W (CGG→TGG) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (135/134);  new base T (10/15);  total (146/149)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

GGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGA  >  CP009273/2604582‑2604880
                                                                                                                                                     |                                                                                                                                                     
gggTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATc                                                                                                                                                       >  1:3430437/1‑150 (MQ=255)
  gTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCgg                                                                                                                                                     >  1:4336600/1‑150 (MQ=255)
   tGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCggg                                                                                                                                                    <  1:1153876/150‑1 (MQ=255)
     aCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGtt                                                                                                                                                  <  1:1402366/150‑1 (MQ=255)
     aCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGtt                                                                                                                                                  >  1:1470823/1‑150 (MQ=255)
      cTATGCAACCGTTGTGGCCGCAGATCCCGCCAGCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTa                                                                                                                                                 >  1:2286773/1‑150 (MQ=255)
       tATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTCCAGGAAATCGGGTTaa                                                                                                                                                >  1:1674067/1‑150 (MQ=255)
           cAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACAcc                                                                                                                                            <  1:4842861/150‑1 (MQ=255)
              ccGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTgg                                                                                                                                         <  2:643875/150‑1 (MQ=255)
                   gtgGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTgcg                                                                                                                                    >  1:1032318/1‑150 (MQ=255)
                       ccGCAGATCCCGCCATCACCGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAACTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGCAATTGGGTTAACACCGGGTTGCGCGAt                                                                                                                                >  2:5004886/1‑150 (MQ=255)
                        cGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATg                                                                                                                               <  2:1942987/150‑1 (MQ=255)
                         gCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGc                                                                                                                              >  1:1087691/1‑150 (MQ=255)
                          cAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGcc                                                                                                                             <  1:619866/150‑1 (MQ=255)
                          cAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGcc                                                                                                                             <  2:2234704/150‑1 (MQ=255)
                                 cGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATc                                                                                                                      >  2:4097453/1‑150 (MQ=255)
                                    cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCTATGCCGACGATCCAc                                                                                                                   >  1:22339/1‑150 (MQ=255)
                                    cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                                   <  1:4749100/150‑1 (MQ=255)
                                    cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                                   >  1:4542/1‑150 (MQ=255)
                                    cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                                   >  1:5121197/1‑150 (MQ=255)
                                       cacaGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGtc                                                                                                                <  1:236346/150‑1 (MQ=255)
                                         caGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGtctc                                                                                                              <  2:1315693/150‑1 (MQ=255)
                                          aGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCa                                                                                                             <  1:4018157/150‑1 (MQ=255)
                                             gTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATa                                                                                                          <  2:3301226/150‑1 (MQ=255)
                                              tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAg                                                                                                         <  2:4282128/150‑1 (MQ=255)
                                              tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAg                                                                                                         <  2:4282553/150‑1 (MQ=255)
                                               aTCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGt                                                                                                        >  1:2832097/1‑150 (MQ=255)
                                                tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                       <  2:4539129/150‑1 (MQ=255)
                                                tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                       >  1:4982386/1‑150 (MQ=255)
                                                tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                       >  1:1850544/1‑150 (MQ=255)
                                                 cGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTAt                                                                                                      >  1:4713118/1‑150 (MQ=255)
                                                  gCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATc                                                                                                     <  1:2587471/150‑1 (MQ=255)
                                                  gCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATc                                                                                                     >  1:4646637/1‑150 (MQ=255)
                                                    cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                   >  2:266434/1‑150 (MQ=255)
                                                     aTTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCGGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                                  <  2:5066833/150‑1 (MQ=255)
                                                     aTTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                                  >  1:1442027/1‑150 (MQ=255)
                                                       tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                >  1:850919/1‑150 (MQ=255)
                                                       tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                <  2:4922262/150‑1 (MQ=255)
                                                         aGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGc                                                                                              >  2:4787752/1‑150 (MQ=255)
                                                           cGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCaa                                                                                            >  2:3374670/1‑150 (MQ=255)
                                                                aGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                       >  2:4682086/1‑150 (MQ=255)
                                                                   gcggcgGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                                    <  2:1200102/150‑1 (MQ=255)
                                                                     ggcggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcgcg                                                                                  >  2:3930277/1‑150 (MQ=255)
                                                                            ggCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCg                                                                           <  1:1029560/150‑1 (MQ=255)
                                                                            ggCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCg                                                                           >  1:1304241/1‑150 (MQ=255)
                                                                             gCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCgg                                                                          <  1:1747919/150‑1 (MQ=255)
                                                                                  cTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGtt                                                                     >  1:3573199/1‑150 (MQ=255)
                                                                                   tATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTa                                                                    <  1:4650804/150‑1 (MQ=255)
                                                                                   tATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTa                                                                    >  1:3166979/1‑150 (MQ=255)
                                                                                   tATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTa                                                                    <  2:2419677/150‑1 (MQ=255)
                                                                                   tATCGTCAGGGGAGAGAAATTTGCTATCGGCTCCTGCGACAATTAATCGACGATCCTACAGGAAATCGGGTTAACACCTGTTTGCTCGATGCCGACGATCCACGTCTCAATAGTATCGTACAGCAACTCTATAGCGTACTCCGGGGGTTa                                                                    >  2:848855/1‑150 (MQ=255)
                                                                                    aTCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTAc                                                                   <  2:3619070/150‑1 (MQ=255)
                                                                                     tCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACa                                                                  <  1:1903264/150‑1 (MQ=255)
                                                                                     tCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACa                                                                  <  1:2914571/150‑1 (MQ=255)
                                                                                         cAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGAc                                                              >  2:785589/1‑150 (MQ=255)
                                                                                          aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                             <  1:1360005/150‑1 (MQ=255)
                                                                                          aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                             >  1:3863045/1‑150 (MQ=255)
                                                                                            gggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGa                                                           >  1:4045696/1‑150 (MQ=255)
                                                                                             ggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGaa                                                          >  2:2978859/1‑150 (MQ=255)
                                                                                             ggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGaa                                                          >  2:383131/1‑150 (MQ=255)
                                                                                                gagaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                       >  2:2715236/1‑150 (MQ=255)
                                                                                                gagaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                       >  1:4199654/1‑150 (MQ=255)
                                                                                                  gaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTg                                                     <  1:4143166/150‑1 (MQ=255)
                                                                                                   aaaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGc                                                    <  1:2022130/150‑1 (MQ=255)
                                                                                                    aaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCg                                                   >  2:4025767/1‑150 (MQ=255)
                                                                                                        tGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGa                                               >  1:4950102/1‑150 (MQ=255)
                                                                                                        tGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGa                                               <  1:3844560/150‑1 (MQ=255)
                                                                                                           gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                            <  1:2735288/150‑1 (MQ=255)
                                                                                                           gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                            <  2:1718890/150‑1 (MQ=255)
                                                                                                            aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                           <  1:2403574/150‑1 (MQ=255)
                                                                                                             tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTg                                          <  1:1739462/150‑1 (MQ=255)
                                                                                                             tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTg                                          >  2:2557394/1‑150 (MQ=255)
                                                                                                                gCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGttt                                       <  2:4595296/150‑1 (MQ=255)
                                                                                                                 cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                      <  2:3163665/150‑1 (MQ=255)
                                                                                                                 cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                      >  1:1237567/1‑150 (MQ=255)
                                                                                                                 cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                      >  1:3413525/1‑150 (MQ=255)
                                                                                                                 cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                      <  1:4820740/150‑1 (MQ=255)
                                                                                                                 cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                      <  2:1654597/150‑1 (MQ=255)
                                                                                                                  tCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTc                                     <  1:3981102/150‑1 (MQ=255)
                                                                                                                      gcgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGa                                 >  2:4627497/1‑150 (MQ=255)
                                                                                                                          cATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCt                             >  1:2109404/1‑150 (MQ=255)
                                                                                                                           aTTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTg                            <  2:3447285/150‑1 (MQ=255)
                                                                                                                              tAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGa                         <  1:1223531/150‑1 (MQ=255)
                                                                                                                                 gCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaa                      >  2:1910401/1‑150 (MQ=255)
                                                                                                                                 gCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaa                      <  2:1442027/150‑1 (MQ=255)
                                                                                                                                  cgacgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaaa                     >  2:2056168/1‑150 (MQ=255)
                                                                                                                                      gaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGtt                 >  1:3089990/1‑150 (MQ=255)
                                                                                                                                      gaTACTACAGGAAATTGGGTTAACACATGGTTGAGCGATGACGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGAGTACTACGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGTACGCTGCGAAAAAAGtt                 >  2:5107279/1‑150 (MQ=255)
                                                                                                                                       aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt                >  1:4761668/1‑150 (MQ=255)
                                                                                                                                        tCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTg               <  2:3566960/150‑1 (MQ=255)
                                                                                                                                        tCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTg               >  1:4471266/1‑150 (MQ=255)
                                                                                                                                          cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc             >  1:4072798/1‑150 (MQ=255)
                                                                                                                                          cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc             >  1:156284/1‑150 (MQ=255)
                                                                                                                                              aGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCa         <  2:3771751/150‑1 (MQ=255)
                                                                                                                                               ggAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAg        <  2:2620631/150‑1 (MQ=255)
                                                                                                                                                 aaaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTa      <  2:1074601/150‑1 (MQ=255)
                                                                                                                                                  aaTTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg     >  1:3653282/1‑150 (MQ=255)
                                                                                                                                                  aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg     <  2:2842859/150‑1 (MQ=255)
                                                                                                                                                    tCGGGTTAACCCCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcg   <  1:3000403/150‑1 (MQ=255)
                                                                                                                                                     cGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcga  <  2:211162/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGA  >  CP009273/2604582‑2604880

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: