Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,676,832 G→T 81.1% intergenic (‑25/+158) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,676,8320GT81.1% 168.4 / 25.0 86intergenic (‑25/+158)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (11/6);  new base T (40/27);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAA  >  CP009273/3676692‑3676966
                                                                                                                                            |                                                                                                                                      
tCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGt                                                                                                                                      >  1:312589/1‑143 (MQ=255)
         ggTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATgg                                                                                                                       <  2:431863‑M1/149‑3 (MQ=255)
           tttCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCCATGGCTAGCGCGGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGtt                                                                                                                           >  1:64437/1‑143 (MQ=255)
             tCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTAt                                                                                                                         <  1:322702/143‑1 (MQ=255)
               aTTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGGCCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTGATgg                                                                                                                       >  1:375140‑M1/1‑141 (MQ=255)
                tttGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCCGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTATTCATAGGGTGTCCTTTAGTAAAACCACATGTAGTGTAGGGTTAtttgttgcct                                                                                                                >  2:253055/1‑140 (MQ=255)
                tttGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggg                                                                                                                      <  1:361457‑M1/143‑4 (MQ=255)
                   gCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcg                                                                                                             <  2:285384‑M1/149‑13 (MQ=255)
                   gCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcg                                                                                                             <  2:292261‑M1/149‑13 (MQ=255)
                     tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGATTATggggtatctgcggg                                                                                                           >  2:419698‑M1/1‑135 (MQ=255)
                       gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggg                                                                                                         >  2:332932‑M1/1‑133 (MQ=255)
                           aTTTCAGGATAGAAATGGCCAAGGAGAATAACAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTACACAGAGAGGTTTTCATAGGGTGTCCTTTAGTCACAACACATGTCTTGTAGGGTTATgggatacctgcggggtacct                                                                                                     >  2:226445‑M1/1‑129 (MQ=255)
                            tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtaccgccggggtacctg                                                                                                    >  2:107847‑M1/1‑128 (MQ=255)
                                aGGATAGAAATTGCCAAGGAGAATACCAATGGCTATCGCGGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggc                                                                                                <  2:314008‑M1/149‑26 (MQ=255)
                                aGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggc                                                                                                <  2:255351‑M1/149‑26 (MQ=255)
                                  gATCGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGCGGTTTTCATAGGGTGTCCTTTAGTAAAACCGCATGTCTTGTAAGGTTATggggtacctgcgggggacctgcggcct                                                                                              >  2:260735‑M1/1‑122 (MQ=255)
                                    tAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcg                                                                                                  >  1:255189‑M1/1‑120 (MQ=255)
                                      gAAGTGGCCAAGGAGGATGCCAATGGCTATCGCTGTGAGGACCTGAAAGTAAAGGCTTTGAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggc                                                                                                <  1:342567‑M1/143‑26 (MQ=255)
                                      gAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggc                                                                                                <  1:162150‑M1/143‑26 (MQ=255)
                                       aaaTGGCCAAGGAGAATACCCATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGGTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtccctgcagagtacctgcggccttaaca                                                                                         >  2:29057‑M1/1‑117 (MQ=255)
                                       aaaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaa                                                                                         <  2:312089‑M1/149‑33 (MQ=255)
                                       aaaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaa                                                                                         <  2:82320‑M1/149‑33 (MQ=255)
                                        aaTGGGCAAGGAGAATACCAATGTCCATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAacac                                                                                              <  1:219591/143‑1 (MQ=255)
                                        aaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaat                                                                                        >  2:105930‑M1/1‑116 (MQ=255)
                                         aTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGGTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                             >  1:321012/1‑143 (MQ=255)
                                          tGGCTAAGGAGAATACCAAGGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataa                                                                                      <  2:173205‑M1/149‑36 (MQ=255)
                                           ggCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataac                                                                                     <  2:370301‑M1/149‑37 (MQ=255)
                                            gCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaa                                                                                          <  1:321659‑M1/143‑32 (MQ=255)
                                             ccAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaa                                                                                         >  1:89659‑M1/1‑111 (MQ=255)
                                              cAAGGGGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacc                                                                                  <  2:333649‑M1/149‑40 (MQ=255)
                                              cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTGTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaat                                                                                        >  1:20820‑M1/1‑110 (MQ=255)
                                               aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacaccc                                                                                 >  2:262159‑M1/1‑109 (MQ=255)
                                               aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacaccc                                                                                 >  2:262191‑M1/1‑109 (MQ=255)
                                               aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaata                                                                                       <  1:348074‑M1/143‑35 (MQ=255)
                                               aaGGAGAATACCAATGGCTATCGCAGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTGAGGTTATggggtacctgcggggtacctgcggccttaaaataacaccc                                                                                 >  2:321339‑M1/1‑109 (MQ=255)
                                                aGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTACGGTTATggggtacctgcggggtacgtgcggcagtaaaataacaccca                                                                                >  2:118733‑M1/1‑108 (MQ=255)
                                                  gagaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGGGGGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccaga                                                                              >  2:230492‑M1/1‑106 (MQ=255)
                                                   agaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagac                                                                             >  2:392252‑M1/1‑105 (MQ=255)
                                                      aTACCTATGGCTATCGCGGTCAGGACCTGAAAGTAAAGGCTTTGAAACAGAGAGGTTTTCATAGTGTGTCCTTTAGTAAAACCACATGTCTTCTAAGGTTATggtgtacctgcggggtacctgcggccttaaaataacaccca                                                                                <  1:60537‑M1/143‑42 (MQ=255)
                                                       tACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtaccttcggggtaccggcggccttaaaatcacaccccg                                                                               >  1:25580‑M1/1‑101 (MQ=255)
                                                          cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                      >  2:61632/1‑149 (MQ=255)
                                                           aaTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACACCACAGaaa                                                                           >  1:206996/1‑143 (MQ=255)
                                                                    cggtgttgggacctgaaaGTAAAGGCTTGTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggttcctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgta                                                            <  2:204202‑M1/141‑62 (MQ=255)
                                                                      cTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGAt                                                                <  1:247604/143‑1 (MQ=255)
                                                                         tCAGGACCTGAAAGTAAAGGCTTTTACACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACCAGTCTTGTAAGGTTATggggtacctgcggggtacctacggccttaaaataacaccaagcgcacatcagagtaaagt                                                             >  1:352165‑M1/1‑83 (MQ=255)
                                                                           aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTGATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggat                                                     >  2:343002‑M1/1‑81 (MQ=255)
                                                                                cTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataa                                                >  2:141511‑M1/1‑76 (MQ=255)
                                                                                cTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaaggtacctggatcataa                                                >  2:28717‑M1/1‑76 (MQ=255)
                                                                                cTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTACTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctgga                                                      >  1:149343‑M1/1‑76 (MQ=255)
                                                                                 tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcaccgaaatgtacctggatcataat                                               >  2:192725‑M1/1‑75 (MQ=255)
                                                                                  ggaagtAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAATACCACATGTCTTGTAAGGTTATggggtacctgcgggggacctgcggccttaaaataacacccagacaacatcacagaaatgtacctcgact                                                    <  1:256779‑M1/141‑70 (MQ=255)
                                                                                  gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggacc                                                    >  1:77875‑M1/1‑74 (MQ=255)
                                                                                      gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccgtaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaa                                          >  2:256226‑M1/1‑70 (MQ=255)
                                                                                       tAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggaccataat                                               >  1:137799‑M1/1‑69 (MQ=255)
                                                                                       tAAAGGCTTTACAACAGAGGGGTTTTCATAGGGGGTCTTTTAGTAAAACCACAGGTCGTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataat                                               <  1:150990‑M1/143‑75 (MQ=255)
                                                                                            gCTTTTAAACAGAGAGGTTTTCATAGGTTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaa                                    <  2:365590‑M1/149‑86 (MQ=255)
                                                                                             cTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAACACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaag                                   <  2:70823‑M1/149‑87 (MQ=255)
                                                                                                ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaa                                      >  1:9689‑M1/1‑60 (MQ=255)
                                                                                                ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacagcacagaaatgtacctggttcataattgaaacaaa                                      >  1:9652‑M1/1‑60 (MQ=255)
                                                                                                  aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctgggtcataattaagaccaaaaggttaa                              >  2:169536‑M1/1‑58 (MQ=255)
                                                                                                  aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggaggacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctgcatcataattgaaacaaaaa                                    >  1:111875‑M1/1‑58 (MQ=255)
                                                                                                   aaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaag                                   <  1:226583‑M1/143‑87 (MQ=255)
                                                                                                         gagGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagttt                       <  2:108832‑M1/149‑99 (MQ=255)
                                                                                                          agGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagacatgtacctggatcataattgaaacaaaacggttaaaaagtttg                      >  2:398582‑M1/1‑50 (MQ=255)
                                                                                                           ggTTTTCATAGGGTGTCCTTTAGGAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttga                     <  2:384402‑M1/149‑101 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgag                    >  2:201816‑M1/1‑48 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgag                    <  2:35496‑M1/149‑102 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaag                          >  1:373379‑M1/1‑48 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaaaatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgag                    <  2:85916‑M1/149‑102 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAGTTGAAACAAAAAGGTAAAAAAGTTTGAGCTGGc                          >  1:194460/1‑143 (MQ=255)
                                                                                                             ttttCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagt                         >  1:291066‑M1/1‑47 (MQ=255)
                                                                                                              tttCATAGGGTGTCCTTTAGTAAAACCAAATGTCTTGTAAGGTTATgtggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataataaaaaaaaaaaggttaaaaagtt                        <  1:296449‑M1/143‑98 (MQ=255)
                                                                                                               ttCATAGGGTGTCCTTTCGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggactgaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggataaaaagttt                       >  1:154537‑M1/1‑45 (MQ=255)
                                                                                                               ttCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaacagttg                       >  1:154543‑M1/1‑45 (MQ=255)
                                                                                                                tCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGATGGCTCGc                      >  1:215710/1‑143 (MQ=255)
                                                                                                                   tAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtaccgggatcataattgaaacaaaaaggttaaaaagtttgagctggctc             >  2:382934‑M1/1‑41 (MQ=255)
                                                                                                                    aGGGTCTCCTTTAGTAAAACTACAGGTCCTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGc            <  2:11181/149‑1 (MQ=255)
                                                                                                                     ggggggcctttgggAAAACCACATGTCTTGTAAGGTTATggggtaccggcgggggacctgcggccttaaaataacacccaggcaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctg                 >  1:342672‑M1/7‑39 (MQ=255)
                                                                                                                      ggTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctgg                >  1:181397‑M1/1‑38 (MQ=255)
                                                                                                                            cctgaggaaaacaacatgTTTCGTAAGGTTAGggggtacctgcggcgtacctgcggccttaaaataacaccaagacaacatcacagaaatgtacctggatcataattgaacaagaaatgttaaaaagtttgagctggctcgca          <  1:376598‑M1/136‑112 (MQ=255)
                                                                                                                             tttAGTAAAACCACATGTCTTGTAAGGTTCTggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcaccgaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcgc         >  1:75090‑M1/1‑31 (MQ=255)
                                                                                                                             tttAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggataaaaagtttgagctggctcgcac         >  1:36797‑M1/1‑31 (MQ=255)
                                                                                                                              ttAGTAATACCCCATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaaccag  <  2:73503‑M1/149‑120 (MQ=255)
                                                                                                                                aGTAAAACCACATGTCTTGTAAGGTTATggcgtaactagggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattacaaaaaacaggttaaaaagtttgagctggccagcccaca      >  1:83109‑M1/1‑28 (MQ=255)
                                                                                                                                 gTAAATCCACAGGTCTTGTAACGTTATggggtacctgcgggggacctgcggccttaaaataacacccagacaacatcacagaaacgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaac     <  1:69211‑M1/143‑117 (MQ=255)
                                                                                                                                  tAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggataaaaagtttgagctgactcgcacaaacc    >  1:55194‑M1/1‑26 (MQ=255)
                                                                                                                                    aaaCCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaaccag  >  1:225949‑M1/1‑24 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAA  >  CP009273/3676692‑3676966

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