Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,421,172 A→C 13.1% D33E (GAT→GAG argT ← lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,421,1720AC13.1% 176.0 / 16.8 80D33E (GAT→GAGargTlysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (38/30);  new base C (6/5);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATTGTGTT  >  CP009273/2421025‑2421317
                                                                                                                                                   |                                                                                                                                                 
tttCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTAt                                                                                                                                                  >  1:676225/1‑149 (MQ=255)
 ttCAGGGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATc                                                                                                                                                 <  1:55243/149‑1 (MQ=255)
  tCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGGACAGGTCGCCTGCATCCGTTTGCACATCTCGATACCGGGATCGATATCAAAGCCAACAAAATCACCGTTAGCATCTTTCGAGAAGGACGGGGCGTAGGGGGTATCg                                                                                                                                                >  2:44242/1‑149 (MQ=255)
    aGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCGTCTTTCGATGAGAACGGTGCGTAGGTGGTGTCGGt                                                                                                                                              >  1:98142/1‑149 (MQ=255)
      tggggggAGCAGCGCGTCAAAGTCACTGGCATCCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                                                                                            <  1:895802/146‑1 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCAGTGGCAACCCAGGTACATTTCAGCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTc                                                                                                                                            >  1:693216/1‑149 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                                                                                            <  1:1023180/149‑1 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGCAGGTGGTATCGGtac                                                                                                                                            >  2:1120233/1‑147 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGGTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                                                                                            <  1:937341/149‑1 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCACTGGCAAACCAGGTACATTTCACCTGGTTCCGTTTGCACATCTCGTTACCGAGCTCGGTATCACAGCCAACAAAATAACCGTTGGGAACTTTCGATGCGAACGGAGCGCAGGGGGTATCGGGTc                                                                                                                                            >  2:386650/1‑149 (MQ=255)
      tGAGGGGATCAGCGCGTCAAAGTCAATGGCAACCCAGGTACATTTCAGCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATGTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTc                                                                                                                                            >  1:693499/1‑149 (MQ=255)
      tGAGGGGATCAGCGAATCAAACTCTCTGCCAATCTGATTACCTTTCACCTGCATGCTTTTTCACATCTCTTTACCTTTATCTATATCAAAGTCAACAAGATCGCCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                                                                                            <  2:107296/149‑1 (MQ=255)
      ggAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                                                                                            <  2:1227480/148‑1 (MQ=255)
       aaGGGTATCACCGCGCCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTcc                                                                                                                                           <  2:403153/148‑1 (MQ=255)
             aTCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACg                                                                                                                                     >  1:901100/1‑149 (MQ=255)
              tCAGCTCGTCAAAGTCACTGGCCACCCAGGTACATTTCACCGGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGAACCGATAcca                                                                                                                                    <  1:139450/149‑3 (MQ=255)
              tCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGAACCGATAcca                                                                                                                                    <  1:1079334/149‑3 (MQ=255)
                aGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGGGCGTAGGTGGTATCGGTTCCGATACGTAc                                                                                                                                  >  1:470894/1‑149 (MQ=255)
                   gcgTCAAAGTCACTGGCAACCCAGGTACATTTCACCAGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTAAGCATCTTACGATGAGAACGGTGCGTAGGTGGTATCAGTTCCGCAACGTACCgc                                                                                                                               >  2:284046/1‑148 (MQ=255)
                   gcgTCAAACTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCCACAAAATCACCTTTAGCATCTTTCGATGAGAACGGCGCGTAGGTGGTATCGGTTCCGGTACGTACCGt                                                                                                                               >  2:334973/1‑149 (MQ=255)
                    cgTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTCCCGtc                                                                                                                              >  2:1042258/1‑149 (MQ=255)
                    cgTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGtc                                                                                                                              >  1:1201171/1‑149 (MQ=255)
                      tCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGtctc                                                                                                                            >  1:750625/1‑149 (MQ=255)
                      tCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGATATCGGTTCCGATACGTACCGtctc                                                                                                                            >  1:756146/1‑149 (MQ=255)
                       ccaattCACTGGCAACCCAGGTACATCTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTcc                                                                                                                           <  2:415703/144‑1 (MQ=255)
                          aGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGt                                                                                                                        <  1:857762/149‑1 (MQ=255)
                            tCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAg                                                                                                                      <  1:1176850/149‑1 (MQ=255)
                                tGGCAACCCATGTACATTTCACTTGCCTCCGTGTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGcgcc                                                                                                                  <  2:707605/149‑1 (MQ=255)
                                 ggCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGcgccg                                                                                                                 <  1:742105/149‑1 (MQ=255)
                                 ggCAACCCAGCTACATTTCACCGGCACCCGTTTGCACATCACGGTATCGAGATCGATATCAAAGCCAACAAAATCACCTTGAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGcgccg                                                                                                                 <  1:139621/149‑1 (MQ=255)
                                        ccggTACATTTCACCTGCCTCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCt                                                                                                          <  1:603704/147‑1 (MQ=255)
                                        cAGGTACATTTCCCCTGCCTCCGTTTGCACATCTCGTCACCGCGACCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGCTACTAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCt                                                                                                          <  1:416890/149‑1 (MQ=255)
                                        cAGGTACATTTCACCTGCATCCGTTTGCCACTCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCt                                                                                                          <  1:567137/149‑1 (MQ=255)
                                        cAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCt                                                                                                          <  1:584541/149‑1 (MQ=255)
                                              cagttCCCCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGc                                                                                                    <  2:946340/146‑1 (MQ=255)
                                                tttCACCTCCATCCGCTTACACATCTCGTTACCGAGATCGATACCAAAGCCAACAAAATCCCCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCg                                                                                                  <  1:180694/149‑1 (MQ=255)
                                                  tCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATAAGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCAGCt                                                                                                >  2:301041/1‑149 (MQ=255)
                                                   cACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCACTTTTCGCTGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCCCCGGAACCGCCGCATAGCTGGAAGCCGccg                                                                                               >  2:194956/1‑147 (MQ=255)
                                                          aTCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGgagag                                                                                        >  1:671168/1‑149 (MQ=255)
                                                            gcggttGCACATCTCGTTACCGAGAACGATATGCAAGCCAACAACATCACCTTTAGGACCGTTCGATGAGCACGCGTCGTAGGTGAAATCGGCTCCCATACGCTCCGTATCCCGTCGCGCTGCGTAGCTGGATGCCGCGGTGGAGAGAc                                                                                      <  2:112190/145‑1 (MQ=255)
                                                               tttGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTGGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTAACGTCACCGCCACCCCCGCACCGCTGCAACCCGCCGTGGAGAGACCgc                                                                                   >  2:259015/1‑148 (MQ=255)
                                                               tttGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGGGAGACCgc                                                                                   >  1:382829/1‑148 (MQ=255)
                                                                  gCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAACCACAGTTAGCATCTTTCGATGAGAACAGGTCGCAGGTGCTAGCGGTTCCGTTACGTTCCGTCTCCGGTAGCGCCCCAAAGCTTGAATCCGCGGTGGAGAGGccaccta                                                                                >  1:609345/1‑144 (MQ=255)
                                                                   cactcctcGTTACCGAGATCGAGATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTaa                                                                               <  2:21912/144‑1 (MQ=255)
                                                                   cacaTCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTaa                                                                               <  1:943788/149‑1 (MQ=255)
                                                                     caTCTCTTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACa                                                                             <  2:1130528/149‑1 (MQ=255)
                                                                       tctcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa                                                                           <  2:693216/149‑1 (MQ=255)
                                                                       tctcGTCACCGAGACCGATATCAAAGCCAACACAATCACCTGTAGCATCTTTCGATCAGAACGGTGCGTAGGTGTTATCTGTTCCGATACGTACCGTCGCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa                                                                           <  1:161352/149‑1 (MQ=255)
                                                                       tccctttACCGAGATCGATATCAAAGCCAACAACATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa                                                                           <  2:693499/144‑1 (MQ=255)
                                                                              aCCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTCCCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACaggg                                                                    >  1:724783/1‑147 (MQ=255)
                                                                                cGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCg                                                                  >  1:301999/1‑149 (MQ=255)
                                                                                cGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCg                                                                  >  1:970572/1‑149 (MQ=255)
                                                                                cGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCg                                                                  >  1:136716/1‑149 (MQ=255)
                                                                                cGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCAACTAACAAAGACAGAGCg                                                                  >  1:1073365/1‑149 (MQ=255)
                                                                                    aTCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTAAACAGCTGGCATTGTATCAGCTACAAACGTTTCCCGGTAGCGCCGAATCTTTTGAGGTCGCTGTGGGGTGGCCGACTAACAAAGACAGAGCGAGaa                                                              >  1:601168/1‑149 (MQ=255)
                                                                                     tCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAAt                                                             <  1:940846/149‑1 (MQ=255)
                                                                                           tCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGGATCGActt                                                       >  1:1126740/1‑149 (MQ=255)
                                                                                           tCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGAGGAGAGACCGACTAACAAAGCCAAAGCGAGAATCGActt                                                       >  2:350246/1‑149 (MQ=255)
                                                                                           tCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGCAACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCCGTGGAGAGACCGACTAACAAAGCCAGACCGAGAATCgccct                                                       >  1:749669/1‑145 (MQ=255)
                                                                                           tCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGCCAGAGCGAGAATCGActt                                                       >  1:826688/1‑149 (MQ=255)
                                                                                           tCAAAGCCAACAAAATCACCTTTAGCATCATTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGActt                                                       >  1:751134/1‑149 (MQ=255)
                                                                                                gCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTCTCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTca                                                  >  1:266132/1‑149 (MQ=255)
                                                                                                   aacaaAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcatac                                               >  1:663680/1‑149 (MQ=255)
                                                                                                     caaAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcatacat                                             >  2:1049566/1‑149 (MQ=255)
                                                                                                       aaaTCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATaa                                           >  1:1200358/1‑149 (MQ=255)
                                                                                                       aaaTCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCGCCGGTAGCGCCGCATAGCTGGCAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGGGAATCGACTTCTTCATACATaa                                           >  1:150109/1‑149 (MQ=255)
                                                                                                          tCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGGGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCAACCATAAccc                                        >  2:373903/1‑149 (MQ=255)
                                                                                                          tCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGCCCGACTAACAAAGACAGAGCGACACTCGACTTCTTCATACATAAccc                                        >  1:249862/1‑149 (MQ=255)
                                                                                                          tCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAccc                                        >  2:554264/1‑149 (MQ=255)
                                                                                                                   gCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaa                               <  1:719122/149‑1 (MQ=255)
                                                                                                                    cATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaaa                              >  1:1022205/1‑149 (MQ=255)
                                                                                                                    cATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaaa                              <  2:901100/149‑1 (MQ=255)
                                                                                                                    cATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGATCCGATACGTACCGTCTCCGGGAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGTCCGACGCACAAAGACAGAGCGAGAATAGCCTTCTCAAAACGAAAACCTCAAGTGaaa                              >  2:679390/1‑149 (MQ=255)
                                                                                                                       cTTCCGATGAGAACGGTGCGTAGGTGGTAGCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATg                           <  1:1006235/149‑1 (MQ=255)
                                                                                                                        tttCGATGAGAACGGTGCGTAGGTGGTAACGGTTGCGATACGTACCGTCTCCGATAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGGGAAAATgg                          >  1:683462/1‑149 (MQ=255)
                                                                                                                            gaggagaACGGTGCGTATGTGGGATGGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTTTGGAGAGACCGACTATCAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCttt                      <  2:394643/146‑1 (MQ=255)
                                                                                                                                 gaacgggtgggtggtggtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCCGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAt                 <  2:1211644/137‑1 (MQ=255)
                                                                                                                                   aCGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAACCTCAAGTGAAAATGGCTTTTTTAtgt               >  1:677655/1‑149 (MQ=255)
                                                                                                                                       tGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTtgtg           >  1:292563/1‑149 (MQ=255)
                                                                                                                                                ggtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAAGGGCTTTTTTATGTTGTGTAttgtgtt  <  2:124036/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATTGTGTT  >  CP009273/2421025‑2421317

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: