Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,138 C→A 5.1% R454S (CGC→AGC)  metG → methionyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,189,1380CA 5.1% 711.7 / 9.4 184R454S (CGC→AGC) metGmethionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (69/105);  new base A (5/5);  total (74/110)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

TTTGACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTG  >  CP009273/2188991‑2189285
                                                                                                                                                   |                                                                                                                                                   
tttGACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAAcc                                                                                                                                                     >  1:463108/1‑148 (MQ=255)
 ttGACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCg                                                                                                                                                    <  2:2803727/148‑1 (MQ=255)
   gACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCt                                                                                                                                                  <  1:1232444/148‑1 (MQ=255)
   gACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCt                                                                                                                                                  <  1:50052/148‑1 (MQ=255)
     cGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTAt                                                                                                                                                <  2:2720427/148‑1 (MQ=255)
     cGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTAt                                                                                                                                                <  2:58595/148‑1 (MQ=255)
       gCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGt                                                                                                                                              <  2:587111/148‑1 (MQ=255)
       gCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAAACTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGACGTGCTCTAAATCATGGCGCTGGCTTATCTGTCTAAAAGCTAttt                                                                                                                                              >  2:1714090/1‑146 (MQ=255)
         gTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCg                                                                                                                                            <  2:3247038/148‑1 (MQ=255)
           gCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGAt                                                                                                                                          <  2:1387328/148‑1 (MQ=255)
           gCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGAt                                                                                                                                          <  2:1387339/148‑1 (MQ=255)
             tGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGa                                                                                                                                        >  2:1162175/1‑148 (MQ=255)
               gCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAAc                                                                                                                                      <  1:1757280/148‑1 (MQ=255)
               gCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAAc                                                                                                                                      <  1:1330508/148‑1 (MQ=255)
                  aGCTAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGACGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAgg                                                                                                                                   >  2:2583730/1‑148 (MQ=255)
                  aGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAgg                                                                                                                                   <  2:1482500/148‑1 (MQ=255)
                   gCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGc                                                                                                                                  <  2:2943600/148‑1 (MQ=255)
                    cGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCt                                                                                                                                 >  1:387035/1‑148 (MQ=255)
                    cGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCt                                                                                                                                 >  1:2992778/1‑148 (MQ=255)
                     gAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTc                                                                                                                                <  2:2200576/148‑1 (MQ=255)
                       aCTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCg                                                                                                                              <  1:3525731/148‑1 (MQ=255)
                           gCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATAATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtggg                                                                                                                          <  1:210605/148‑1 (MQ=255)
                                    cAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGa                                                                                                                 <  1:2719136/148‑1 (MQ=255)
                                    cAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGa                                                                                                                 <  1:2496313/148‑1 (MQ=255)
                                      gTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGaaa                                                                                                               >  2:3306452/1‑148 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGAGGGTGGTGGCGAAAc                                                                                                              >  1:2743059/1‑148 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAAc                                                                                                              >  1:2067463/1‑148 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAAc                                                                                                              >  2:3399762/1‑148 (MQ=255)
                                        tGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACa                                                                                                             >  1:1929934/1‑148 (MQ=255)
                                        tGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACa                                                                                                             >  1:58929/1‑148 (MQ=255)
                                               aaCCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGc                                                                                                      <  2:847485/148‑1 (MQ=255)
                                               aaCCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGc                                                                                                      <  2:537234/148‑1 (MQ=255)
                                                aCCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGcc                                                                                                     >  2:3064081/1‑148 (MQ=255)
                                                 ccTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCcg                                                                                                    <  2:745306/148‑1 (MQ=255)
                                                 ccTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCcg                                                                                                    <  2:2532139/148‑1 (MQ=255)
                                                   ttCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCcgcg                                                                                                  <  2:2099564/148‑1 (MQ=255)
                                                    tCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGa                                                                                                 <  1:2633069/148‑1 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               >  1:790471/1‑148 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               >  2:1572241/1‑148 (MQ=255)
                                                        tGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGcc                                                                                             <  2:3435895/148‑1 (MQ=255)
                                                         gATGCCGCTGTAGTTATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCg                                                                                            <  1:2717548/148‑1 (MQ=255)
                                                         gATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCg                                                                                            >  2:2853129/1‑148 (MQ=255)
                                                         gATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCg                                                                                            >  1:1008189/1‑148 (MQ=255)
                                                         gATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCg                                                                                            <  1:1371400/148‑1 (MQ=255)
                                                          aTGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGa                                                                                           >  2:1248893/1‑148 (MQ=255)
                                                                 tGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAg                                                                                    >  1:2080617/1‑148 (MQ=255)
                                                                    aGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCa                                                                                 >  1:3363224/1‑148 (MQ=255)
                                                                    aGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCa                                                                                 <  1:1668214/148‑1 (MQ=255)
                                                                      tGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAAt                                                                               <  2:3019791/148‑1 (MQ=255)
                                                                       gATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAAtt                                                                              <  2:3456605/148‑1 (MQ=255)
                                                                              gAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCaa                                                                       <  1:1407843/148‑1 (MQ=255)
                                                                              gAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCaa                                                                       <  1:2415680/148‑1 (MQ=255)
                                                                               aaGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAAt                                                                      >  1:411661/1‑148 (MQ=255)
                                                                               aaGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAAt                                                                      >  1:3200275/1‑148 (MQ=255)
                                                                                   gTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGc                                                                  <  1:635946/148‑1 (MQ=255)
                                                                                   gTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGc                                                                  <  1:635962/148‑1 (MQ=255)
                                                                                    tGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCa                                                                 >  2:3286614/1‑148 (MQ=255)
                                                                                       gAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCa                                                              <  2:2135207/148‑1 (MQ=255)
                                                                                         aaGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCTCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAAc                                                            <  1:551975/148‑1 (MQ=255)
                                                                                              gTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGtt                                                       >  2:1050417/1‑148 (MQ=255)
                                                                                              gTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGtt                                                       >  2:2830256/1‑148 (MQ=255)
                                                                                                  aTTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCcgc                                                   >  2:839113/1‑148 (MQ=255)
                                                                                                   tttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCcgcg                                                  >  1:3372707/1‑148 (MQ=255)
                                                                                                    ttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGt                                                 >  1:1592534/1‑148 (MQ=255)
                                                                                                    ttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGt                                                 >  1:1592537/1‑148 (MQ=255)
                                                                                                    ttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGGGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCCACCTGTTCCGCGt                                                 >  2:669008/1‑148 (MQ=255)
                                                                                                     tGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTg                                                <  2:280143/148‑1 (MQ=255)
                                                                                                      ggTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGc                                               <  1:3321686/148‑1 (MQ=255)
                                                                                                        tAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtg                                             >  2:2949659/1‑148 (MQ=255)
                                                                                                          aaGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgat                                           >  1:2550509/1‑148 (MQ=255)
                                                                                                           aGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatg                                          >  1:430517/1‑148 (MQ=255)
                                                                                                           aGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatg                                          <  2:2385726/148‑1 (MQ=255)
                                                                                                           aGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatg                                          <  2:3185328/148‑1 (MQ=255)
                                                                                                            gCCGTGCGCGAAATCCTGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatga                                         <  1:2589038/148‑1 (MQ=255)
                                                                                                            gCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatga                                         >  1:1453519/1‑148 (MQ=255)
                                                                                                            gCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatga                                         <  1:3564667/148‑1 (MQ=255)
                                                                                                            gCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatga                                         <  1:3564648/148‑1 (MQ=255)
                                                                                                             ccGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGAc                                        <  2:2041906/148‑1 (MQ=255)
                                                                                                             ccGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGAc                                        <  2:463108/148‑1 (MQ=255)
                                                                                                                tGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTa                                     <  2:2730745/148‑1 (MQ=255)
                                                                                                                 gcgcGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTAc                                    <  1:2395063/148‑1 (MQ=255)
                                                                                                                  cgcgAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTAcc                                   <  2:3278258/148‑1 (MQ=255)
                                                                                                                     gAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGa                                <  1:3220783/148‑1 (MQ=255)
                                                                                                                         tCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGcc                            >  2:532999/1‑148 (MQ=255)
                                                                                                                           aTGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCgg                          <  1:2949659/148‑1 (MQ=255)
                                                                                                                            tGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGt                         <  2:1972686/148‑1 (MQ=255)
                                                                                                                             ggCTCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTa                        >  1:1967042/1‑148 (MQ=255)
                                                                                                                             ggCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTa                        >  1:1997088/1‑148 (MQ=255)
                                                                                                                              gcgcTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTAc                       >  2:2473886/1‑148 (MQ=255)
                                                                                                                                gctggctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACATGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTg                     <  1:1878129/148‑1 (MQ=255)
                                                                                                                                  tggctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGcc                   <  1:2253280/148‑1 (MQ=255)
                                                                                                                                   ggctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCg                  <  2:2720329/148‑1 (MQ=255)
                                                                                                                                     ctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCACTGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGaa                >  1:1971236/1‑148 (MQ=255)
                                                                                                                                     ctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGaa                <  2:316186/148‑1 (MQ=255)
                                                                                                                                      tgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGaaa               <  1:532649/148‑1 (MQ=255)
                                                                                                                                        aTCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACt             >  1:3409202/1‑148 (MQ=255)
                                                                                                                                             gCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCg        <  1:1904799/148‑1 (MQ=255)
                                                                                                                                               tAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAg      <  1:673596/148‑1 (MQ=255)
                                                                                                                                               tAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAg      <  1:2460041/148‑1 (MQ=255)
                                                                                                                                                 aCCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCg    <  1:1162175/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTGACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTG  >  CP009273/2188991‑2189285

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: