Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 635,464 G→T 16.3% E86* (GAG→TAG)  ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273635,4640GT16.3% 1065.1 / 152.0 337E86* (GAG→TAG) ahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (127/155);  new base T (18/37);  total (145/192)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 338 total aligned reads displayed.

AAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGT  >  CP009273/635317‑635612
                                                                                                                                                   |                                                                                                                                                    
   gAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGtt                                                                                                                                                  >  1:3146060/1‑149 (MQ=255)
    aaCTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTc                                                                                                                                                 >  1:5062005/1‑149 (MQ=255)
     aCTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCa                                                                                                                                                >  2:3304390/1‑149 (MQ=255)
     aCTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCa                                                                                                                                                >  2:2549016/1‑149 (MQ=255)
       tGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCAcc                                                                                                                                              <  2:2594722/149‑1 (MQ=255)
       tGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCAcc                                                                                                                                              >  1:2402778/1‑149 (MQ=255)
        gTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCt                                                                                                                                             <  1:1587959/149‑1 (MQ=255)
              tGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTgg                                                                                                                                       >  2:4202367/1‑149 (MQ=255)
                 aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTAc                                                                                                                                    <  1:658318/149‑1 (MQ=255)
                  aTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACt                                                                                                                                   >  1:4216783/1‑149 (MQ=255)
                   tCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTg                                                                                                                                  <  2:2643949/149‑1 (MQ=255)
                    cGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTgg                                                                                                                                 <  1:3653520/149‑1 (MQ=255)
                       aGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                                                                                                              <  1:4432614/149‑1 (MQ=255)
                        gAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGtt                                                                                                                             >  2:491445/1‑149 (MQ=255)
                        gAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGtt                                                                                                                             >  2:3961525/1‑149 (MQ=255)
                          aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGc                                                                                                                           >  2:151815/1‑149 (MQ=255)
                              tCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCtgtg                                                                                                                       <  1:3917075/149‑1 (MQ=255)
                                aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGa                                                                                                                     <  2:3042573/149‑1 (MQ=255)
                                 gACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAc                                                                                                                    >  2:2925199/1‑149 (MQ=255)
                                  aCAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAcc                                                                                                                   <  2:2354901/149‑1 (MQ=255)
                                     aaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCggt                                                                                                                <  1:460448/149‑1 (MQ=255)
                                      aGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTCCTGGCGTTGCTGTGGACCggtg                                                                                                               >  2:4059632/1‑149 (MQ=255)
                                      aGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCggtg                                                                                                               <  1:3560376/149‑1 (MQ=255)
                                       gTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCggtgg                                                                                                              <  1:3408185/149‑1 (MQ=255)
                                        tCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCggtggt                                                                                                             >  1:1951658/1‑149 (MQ=255)
                                         cACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTc                                                                                                            >  1:259316/1‑149 (MQ=255)
                                         cACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTc                                                                                                            <  1:2321320/149‑1 (MQ=255)
                                          aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCa                                                                                                           <  2:160361/149‑1 (MQ=255)
                                          aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCa                                                                                                           <  2:2763680/149‑1 (MQ=255)
                                              ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCg                                                                                                       <  2:597210/149‑1 (MQ=255)
                                               tAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGt                                                                                                      >  2:323228/1‑149 (MQ=255)
                                                 aagaagATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCg                                                                                                    >  1:214712/1‑149 (MQ=255)
                                                 aagaagATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCg                                                                                                    <  2:791553/149‑1 (MQ=255)
                                                  agaagaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGa                                                                                                   >  2:680148/1‑149 (MQ=255)
                                                    aagaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGaaa                                                                                                 <  2:2303210/149‑1 (MQ=255)
                                                     agaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaag                                                                                                <  2:1009024/149‑1 (MQ=255)
                                                      gaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaaga                                                                                               <  2:2038680/149‑1 (MQ=255)
                                                       aTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaagaa                                                                                              <  2:1243255/149‑1 (MQ=255)
                                                        tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaagaag                                                                                             >  1:804751/1‑149 (MQ=255)
                                                         aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGACAGAAgc                                                                                            >  2:2208422/1‑149 (MQ=255)
                                                         aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAgc                                                                                            >  1:3833798/1‑149 (MQ=255)
                                                          aCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAgcg                                                                                           <  2:4604525/149‑1 (MQ=255)
                                                           cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAgcgc                                                                                          >  1:3027520/1‑149 (MQ=255)
                                                             gCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAg                                                                                        >  2:3862201/1‑149 (MQ=255)
                                                               ttGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGtc                                                                                      >  1:456046/1‑149 (MQ=255)
                                                                tGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGtct                                                                                     <  1:4946357/149‑1 (MQ=255)
                                                                    ggTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTctgc                                                                                 <  2:1718110/149‑1 (MQ=255)
                                                                          tAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGc                                                                           >  2:3599542/1‑149 (MQ=255)
                                                                          tAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGc                                                                           <  2:5166113/149‑1 (MQ=255)
                                                                          tAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGc                                                                           <  1:2000141/149‑1 (MQ=255)
                                                                              ccGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGAt                                                                       <  2:3891107/149‑1 (MQ=255)
                                                                               cGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGAtt                                                                      >  2:3641926/1‑149 (MQ=255)
                                                                               cGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGAtt                                                                      <  2:1800027/149‑1 (MQ=255)
                                                                                gTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTc                                                                     <  1:4293006/149‑1 (MQ=255)
                                                                                    ttCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCa                                                                 <  2:1382895/149‑1 (MQ=255)
                                                                                      ccTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCata                                                               <  2:781444/149‑1 (MQ=255)
                                                                                          aTCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGa                                                           >  1:4017791/1‑149 (MQ=255)
                                                                                              ccAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGt                                                       <  2:4816671/149‑1 (MQ=255)
                                                                                              ccAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGt                                                       >  2:572870/1‑149 (MQ=255)
                                                                                                 aCCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGAt                                                    <  2:1832624/149‑1 (MQ=255)
                                                                                                 aCCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGAt                                                    <  1:4402152/149‑1 (MQ=255)
                                                                                                    cAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGAtttt                                                 >  2:408789/1‑149 (MQ=255)
                                                                                                      ggTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGa                                               <  2:1947162/149‑1 (MQ=255)
                                                                                                      ggTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGa                                               <  1:3569042/149‑1 (MQ=255)
                                                                                                      ggTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGa                                               <  1:1037261/149‑1 (MQ=255)
                                                                                                         tCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAAtt                                            <  1:3304390/149‑1 (MQ=255)
                                                                                                         tCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAAtt                                            <  1:2549011/149‑1 (MQ=255)
                                                                                                         tCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAAtt                                            >  2:1857985/1‑149 (MQ=255)
                                                                                                         tCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAAtt                                            <  2:121556/149‑1 (MQ=255)
                                                                                                          ccAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTc                                           <  2:2415151/149‑1 (MQ=255)
                                                                                                           cAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCg                                          >  2:4495627/1‑149 (MQ=255)
                                                                                                            aaCCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGa                                         >  2:3173248/1‑149 (MQ=255)
                                                                                                              ccAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGaaa                                       <  1:1639165/149‑1 (MQ=255)
                                                                                                               cAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAAc                                      >  1:3811665/1‑149 (MQ=255)
                                                                                                                aGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAAcc                                     <  2:3742/149‑1 (MQ=255)
                                                                                                                aGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAAcc                                     >  2:2577099/1‑149 (MQ=255)
                                                                                                                  gggCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTa                                   <  2:1703117/149‑1 (MQ=255)
                                                                                                                   ggCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAt                                  <  1:3479570/149‑1 (MQ=255)
                                                                                                                    gCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAtt                                 <  1:1504236/149‑1 (MQ=255)
                                                                                                                      cACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTAc                               >  2:3466943/1‑149 (MQ=255)
                                                                                                                      cACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTAc                               >  2:1072218/1‑149 (MQ=255)
                                                                                                                           tttGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGct                          <  2:456068/149‑1 (MQ=255)
                                                                                                                              gCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGctctc                       <  2:1371990/149‑1 (MQ=255)
                                                                                                                              gCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGctctc                       <  1:4855327/149‑1 (MQ=255)
                                                                                                                              gCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGctctc                       <  2:2137008/149‑1 (MQ=255)
                                                                                                                                aGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCtt                     >  1:3572670/1‑149 (MQ=255)
                                                                                                                                  gCTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGc                   >  2:1268935/1‑149 (MQ=255)
                                                                                                                                  gCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTAGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGc                   <  2:852770/149‑1 (MQ=255)
                                                                                                                                   cTCCCCGCTGGGCCACTAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGcc                  >  2:3575084/1‑149 (MQ=255)
                                                                                                                                   cTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGcc                  <  2:4106944/149‑1 (MQ=255)
                                                                                                                                     ccccGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCcac                >  1:1252135/1‑149 (MQ=255)
                                                                                                                                      cccGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCcaca               >  1:2846632/1‑149 (MQ=255)
                                                                                                                                      cccGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCcaca               >  1:1043431/1‑149 (MQ=255)
                                                                                                                                        cGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAAc             <  1:4088760/149‑1 (MQ=255)
                                                                                                                                            gggCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGcc         >  1:2814732/1‑149 (MQ=255)
                                                                                                                                             ggCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGccc        <  1:1450247/149‑1 (MQ=255)
                                                                                                                                                cACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGa     <  2:4402340/149‑1 (MQ=255)
                                                                                                                                                cACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGa     >  1:925315/1‑149 (MQ=255)
                                                                                                                                                 aCGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGAc    >  1:5118971/1‑149 (MQ=255)
                                                                                                                                                   gAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACgt  <  2:2929942/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGT  >  CP009273/635317‑635612

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: