Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,033,606 G→C 15.1% P216P (CCG→CCC yedY → membrane‑anchored, periplasmic TMAO, DMSO reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,033,6060GC15.1% 168.3 / 3.1 68P216P (CCG→CCCyedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base C (8/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAAC  >  CP009273/2033463‑2033752
                                                                                                                                               |                                                                                                                                                  
gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTcaca                                                                                                                                              <  1:1052979/150‑1 (MQ=255)
gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTcaca                                                                                                                                              <  1:911102/150‑1 (MQ=255)
 aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCacac                                                                                                                                             <  1:782710/150‑1 (MQ=255)
   ccTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGAAGGGATTCGTTTCCACGAAGCCAAGCATCCCCCCCCACtg                                                                                                                                           >  2:331252/1‑150 (MQ=255)
       ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATTTCGAAGGATTTCGTCTCGACGAAGCAAAGCATCCGCTCACCCttatga                                                                                                                                       >  1:935682/1‑150 (MQ=255)
       ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgat                                                                                                                                         <  1:723532/148‑1 (MQ=255)
         aGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATGTTGAAGGGTTTTGTCTCGACGAAGCAAAGGATTCCCCCACACtgatgaac                                                                                                                                     >  1:378423/1‑148 (MQ=255)
           cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGt                                                                                                                                   <  1:407026/150‑1 (MQ=255)
              tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTATTTCGAAGGATTGCGGCTCGACGAAGCAAAGCAACCCCTCACCCTGATGACCGGAgg                                                                                                                                >  2:603428/1‑150 (MQ=255)
                cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgtg                                                                                                                              <  2:530458/150‑1 (MQ=255)
                cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                                                                               <  2:1169338/149‑1 (MQ=255)
                cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                                                                               <  2:809158/149‑1 (MQ=255)
                  aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGtt                                                                                                                            <  1:102510/150‑1 (MQ=255)
                  aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCCCTATTTCGAAGGATTGCGGCTCGGCGAAGCAAAGCATCCGCTCCCACTGATGACCGGAGGGGtt                                                                                                                            >  1:1066631/1‑150 (MQ=255)
                     tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTAt                                                                                                                         <  1:100793/150‑1 (MQ=255)
                       tGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATgg                                                                                                                       <  1:158672/150‑1 (MQ=255)
                         tCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCa                                                                                                                     <  2:921912/150‑1 (MQ=255)
                               tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGc                                                                                                                 <  1:942704/148‑1 (MQ=255)
                                    aCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCg                                                                                                          <  2:608718/150‑1 (MQ=255)
                                       aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCCTTTTGTCGAAGGAATTCGTCCCGACGAAGCAATGCATCCCCTCCCACTGATGAGCGTAGGGGTTTATAGCAAGGCGGTACCGccc                                                                                                       >  1:1087974/1‑149 (MQ=255)
                                          tATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCGAAATATCCTTATGTCGAAGGATTGGGTCCCGACGAAGCAAAGCATCCGCCCACCCCGATGACCGTAGGGGTTTATGGCAAGGCGGTTCCGCccca                                                                                                     >  1:639800/1‑149 (MQ=255)
                                           aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCCCACACTGATGACCGTAGGGGTTTATGGCAAGGCGTTACCGCCACaa                                                                                                    >  1:8802/1‑149 (MQ=255)
                                           aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTTAAATATTCTTATGTTGTAGGGTTTCGTCTCCACGGAGCAAAGCATCCGCTCACACTGATGACCGCAGGGGGTTATTGCAAGGCGTTACCGCCACaa                                                                                                    >  2:402720/1‑149 (MQ=255)
                                               aCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCCTATTTCGGAGGATTTGGGCTTGACGCAACAAAGCATTCCCCCCCACCGATGACCCTAGGGGGTTATTGCAAAGCGGTACCCCCCCAAAAtag                                                                                               >  1:488677/1‑148 (MQ=255)
                                               aCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATgg                                                                                               <  2:1019296/150‑1 (MQ=255)
                                                        gATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATTTCGAAGGGTTTCGTCTCGACGAAACAATGCATCCCCCCACACCGATGACCGGAGGGGGTTATGGCAAAGCGTTACCGCCACAAAAAGGCGCGCCggg                                                                                      >  1:134051/1‑149 (MQ=255)
                                                         aTGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAGGGCTTTGCCGCCACAAAATGGCGCGCCGGTg                                                                                     <  1:877759/150‑1 (MQ=255)
                                                           gccaggccagCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATTCTTATGTCGGAGGATTGCGTCTCGGCGGAGGAAAGCATCCGCTCCCACTGATGAACGTAGGGGGTTTTTGGAAGGGGGTACCGCCACAAAAAGGCGGGCCGGGGCg                                                                                   >  2:78343/1‑150 (MQ=255)
                                                            ccaggccagCAGGACCGCTTTATCGGCGGTGGGGCGAAAAATCCTTTTTTCGGAGGGGTGCGTCTCGACGGAGCAAAGCATCCCCCCCCACCGATGGCCGTAGGGGGTTTTTGCAAGGGGGTTCCCCCCCAAAAAAGGGCGCCGGGGCgg                                                                                  >  2:225487/1‑149 (MQ=255)
                                                             caggccagCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGAc                                                                                 <  1:405722/150‑1 (MQ=255)
                                                              aggccagCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTATTTCGAAGGGTTGCGGCTCCACGAAGCAATGCATCCCCCCCCACCGATGACCGGAGGGGGTTATTGCAAAGCGGTACCCCCCCCAAAAGGCGCGCCgggggggcc                                                                                >  1:957756/1‑143 (MQ=255)
                                                                  cagcagGACCGCTTTATCGGCGGTGGGCTGAAATTTCCTTATGTCGAAGGATTGCGTCTCGAGGAAGCAATGCATCCGCTCACACGGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAtt                                                                            <  2:101632/150‑1 (MQ=255)
                                                                      agGCCCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCCCCCTGATGCCCGGAGGTGTTTATGGAAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGc                                                                        <  2:74037/150‑1 (MQ=255)
                                                                      agGACCGCTTTATCGGCGGTGGGGTGAAATATTCTTATTTCGAAGGATTTCGTCTCCACGAAGCAATGCATCCGCTCACACCGCTGACCGTAGGTGTTTTTTGCAAAGCGGTACCCCCCCAAAAAAGCGCGCCCGGGCGGCTGGTtttgg                                                                        >  2:769636/1‑149 (MQ=255)
                                                                                tgtCGGCGATGGTCTAAAATATCCTTATGTAGAGGGATGCCGTCTCGGCGAAGCAATGAATCCCCTCCCACTTATGACCGTGGGTTTTTATGGAAAGCGTTTGCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat                                                               <  1:991163/147‑1 (MQ=255)
                                                                                 aTCGGCGGTGGGCTAAAATATCTTTATGTCGAAGGATTGCGTCCGGACGAAGCAATGCACCCGCTCCCACTGATGACCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata                                                              <  1:645901/149‑1 (MQ=255)
                                                                                      gggggggCTAAAATATCCTTATGTAGAAGGATTGCGTCCCGACGAAGCAATGCTTCCGCTCCCACTGATGCCCGTAGGTTTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt                                                        <  2:843539/146‑1 (MQ=255)
                                                                                        gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCttt                                                       >  1:20114/1‑149 (MQ=255)
                                                                                        gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCttt                                                       >  2:9680/1‑149 (MQ=255)
                                                                                           ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg                                                   >  2:498305/1‑150 (MQ=255)
                                                                                           ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg                                                   >  2:454236/1‑150 (MQ=255)
                                                                                           ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg                                                   >  2:421465/1‑150 (MQ=255)
                                                                                             ctaaaaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAggg                                                 <  2:993983/147‑1 (MQ=255)
                                                                                             cTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg                                                  >  2:396283/1‑149 (MQ=255)
                                                                                                aaaTATCCTTATGTGGAGGGATGGCGTCCCGACGAAACAATCCACCCGCTCACCCGGTAGACGGTAGGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAtt                                              <  2:204757/150‑1 (MQ=255)
                                                                                                aaaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAt                                               >  1:171969/1‑149 (MQ=255)
                                                                                                aaaTATCCTTATGCCGAAGGATTGCGTCTCGACGAAGCAATGCATCCCCTCACACTGATGACCGTAGTTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAt                                               <  2:425025/149‑1 (MQ=255)
                                                                                                 aaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAtt                                              >  2:496771/1‑149 (MQ=255)
                                                                                                   tatCCTTATGTCGCAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGTTGACCGTAGGTGTTTATGGCAAGTCGTTACCTCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTaaa                                           >  2:846584/1‑150 (MQ=255)
                                                                                                         tATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTGGGTGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGAt                                      <  2:1087974/149‑1 (MQ=255)
                                                                                                          acggggAAGGTTTGCGTCTCGAGGAGGAAAGGCATCCGCTCACACTGATGACCGTGGGTTTTTATGGCAAGGCATTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATa                                     <  2:260242/144‑1 (MQ=255)
                                                                                                          aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGCGGAAATATGGCTTTAAAGGGATTAAATCGATAg                                    >  2:593351/1‑150 (MQ=255)
                                                                                                           tGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACGGTAGGTGTTTATGGCAAGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                   <  1:774458/150‑1 (MQ=255)
                                                                                                           tGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                   >  1:286070/1‑150 (MQ=255)
                                                                                                             tCGAAGGATTGCGTCCGGACAAAGCAATGCACCCCCTCACACTGATGCCCGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCa                                 <  2:8802/150‑1 (MQ=255)
                                                                                                             tCGAAGGATCGCGTGCGGCCAAAGCAATGTGTCCGCTCGCGCGGACGACCGTAGATTTTTATGAAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTc                                  <  2:5394/149‑1 (MQ=255)
                                                                                                                  ggATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTAt                             >  1:63810/1‑149 (MQ=255)
                                                                                                                   gATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTAtt                            >  1:566861/1‑149 (MQ=255)
                                                                                                                   gATTGCGTCTCGACGAAGAAATGCATCCGCTCACCCTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTa                           <  2:716322/150‑1 (MQ=255)
                                                                                                                    aTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTAtt                            >  2:63829/1‑148 (MQ=255)
                                                                                                                      tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                        >  1:485035/1‑150 (MQ=255)
                                                                                                                          tctcGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGAAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGAc                    <  1:78343/150‑1 (MQ=255)
                                                                                                                             cgacgaAGCAATGCATCCGCTCACACTGATGCCCGTAGGTTTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcg                 <  1:331252/150‑1 (MQ=255)
                                                                                                                             cgacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTTTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcg                 <  1:1139693/150‑1 (MQ=255)
                                                                                                                             cgacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcg                 >  2:263246/1‑150 (MQ=255)
                                                                                                                                    gCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGt          >  2:945739/1‑150 (MQ=255)
                                                                                                                                      aaTGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAACTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTcc        >  2:463680/1‑150 (MQ=255)
                                                                                                                                            tCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac  >  2:660833/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAAC  >  CP009273/2033463‑2033752

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