Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,541,648 A→G R62R (CGT→CGC RPPX_20280 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,541,6480AG99.6% 638.9 / ‑5.1 238R62R (CGT→CGCRPPX_20280hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (128/109);  total (128/110)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 243 total aligned reads displayed.

CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGCA  >  CP009974/4541500‑4541795
                                                                                                                                                    |                                                                                                                                                   
cGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgcg                                                                                                                                                   <  2:1116855/151‑1 (MQ=255)
 ggataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCgg                                                                                                                                                  <  1:2321055/147‑1 (MQ=255)
  gataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGa                                                                                                                                                 <  1:1019611/148‑1 (MQ=255)
    taGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAgg                                                                                                                                               <  1:1715658/150‑1 (MQ=255)
     aGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAggg                                                                                                                                              >  2:4509013/1‑151 (MQ=19)
      ggTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCCAAGCGCGGAGGgt                                                                                                                                             >  2:5406636/1‑151 (MQ=255)
       gTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtg                                                                                                                                            <  1:1123107/151‑1 (MQ=255)
       gTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtg                                                                                                                                            >  1:4240755/1‑151 (MQ=25)
         ggcgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTa                                                                                                                                          <  1:3503205/150‑1 (MQ=255)
         aGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTa                                                                                                                                          <  1:3601496/151‑1 (MQ=255)
           cgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAgg                                                                                                                                        >  2:6101895/1‑151 (MQ=25)
           cgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAgg                                                                                                                                        >  1:3674423/1‑151 (MQ=255)
            gcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGc                                                                                                                                       >  2:2153652/1‑151 (MQ=25)
            gcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGc                                                                                                                                       <  1:5945539/151‑1 (MQ=19)
             cTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGcc                                                                                                                                      >  2:2370174/1‑151 (MQ=255)
              tCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCg                                                                                                                                     >  2:6398736/1‑151 (MQ=25)
                 cggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGt                                                                                                                                  <  2:6595773/151‑1 (MQ=19)
                 cggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGt                                                                                                                                  <  2:820224/151‑1 (MQ=255)
                    cggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgct                                                                                                                               >  1:1432796/1‑151 (MQ=21)
                    cggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgct                                                                                                                               >  2:1225574/1‑151 (MQ=255)
                     ggtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctg                                                                                                                              <  1:5197859/151‑1 (MQ=19)
                      gtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctgc                                                                                                                             >  1:5509589/1‑151 (MQ=25)
                      gtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctgc                                                                                                                             >  2:6639328/1‑151 (MQ=255)
                         gtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGc                                                                                                                          <  2:3141835/151‑1 (MQ=16)
                         gtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTGTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGc                                                                                                                          <  2:4824675/151‑1 (MQ=255)
                           aGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGACGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCATTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGt                                                                                                                        <  1:1283757/151‑1 (MQ=255)
                              cctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTggg                                                                                                                     <  1:4448133/151‑1 (MQ=255)
                              cctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTggg                                                                                                                     <  1:5127549/151‑1 (MQ=255)
                                tccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGa                                                                                                                   <  1:441874/151‑1 (MQ=255)
                                tccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGa                                                                                                                   <  1:3783269/151‑1 (MQ=255)
                                tccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGa                                                                                                                   <  2:3902708/151‑1 (MQ=21)
                                  cAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACCCGTAGACATAGCCTGGGCCCAAACGGCGGCGGGGGTGGGCCAAGGGCGGCAGCGCGGGGAAg                                                                                                                 >  2:5163636/1‑151 (MQ=11)
                                    gCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGcc                                                                                                               >  2:1845583/1‑151 (MQ=255)
                                     cccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCt                                                                                                              >  1:4344110/1‑151 (MQ=255)
                                      ccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                             <  1:2056501/151‑1 (MQ=255)
                                       cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTg                                                                                                            <  1:581856/151‑1 (MQ=19)
                                        ttGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGt                                                                                                           >  1:1655730/1‑151 (MQ=19)
                                         tGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGtt                                                                                                          <  1:5049355/151‑1 (MQ=16)
                                          ggTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTc                                                                                                         <  1:6492942/151‑1 (MQ=16)
                                           gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                                        <  1:278590/151‑1 (MQ=17)
                                           gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                                        >  1:623750/1‑151 (MQ=16)
                                           gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGGGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCGTGTTcc                                                                                                        >  1:939829/1‑151 (MQ=255)
                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGATTGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAACGCGCGGAGGGGGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTccc                                                                                                       >  2:5358099/1‑150 (MQ=11)
                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCa                                                                                                       <  1:3908844/151‑1 (MQ=14)
                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCa                                                                                                       >  2:415744/1‑151 (MQ=9)
                                              gctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAgg                                                                                                     >  1:3326788/1‑151 (MQ=9)
                                               ctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGt                                                                                                    >  1:5019884/1‑151 (MQ=9)
                                                tcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGtt                                                                                                   >  1:5346684/1‑151 (MQ=2)
                                                            agcggCCCGCGCCGGCGTGCTCATGAATGACCAGCGTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATc                                                                                       <  2:4910340/148‑1 (MQ=0)
                                                                                           cAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTg                                                        <  1:4189875/151‑1 (MQ=9)
                                                                                            aGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGc                                                       <  1:6475456/151‑1 (MQ=9)
                                                                                             gCTTTTGGCCCTGAGGGCCTGTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCg                                                      <  2:5748783/151‑1 (MQ=9)
                                                                                                 ttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGc                                                  >  1:2219495/1‑151 (MQ=14)
                                                                                                    gCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACg                                               <  2:2667158/151‑1 (MQ=12)
                                                                                                     cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACgg                                              >  2:5499258/1‑150 (MQ=9)
                                                                                                      ccggggggCCTTTCATGTACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGAt                                             <  2:711066/146‑1 (MQ=255)
                                                                                                       cTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGGAGGCCGAGTGCTGCCGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATg                                            >  2:2721871/1‑151 (MQ=9)
                                                                                                        tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGTAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                                           <  1:3148867/151‑1 (MQ=9)
                                                                                                         gAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                                          <  2:1438510/151‑1 (MQ=12)
                                                                                                         gAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                                          <  2:2195965/151‑1 (MQ=17)
                                                                                                          aGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCg                                         >  1:499410/1‑151 (MQ=12)
                                                                                                           gggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGt                                        >  2:5975484/1‑151 (MQ=9)
                                                                                                           gggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGt                                        <  1:1773156/151‑1 (MQ=14)
                                                                                                           gggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGt                                        <  1:6992548/151‑1 (MQ=17)
                                                                                                            ggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                       >  2:3752192/1‑151 (MQ=9)
                                                                                                             gCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCCCGATGGCGTAg                                      >  2:4643250/1‑151 (MQ=9)
                                                                                                             gCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAg                                      >  2:6285607/1‑151 (MQ=14)
                                                                                                             gCCTGTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAGGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAg                                      <  1:2926638/151‑1 (MQ=9)
                                                                                                              ccTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgc                                     <  1:2425582/151‑1 (MQ=12)
                                                                                                               cTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcg                                    <  1:4817954/151‑1 (MQ=17)
                                                                                                                tttCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcgc                                   >  2:4764303/1‑151 (MQ=9)
                                                                                                                  tCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAc                                 >  1:4493123/1‑151 (MQ=14)
                                                                                                                   cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAcc                                <  1:6653680/151‑1 (MQ=16)
                                                                                                                   cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGGGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTCGCGCAcc                                >  2:5818834/1‑151 (MQ=255)
                                                                                                                    aTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCg                               >  1:4078684/1‑151 (MQ=18)
                                                                                                                    aTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCg                               <  2:542027/151‑1 (MQ=18)
                                                                                                                     tGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGTGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCgg                              <  2:5248868/151‑1 (MQ=255)
                                                                                                                      gAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGc                             >  2:2274056/1‑151 (MQ=16)
                                                                                                                       aaCACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGCGTGCTGCAGCGTGGCGAAGCATTGCACCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCa                            >  2:5109190/1‑151 (MQ=11)
                                                                                                                       aaCACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCa                            >  1:5430365/1‑151 (MQ=19)
                                                                                                                        acacGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAg                           >  1:5977627/1‑151 (MQ=255)
                                                                                                                          acGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCTGTGCGGGGCACGATGGCGTAGCGCCCCGGCAGGa                         >  1:2492713/1‑151 (MQ=255)
                                                                                                                              aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg                     <  1:717321/151‑1 (MQ=255)
                                                                                                                              aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg                     <  2:3347786/151‑1 (MQ=255)
                                                                                                                              aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg                     <  1:5323631/151‑1 (MQ=255)
                                                                                                                              aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGCGCTGCAGCGTGGGGAAGCCTGGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg                     <  1:6717209/151‑1 (MQ=255)
                                                                                                                               gACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGc                    >  1:776431/1‑151 (MQ=255)
                                                                                                                                aCATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCa                   >  1:4908232/1‑151 (MQ=255)
                                                                                                                                aCATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCa                   <  1:1347876/151‑1 (MQ=255)
                                                                                                                                  aTAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcg                 >  2:1250950/1‑151 (MQ=255)
                                                                                                                                   tAGCCTGTGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGGAGCGTGGGGAAGCCTTGCTCCCGGTTGAGGCCTGCATCGCCATAGCGGCAGGCGGGGGCGACGCCGGTGGGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc                <  2:6751855/151‑1 (MQ=11)
                                                                                                                                   tAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc                >  2:2956967/1‑151 (MQ=255)
                                                                                                                                    aGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgcg               >  2:2002393/1‑151 (MQ=255)
                                                                                                                                       gtggggGCAAAGCGCGGAGGGTGTAGGGCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgc            <  2:7082689/149‑1 (MQ=255)
                                                                                                                                       ctgggCGCAAAGCGCGGAGGGTGTAGGCCGGGTGCTGCAGCGGGGGGAAGCCTTGTTGCGGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgc            <  2:6909128/149‑1 (MQ=11)
                                                                                                                                       ctgggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgc            >  2:1320254/3‑151 (MQ=255)
                                                                                                                                           gcgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTg        >  1:1002709/1‑151 (MQ=255)
                                                                                                                                             gcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgca      <  1:6369463/151‑1 (MQ=255)
                                                                                                                                             gcAAAGCGCGGAGGGTGTAGGCCGAGAGCTGCAGCGTGGGGAAGCCTCGATCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTAACCGGTGCGCGGCACGAGGGCGTAGCGCAACGGCAGGATAGGCACGCGTGCGCTgac      >  2:4419189/1‑149 (MQ=11)
                                                                                                                                                 aGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcaggca  >  1:2223961/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGCA  >  CP009974/4541500‑4541795

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