Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,138 C→A 11.2% 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,1380CA11.2% 541.5 / 30.4 160R454S (CGC→AGC) metGmethionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (61/81);  new base A (7/11);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 163 total aligned reads displayed.

TTTGACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTG  >  CP009273/2188991‑2189285
                                                                                                                                                   |                                                                                                                                                   
   gACGGCGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCt                                                                                                                                                  <  1:1220269/148‑1 (MQ=255)
        cGTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTc                                                                                                                                             >  1:3093876/1‑148 (MQ=255)
         gTGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCg                                                                                                                                            <  1:319849/148‑1 (MQ=255)
          tGCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGa                                                                                                                                           <  1:1745519/148‑1 (MQ=255)
           gCTGGCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGAt                                                                                                                                          >  2:2300079/1‑148 (MQ=255)
              ggCAAGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGaa                                                                                                                                       >  1:854192/1‑148 (MQ=255)
                 aaGCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAg                                                                                                                                    <  1:391282/148‑1 (MQ=255)
                   gCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGc                                                                                                                                  >  2:496001/1‑148 (MQ=255)
                   gCGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGc                                                                                                                                  <  1:2005134/148‑1 (MQ=255)
                    cGAACTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCt                                                                                                                                 <  2:1579567/148‑1 (MQ=255)
                       aCTGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTTAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCg                                                                                                                              <  1:1781682/148‑1 (MQ=255)
                         tGGCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtg                                                                                                                            >  1:1891740/1‑148 (MQ=255)
                          ggCTGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtgg                                                                                                                           <  1:1025244/148‑1 (MQ=255)
                             tGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtgggtg                                                                                                                        <  1:2020568/148‑1 (MQ=255)
                             tGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtgggtg                                                                                                                        <  1:2050250/148‑1 (MQ=255)
                             tGACCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCgtgggtg                                                                                                                        <  2:2427134/148‑1 (MQ=255)
                               aCCCGCAGTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGggtggt                                                                                                                      >  1:2682961/1‑148 (MQ=255)
                                      gTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGaaa                                                                                                               <  2:2195870/148‑1 (MQ=255)
                                      gTTGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGaaa                                                                                                               >  2:1289793/1‑148 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGACTTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAAc                                                                                                              >  2:1483063/1‑148 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAAc                                                                                                              <  1:2976337/148‑1 (MQ=255)
                                       ttGTACAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAAc                                                                                                              >  2:536940/1‑148 (MQ=255)
                                           aCAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGa                                                                                                          >  1:2070125/1‑148 (MQ=255)
                                           aCAAAACCTTCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGa                                                                                                          >  1:168241/1‑148 (MQ=255)
                                                   ttCACTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCcgcg                                                                                                  <  1:2199283/148‑1 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               <  2:2055373/148‑1 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               <  2:1401436/148‑1 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               <  2:3129151/148‑1 (MQ=255)
                                                      aCTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATg                                                                                               <  2:3129440/148‑1 (MQ=255)
                                                       cTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGc                                                                                              <  2:341761/148‑1 (MQ=255)
                                                       cTGATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGc                                                                                              >  2:3293230/1‑148 (MQ=255)
                                                         gATGCCGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCg                                                                                            <  1:679177/148‑1 (MQ=255)
                                                              cGCTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTg                                                                                       >  2:1427697/1‑148 (MQ=255)
                                                               gCTGAAGTGATTGGTGAAGCGTGGGGAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGc                                                                                      <  2:1311804/148‑1 (MQ=255)
                                                                cTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCa                                                                                     <  1:1743429/148‑1 (MQ=255)
                                                                cTGAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCa                                                                                     <  1:356197/148‑1 (MQ=255)
                                                                  gAAGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAgg                                                                                   <  2:1828623/148‑1 (MQ=255)
                                                                    aGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCa                                                                                 <  2:1122481/148‑1 (MQ=255)
                                                                    aGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCa                                                                                 <  2:1122475/148‑1 (MQ=255)
                                                                    aGTGATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCa                                                                                 <  2:3227085/148‑1 (MQ=255)
                                                                       gATTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAAtt                                                                              <  2:2308944/148‑1 (MQ=255)
                                                                        aTTGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAAttt                                                                             >  2:668603/1‑148 (MQ=255)
                                                                         ttGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTg                                                                            >  1:508228/1‑148 (MQ=255)
                                                                         ttGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTg                                                                            >  1:508276/1‑148 (MQ=255)
                                                                         ttGGTGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTg                                                                            <  1:1427697/148‑1 (MQ=255)
                                                                             tGAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCa                                                                        <  2:3187388/148‑1 (MQ=255)
                                                                              gAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCaa                                                                       <  2:1652423/148‑1 (MQ=255)
                                                                              gAAGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCaa                                                                       >  1:2421625/1‑148 (MQ=255)
                                                                               aaGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAAt                                                                      <  1:365052/148‑1 (MQ=255)
                                                                               aaGCGTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAAt                                                                      <  1:1199161/148‑1 (MQ=255)
                                                                                   gTGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGc                                                                  <  1:476335/148‑1 (MQ=255)
                                                                                    tGGGAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCa                                                                 >  1:1152140/1‑148 (MQ=255)
                                                                                     gggAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCAt                                                                <  2:3265561/148‑1 (MQ=255)
                                                                                      ggAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATc                                                               <  1:66683/148‑1 (MQ=255)
                                                                                      ggAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATc                                                               <  1:66656/148‑1 (MQ=255)
                                                                                       gAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCa                                                              >  1:2334093/1‑148 (MQ=255)
                                                                                       gAAAGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCa                                                              >  1:2333929/1‑148 (MQ=255)
                                                                                        aaaGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCaa                                                             <  2:1878945/148‑1 (MQ=255)
                                                                                        aaaGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCaa                                                             <  2:166367/148‑1 (MQ=255)
                                                                                         aaGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAAc                                                            <  1:1709008/148‑1 (MQ=255)
                                                                                         aaGCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAAc                                                            <  1:32582/148‑1 (MQ=255)
                                                                                           gCCGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCt                                                          >  1:2867196/1‑148 (MQ=255)
                                                                                             cGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCTATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGAATCAACATGt                                                        >  1:1905253/1‑148 (MQ=255)
                                                                                             cGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGt                                                        <  2:2179579/148‑1 (MQ=255)
                                                                                             cGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGt                                                        <  1:34304/148‑1 (MQ=255)
                                                                                             cGTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGt                                                        <  2:3029019/148‑1 (MQ=255)
                                                                                              gTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGtt                                                       <  2:1089885/148‑1 (MQ=255)
                                                                                              gTGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGtt                                                       <  1:2967971/148‑1 (MQ=255)
                                                                                               tGAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTc                                                      >  1:199213/1‑148 (MQ=255)
                                                                                                gAATTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTcc                                                     >  1:2863300/1‑148 (MQ=255)
                                                                                                  aTTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCcgc                                                   >  2:2853292/1‑148 (MQ=255)
                                                                                                  aTTTGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCcgc                                                   >  2:2853154/1‑148 (MQ=255)
                                                                                                   tttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCcgcg                                                  <  1:1518097/148‑1 (MQ=255)
                                                                                                    ttGGTAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGt                                                 >  2:2066926/1‑148 (MQ=255)
                                                                                                        tAAAGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtg                                             <  2:3246980/148‑1 (MQ=255)
                                                                                                         aaaGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtga                                            <  1:1818486/148‑1 (MQ=255)
                                                                                                           aGCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCTCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatg                                          >  1:1619230/1‑148 (MQ=255)
                                                                                                            gCCGTGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCtgatga                                         >  2:984255/1‑148 (MQ=255)
                                                                                                                tGCGCGAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTa                                     <  1:2341884/148‑1 (MQ=255)
                                                                                                                 gcgcGAAATCATGGCGCTGGCTGATCTGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTAc                                    <  1:3252072/148‑1 (MQ=255)
                                                                                                                   gcgAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCt                                  <  1:275123/148‑1 (MQ=255)
                                                                                                                     gAAATCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGa                                <  2:1892075/148‑1 (MQ=255)
                                                                                                                       aaTCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAg                              >  1:1118709/1‑148 (MQ=255)
                                                                                                                        aTCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTGCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGc                             <  2:2134762/148‑1 (MQ=255)
                                                                                                                        aTCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGc                             <  2:1715405/148‑1 (MQ=255)
                                                                                                                         tCATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGcc                            >  2:201577/1‑148 (MQ=255)
                                                                                                                          cATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCg                           >  2:125404/1‑148 (MQ=255)
                                                                                                                          cATGGCGCTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCg                           >  2:1157831/1‑148 (MQ=255)
                                                                                                                               cgcTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTCCt                      >  2:84430/1‑148 (MQ=255)
                                                                                                                               cgcTGGCTGATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACt                      >  2:213213/1‑148 (MQ=255)
                                                                                                                                   ggctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCg                  >  2:2219030/1‑148 (MQ=255)
                                                                                                                                    gctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGa                 >  2:2600344/1‑148 (MQ=255)
                                                                                                                                    gctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGa                 >  2:2600317/1‑148 (MQ=255)
                                                                                                                                     ctgATCTGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGaa                <  1:842662/148‑1 (MQ=255)
                                                                                                                                           tGGCTAACCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGAc          <  2:77042/148‑1 (MQ=255)
                                                                                                                                           tGGCTAACAGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGAc          <  1:1981881/148‑1 (MQ=255)
                                                                                                                                             gCTAACAGCTATGTCGATGAGCAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCg        >  1:1020582/1‑148 (MQ=255)
                                                                                                                                                aaCCGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGc     >  1:3111196/1‑148 (MQ=255)
                                                                                                                                                 aCCGATATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCg    <  1:582716/148‑1 (MQ=255)
                                                                                                                                                   cGCTATGTCGATGAACAGGCTCCGTGGGTGGTGGCGAAACAGGAAGGCCGCGATGCCGACCTGCAGGCAATTTGCTCAATGGGCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTg  >  1:2414730/1‑148 (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: