Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 754,572 G→T 100% 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,5720GT99.4% 506.9 / ‑4.0 158E137D (GAG→GATsucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (82/75);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 158 total aligned reads displayed.

AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAA  >  CP009273/754427‑754716
                                                                                                                                                 |                                                                                                                                                
   ttAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCtt                                                                                                                                             >  2:2122857/1‑148 (MQ=255)
   ttAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCtt                                                                                                                                             >  1:1197876/1‑148 (MQ=255)
     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCa                                                                                                                                           <  2:2676725/148‑1 (MQ=255)
     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCa                                                                                                                                           <  2:1393255/148‑1 (MQ=255)
      aaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCaa                                                                                                                                          <  1:2441182/148‑1 (MQ=255)
       aGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAc                                                                                                                                         >  1:1349963/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    >  1:2100543/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    >  1:2100876/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    <  1:1938493/148‑1 (MQ=255)
                gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                <  2:1405018/148‑1 (MQ=255)
                gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                >  2:1952097/1‑148 (MQ=255)
                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                           >  2:245653/1‑148 (MQ=255)
                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                           <  2:1080454/148‑1 (MQ=255)
                        aCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAg                                                                                                                        >  2:1099492/1‑148 (MQ=255)
                           cATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCgg                                                                                                                     >  2:2815478/1‑148 (MQ=255)
                              aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGcca                                                                                                                  >  2:1764781/1‑146 (MQ=255)
                              aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCaa                                                                                                                  >  1:60220/1‑148 (MQ=255)
                                 gCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGa                                                                                                               >  1:1283971/1‑148 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  1:1094899/148‑1 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             >  1:436856/1‑148 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  2:248523/148‑1 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  2:1918461/148‑1 (MQ=255)
                                     ccGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAAcc                                                                                                           <  2:174358/148‑1 (MQ=255)
                                         ggTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGa                                                                                                       <  1:2059521/148‑1 (MQ=255)
                                          gTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaa                                                                                                      <  1:1068495/148‑1 (MQ=255)
                                            cACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAAc                                                                                                    >  1:371045/1‑148 (MQ=255)
                                              ccAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTc                                                                                                  >  1:2764555/1‑148 (MQ=255)
                                              ccAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTc                                                                                                  >  1:3213059/1‑148 (MQ=255)
                                                aGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCgg                                                                                                <  2:458382/148‑1 (MQ=255)
                                                  cATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCg                                                                                              <  1:284670/148‑1 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           >  2:2270640/1‑148 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           >  2:145309/1‑148 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           <  2:1068067/148‑1 (MQ=255)
                                                       gAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctg                                                                                         <  2:2674649/148‑1 (MQ=255)
                                                        aaTCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCCAAGTAACCATGAAACTCGGCGAgctgc                                                                                        >  1:2689005/1‑148 (MQ=255)
                                                           ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTgg                                                                                     >  2:457742/1‑148 (MQ=255)
                                                           ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTgg                                                                                     >  2:431021/1‑148 (MQ=255)
                                                              gATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAg                                                                                  <  2:3197164/148‑1 (MQ=255)
                                                               aTCCGCTGGGACTGTGTCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGc                                                                                 <  1:2087117/148‑1 (MQ=255)
                                                                tCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                >  1:866295/1‑148 (MQ=255)
                                                                 ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGACACTCGGCGAGCTGCTGGAAGccc                                                                               >  2:3212343/1‑148 (MQ=255)
                                                                 ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGccc                                                                               <  1:2664452/148‑1 (MQ=255)
                                                                  cGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCt                                                                              >  2:68603/1‑148 (MQ=255)
                                                                       ggACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGc                                                                         >  2:509817/1‑148 (MQ=255)
                                                                           tgtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAc                                                                     >  1:2850958/1‑148 (MQ=255)
                                                                            gtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAcc                                                                    <  1:3276540/148‑1 (MQ=255)
                                                                               gcagcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTAc                                                                 >  1:1613760/1‑148 (MQ=255)
                                                                                  gcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGc                                                              <  2:2224751/148‑1 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGTTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1373733/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  2:741316/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          <  2:2047873/148‑1 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1771202/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1706348/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  2:2369901/1‑148 (MQ=255)
                                                                                       aTAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGccc                                                         >  1:119720/1‑148 (MQ=255)
                                                                                        tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCg                                                        <  1:3108776/148‑1 (MQ=255)
                                                                                        tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCg                                                        <  1:3108778/148‑1 (MQ=255)
                                                                                          aaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAt                                                      >  2:737865/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    <  1:837656/148‑1 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  2:1624899/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  2:568780/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  1:2767004/1‑148 (MQ=255)
                                                                                             tGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTgg                                                   >  1:2163460/1‑148 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTg                                                 <  2:920614/148‑1 (MQ=255)
                                                                                                ccGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGc                                                >  1:2637750/1‑148 (MQ=255)
                                                                                                  gATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCg                                              <  1:3186084/148‑1 (MQ=255)
                                                                                                   aTCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                             >  2:2006858/1‑148 (MQ=255)
                                                                                                   aTCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                             >  2:338715/1‑148 (MQ=255)
                                                                                                    tCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                            >  2:2692731/1‑148 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       >  2:1867605/1‑148 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       <  2:736452/148‑1 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       >  2:2400062/1‑148 (MQ=255)
                                                                                                          tCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                      >  2:3268624/1‑148 (MQ=255)
                                                                                                          tCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGAGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                      >  2:605387/1‑148 (MQ=255)
                                                                                                           ccGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                     <  2:2596878/148‑1 (MQ=255)
                                                                                                           ccGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                     <  2:2596885/148‑1 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               >  2:1293343/1‑148 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               <  2:1813196/148‑1 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               <  2:1813226/148‑1 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                              <  1:2495920/148‑1 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                              <  1:1154378/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCATGATACCTTCAACGTCGGTTCTTTTGACAGCGGCAAAGAAACCATGAAACTCGGCGAGCTTCTGGAAGACCTCAAGCAAACCTACTGCGGCCCGTTTGGTGCCGAGTATATGGACATTACCa                            <  2:2516814/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            <  2:2673322/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            >  1:1655725/1‑148 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            >  1:3293892/1‑148 (MQ=255)
                                                                                                                        aTCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCTGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCTAGCTGCTGGAAGCCCTCACGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAc                        >  2:3033541/1‑148 (MQ=255)
                                                                                                                             aCCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                   >  1:69351/1‑148 (MQ=255)
                                                                                                                              ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                  >  2:1431864/1‑148 (MQ=255)
                                                                                                                                 aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa               <  1:3196302/148‑1 (MQ=255)
                                                                                                                                 aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa               >  2:1002110/1‑148 (MQ=255)
                                                                                                                                  aGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa              <  1:114506/148‑1 (MQ=255)
                                                                                                                                    cAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg            <  2:1197876/148‑1 (MQ=255)
                                                                                                                                    cAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg            <  2:1483320/148‑1 (MQ=255)
                                                                                                                                      gACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCt          <  1:68603/148‑1 (MQ=255)
                                                                                                                                       aCTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg         >  1:2127034/1‑148 (MQ=255)
                                                                                                                                       aCTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACAGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg         >  2:1275017/1‑148 (MQ=255)
                                                                                                                                         ttCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGa       >  1:933138/1‑148 (MQ=255)
                                                                                                                                           ccAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc     <  1:769982/148‑1 (MQ=255)
                                                                                                                                           ccAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc     <  1:2562586/148‑1 (MQ=255)
                                                                                                                                              ggATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCaa  <  2:297055/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAA  >  CP009273/754427‑754716

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: