Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 754,572 G→T 83.1% E137D (GAG→GAT sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273754,5720GT83.1% 650.9 / 151.8 271E137D (GAG→GATsucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (17/29);  new base T (122/103);  total (139/132)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 274 total aligned reads displayed.

GCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAG  >  CP009273/754425‑754719
                                                                                                                                                   |                                                                                                                                                   
  aGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATAcc                                                                                                                                                  >  2:4183664/1‑149 (MQ=255)
        aaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAc                                                                                                                                            >  1:172361/1‑149 (MQ=255)
              tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGt                                                                                                                                      <  1:138810/149‑1 (MQ=255)
               gCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGtt                                                                                                                                     <  2:1984427/149‑1 (MQ=255)
                 aGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                   >  1:430178/1‑149 (MQ=255)
                 aGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                   >  1:1749490/1‑149 (MQ=255)
                  gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAt                                                                                                                                  <  2:2695686/149‑1 (MQ=255)
                  gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAt                                                                                                                                  >  1:3229164/1‑149 (MQ=255)
                    tCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAttt                                                                                                                                >  2:1853879/1‑149 (MQ=255)
                    tCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAttt                                                                                                                                >  1:316953/1‑149 (MQ=255)
                      aTTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAGTCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                              <  2:4592020/149‑1 (MQ=255)
                      aTTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                              >  1:1875184/1‑149 (MQ=255)
                       ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGcc                                                                                                                             <  1:3878217/149‑1 (MQ=255)
                       ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGcc                                                                                                                             <  1:1791332/149‑1 (MQ=255)
                       ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGcc                                                                                                                             <  1:547170/149‑1 (MQ=255)
                          aCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGc                                                                                                                          >  2:2414081/1‑149 (MQ=255)
                          aCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGc                                                                                                                          >  2:3841243/1‑149 (MQ=255)
                          aCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGc                                                                                                                          >  2:372481/1‑149 (MQ=255)
                              aTACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCa                                                                                                                      >  1:2781583/1‑149 (MQ=255)
                                aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCaaa                                                                                                                    >  1:3502150/1‑149 (MQ=255)
                                aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCaaa                                                                                                                    >  1:523480/1‑149 (MQ=255)
                                 ccGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAg                                                                                                                   >  2:552027/1‑149 (MQ=255)
                                   gCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaa                                                                                                                 <  2:1756794/149‑1 (MQ=255)
                                   gCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaa                                                                                                                 >  2:1009504/1‑149 (MQ=255)
                                       ccGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCa                                                                                                             >  1:3889302/1‑149 (MQ=255)
                                       cagTGGTCTCCAGCATGCGAATCTCGATGCGCTGGGACTGTGGCAGCAAGATTAAGTGGCCGTTCTGTATCCGTCTTTCCTCGATCTGACCGTAGCAGTCTTCAAGGTGACCTTGTACTTCTGTTCTTTTGCCTGCGGCTAAGAATCCa                                                                                                             <  2:2731600/147‑1 (MQ=255)
                                         gTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATg                                                                                                           >  2:1150830/1‑149 (MQ=255)
                                         gTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATg                                                                                                           <  2:3268169/149‑1 (MQ=255)
                                            gTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaaa                                                                                                        <  1:344645/149‑1 (MQ=255)
                                            gTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaaa                                                                                                        >  1:4265823/1‑149 (MQ=255)
                                             tCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAAc                                                                                                       <  1:2972713/149‑1 (MQ=255)
                                               aCCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTc                                                                                                     >  1:4660026/1‑149 (MQ=255)
                                                  aGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGc                                                                                                  <  1:1067824/149‑1 (MQ=255)
                                                      tGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                              >  2:82909/1‑149 (MQ=255)
                                                      tGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCCGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                              <  1:2043956/149‑1 (MQ=255)
                                                         gAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctgc                                                                                           <  1:2446263/149‑1 (MQ=255)
                                                           aTCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGctgctg                                                                                         <  1:2211206/149‑1 (MQ=255)
                                                             ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGa                                                                                       <  2:475158/149‑1 (MQ=255)
                                                              tcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGaa                                                                                      <  2:1668520/149‑1 (MQ=255)
                                                              tcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGaa                                                                                      >  2:3869521/1‑149 (MQ=255)
                                                                 aTCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                   >  1:2160476/1‑149 (MQ=255)
                                                                  tCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGccc                                                                                  <  1:4200347/149‑1 (MQ=255)
                                                                  tCCGCTGGGACTGTGGCAGCAAGAGAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATGTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGccc                                                                                  <  2:4392856/149‑1 (MQ=255)
                                                                   ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCt                                                                                 >  2:3073022/1‑149 (MQ=255)
                                                                        gggACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGc                                                                            >  2:3104504/1‑149 (MQ=255)
                                                                        gggACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGc                                                                            >  2:3912034/1‑149 (MQ=255)
                                                                            cTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAc                                                                        <  2:1275003/149‑1 (MQ=255)
                                                                             tgtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAcc                                                                       >  1:743270/1‑149 (MQ=255)
                                                                             tgtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAcc                                                                       >  2:2858201/1‑149 (MQ=255)
                                                                              gtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCt                                                                      >  1:3705375/1‑149 (MQ=255)
                                                                               tgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTa                                                                     <  1:4115177/149‑1 (MQ=255)
                                                                                ggCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTAc                                                                    >  1:2296758/1‑149 (MQ=255)
                                                                                 gcagcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACt                                                                   <  2:2259274/149‑1 (MQ=255)
                                                                                   agcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGc                                                                 <  1:588709/149‑1 (MQ=255)
                                                                                     caAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCgg                                                               >  1:301310/1‑149 (MQ=255)
                                                                                        gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGccc                                                            >  1:2861814/1‑149 (MQ=255)
                                                                                          tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCTGGATCCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAACCCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGa                                                          >  1:1027282/1‑149 (MQ=255)
                                                                                          tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGa                                                          >  1:4352560/1‑149 (MQ=255)
                                                                                           acaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAt                                                         <  1:2012968/147‑1 (MQ=255)
                                                                                           aaaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAt                                                         <  1:2109176/149‑1 (MQ=255)
                                                                                            aaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAtt                                                        <  2:3155384/149‑1 (MQ=255)
                                                                                             aGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                       <  1:3833292/149‑1 (MQ=255)
                                                                                             aGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                       >  2:2133593/1‑149 (MQ=255)
                                                                                              gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTgg                                                      <  1:996308/149‑1 (MQ=255)
                                                                                               tGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGt                                                     <  1:1187959/149‑1 (MQ=255)
                                                                                                 gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGc                                                   >  1:3605009/1‑149 (MQ=255)
                                                                                                  ccGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                  <  2:4159080/149‑1 (MQ=255)
                                                                                                    gATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                                <  1:4627041/149‑1 (MQ=255)
                                                                                                       cTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGta                                             <  2:1130923/149‑1 (MQ=255)
                                                                                                           aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                         >  1:489504/1‑149 (MQ=255)
                                                                                                              cGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGcac                                      >  2:1466591/1‑149 (MQ=255)
                                                                                                               gTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGcaca                                     >  1:1330320/1‑149 (MQ=255)
                                                                                                                tCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAt                                    <  2:3004338/149‑1 (MQ=255)
                                                                                                                 cTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAtt                                   <  1:923965/149‑1 (MQ=255)
                                                                                                                  tttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                                  <  1:3354389/149‑1 (MQ=255)
                                                                                                                  tttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                                  <  1:4183302/149‑1 (MQ=255)
                                                                                                                   ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                                 <  2:409365/149‑1 (MQ=255)
                                                                                                                         gATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAc                           <  2:3446439/149‑1 (MQ=255)
                                                                                                                          aTCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                          >  2:3603057/1‑149 (MQ=255)
                                                                                                                           tCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCg                         >  2:571382/1‑149 (MQ=255)
                                                                                                                           tCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCg                         >  2:3754684/1‑149 (MQ=255)
                                                                                                                            cTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAACCCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCga                        <  2:2612545/149‑1 (MQ=255)
                                                                                                                             tGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                       <  1:2639872/149‑1 (MQ=255)
                                                                                                                              gACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                      >  1:2295097/1‑149 (MQ=255)
                                                                                                                               aCCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                     <  2:4368919/149‑1 (MQ=255)
                                                                                                                                ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                    >  2:3283534/1‑149 (MQ=255)
                                                                                                                                ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                    >  2:3283512/1‑149 (MQ=255)
                                                                                                                                  gAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa                  >  2:625155/1‑149 (MQ=255)
                                                                                                                                   aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                 >  1:1657740/1‑149 (MQ=255)
                                                                                                                                   aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                 <  2:1548687/149‑1 (MQ=255)
                                                                                                                                      cAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGc              <  2:4541581/149‑1 (MQ=255)
                                                                                                                                        gACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg            <  2:1405330/149‑1 (MQ=255)
                                                                                                                                        gACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg            >  1:3900670/1‑149 (MQ=255)
                                                                                                                                        gACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg            >  2:1593044/1‑149 (MQ=255)
                                                                                                                                         aCTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg           <  1:2386462/149‑1 (MQ=255)
                                                                                                                                            tCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc        <  1:4329519/149‑1 (MQ=255)
                                                                                                                                             ccAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATcc       <  2:3935339/149‑1 (MQ=255)
                                                                                                                                              cAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCa      >  2:912788/1‑149 (MQ=255)
                                                                                                                                               aGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCaa     >  1:2662650/1‑149 (MQ=255)
                                                                                                                                                  agaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCTGCACCGAAGAAAAACGCTGGATCCAACAg  >  2:4268190/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAG  >  CP009273/754425‑754719

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: