Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,597,121 C→T T614T (ACC→ACT RPPX_24915 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,597,1210CT99.8% 1407.6 / NA 482T614T (ACC→ACTRPPX_24915hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (229/252);  minor base G (0/1);  total (229/253)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

TCACCGAGCAGGCCATTACGGAAAATCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5596985‑5597268
                                                                                                                                        |                                                                                                                                                   
                            cTTGAGCCAGGAGGTGCAGCGGTATTTGAAGCGCGGCCTGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTcc                                                                                                           <  2:4268812/151‑1 (MQ=11)
                             tttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcg                                                                                                          <  2:2848385/151‑1 (MQ=14)
                              ttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCGCTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcgg                                                                                                         >  2:5163057/1‑151 (MQ=9)
                               tccccAGGCGGTGCAGCGGAATTTGAAGAGCGGCATAGCACAGGTTGAGGAGATGAAAGCGTGGCGCATCTCGGGTAGCCGCACGGTGTGGAGCGAGCTGACTACTACTGCGCGCGGGAAACCCAACTTCCTTGGGTGGCTCACCTCcggc                                                                                                        <  2:2731569/148‑1 (MQ=11)
                                aGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcg                                                                                                       <  2:4138272/151‑1 (MQ=14)
                                 gCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcgg                                                                                                      >  1:4514711/1‑151 (MQ=12)
                                   cAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtg                                                                                                    <  2:5795733/151‑1 (MQ=12)
                                    aggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtgg                                                                                                   >  2:3129358/1‑151 (MQ=12)
                                    aggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtgg                                                                                                   >  2:5350950/1‑151 (MQ=12)
                                     ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggt                                                                                                  <  1:549297/151‑1 (MQ=14)
                                       aggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtgg                                                                                                <  2:1617807/151‑1 (MQ=12)
                                        ggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtggt                                                                                               >  1:2402757/1‑151 (MQ=12)
                                        ggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtgggggt                                                                                               >  1:142904/1‑151 (MQ=12)
                                         gTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCgggggtggtg                                                                                              >  2:5213564/1‑151 (MQ=9)
                                            cAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGtt                                                                                           <  2:4078567/151‑1 (MQ=12)
                                               cGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaa                                                                                        <  1:6641141/151‑1 (MQ=12)
                                                  tATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTCTAAACg                                                                                     >  2:6104676/1‑151 (MQ=12)
                                                      tGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAgg                                                                                 <  2:5104928/151‑1 (MQ=21)
                                                       gAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAggg                                                                                >  1:3240386/1‑151 (MQ=18)
                                                          gagCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtg                                                                             >  1:5125241/1‑151 (MQ=255)
                                                           agCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgc                                                                            >  2:7051439/1‑150 (MQ=16)
                                                            gCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgct                                                                           <  1:1585163/151‑3 (MQ=255)
                                                             cGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgttg                                                                          >  1:1775427/1‑151 (MQ=255)
                                                             cGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgctg                                                                          >  1:1026213/1‑148 (MQ=12)
                                                               gCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTtgtg                                                                        <  2:661846/151‑1 (MQ=255)
                                                               gCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCtgtg                                                                        <  2:3036571/151‑1 (MQ=16)
                                                                  tGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGtt                                                                     >  1:2986162/1‑151 (MQ=255)
                                                                    gCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCa                                                                   <  2:261210/151‑1 (MQ=255)
                                                                    gCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCa                                                                   <  1:789475/151‑1 (MQ=255)
                                                                    gCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCa                                                                   >  1:2777125/1‑151 (MQ=255)
                                                                     cacaGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCaa                                                                  <  1:3824461/151‑1 (MQ=255)
                                                                     cacaGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCaa                                                                  >  2:2530479/1‑151 (MQ=16)
                                                                      acaGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAAc                                                                 >  2:3308825/1‑151 (MQ=255)
                                                                       caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACa                                                                <  1:3145447/151‑1 (MQ=255)
                                                                       caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACa                                                                <  1:3472653/151‑1 (MQ=255)
                                                                       caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACa                                                                >  2:109779/1‑151 (MQ=255)
                                                                        aGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGGTGTTGTGGTTCAAcct                                                               >  1:6370914/1‑149 (MQ=255)
                                                                         ggTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATc                                                              <  1:3188766/151‑1 (MQ=255)
                                                                         ggTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATc                                                              <  2:4629083/151‑1 (MQ=255)
                                                                         ggTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATc                                                              >  1:979402/1‑151 (MQ=255)
                                                                          gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCt                                                             <  1:4868362/151‑1 (MQ=255)
                                                                           ttGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCtt                                                            >  1:4292138/1‑151 (MQ=255)
                                                                           gtGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCtt                                                            <  2:6356631/150‑1 (MQ=255)
                                                                            tGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTa                                                           >  2:1551280/1‑151 (MQ=255)
                                                                            tGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTa                                                           >  1:7168931/1‑151 (MQ=255)
                                                                             gAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTAt                                                          >  1:214323/1‑151 (MQ=255)
                                                                               agagATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCg                                                        <  1:6857578/151‑1 (MQ=255)
                                                                                gagaTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGGTCAACATCTCATCGc                                                       >  1:4224169/1‑151 (MQ=255)
                                                                                  gaTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTc                                                     >  2:371218/1‑151 (MQ=255)
                                                                                   aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCa                                                    <  2:1906876/151‑1 (MQ=255)
                                                                                    tGAAAGCGATGCGCAACTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGAGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaa                                                   >  2:6826050/1‑151 (MQ=11)
                                                                                     gAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaac                                                  >  1:6009757/1‑150 (MQ=255)
                                                                                          cGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGc                                             <  2:3894854/151‑1 (MQ=255)
                                                                                           gTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGcc                                            >  1:5950709/1‑151 (MQ=255)
                                                                                           gTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACGACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGcc                                            <  1:2248889/151‑1 (MQ=255)
                                                                                            ttGCGCATGTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCt                                           <  2:6478999/151‑1 (MQ=255)
                                                                                             tGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTa                                          <  2:7185601/151‑1 (MQ=11)
                                                                                               cgcATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTaca                                        >  1:5682230/1‑151 (MQ=255)
                                                                                                gcATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaa                                       <  1:1810345/151‑1 (MQ=255)
                                                                                                gcATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaa                                       >  1:526738/1‑151 (MQ=255)
                                                                                                  aTCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTacaaca                                     >  2:4013138/1‑151 (MQ=255)
                                                                                                    ctcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                   >  2:4845795/1‑151 (MQ=255)
                                                                                                    ctcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGc                                   <  1:525154/151‑1 (MQ=255)
                                                                                                     tcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGcc                                  <  2:3388599/151‑1 (MQ=14)
                                                                                                     tcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGcc                                  >  1:545154/1‑151 (MQ=255)
                                                                                                      cGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCt                                 >  2:2453404/1‑151 (MQ=255)
                                                                                                       gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTg                                >  1:1193536/1‑151 (MQ=255)
                                                                                                       gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTg                                >  1:1988092/1‑151 (MQ=255)
                                                                                                        ggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                               <  2:5495366/151‑1 (MQ=255)
                                                                                                        ggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGc                               >  2:4407100/1‑151 (MQ=11)
                                                                                                         gTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGCGGTGGTGGTTCAAACGCAGGGATGTTGTGGTCCCACACCGTATCGCTCAAGACGGCCTCCAACAGCCTCCa                              >  2:5927612/1‑151 (MQ=14)
                                                                                                          tAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCag                             >  1:262652/1‑151 (MQ=255)
                                                                                                          tAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCag                             >  2:4747678/1‑151 (MQ=255)
                                                                                                          tAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCag                             >  1:4173100/1‑151 (MQ=14)
                                                                                                           aGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCaga                            >  1:2044599/1‑151 (MQ=255)
                                                                                                            gCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaa                           <  1:390962/151‑1 (MQ=255)
                                                                                                             ccGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCagaag                          <  2:600177/151‑1 (MQ=14)
                                                                                                              cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                         <  1:6952103/151‑1 (MQ=255)
                                                                                                              cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCagaaga                         >  2:379245/1‑151 (MQ=255)
                                                                                                                 gcGGTGTCTATCGCGCGGACCACTGCTGCGCGCGTGAAACCCACCGTCCTGGGGTTGGTCACCTCCGGCGGTGGTGGTGGTTTAACCGCAGGGATGTTGTGGTTCAACATCTTATCGCTGAAGCCGGCCTACAACAGCCTGCAGAagagcg                      <  2:6172861/150‑3 (MQ=14)
                                                                                                                 aCGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTg                      <  2:2027078/151‑1 (MQ=255)
                                                                                                                   ggTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGAt                    <  2:1710268/151‑1 (MQ=255)
                                                                                                                    gtgtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcgacg                   >  2:3603358/1‑146 (MQ=11)
                                                                                                                     tgtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGc                  >  2:2975615/1‑151 (MQ=255)
                                                                                                                      gtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcgacgca                 >  1:6190799/1‑144 (MQ=11)
                                                                                                                        cTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAACCGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCAcc               <  1:5827727/151‑1 (MQ=11)
                                                                                                                         tATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCg              <  1:4665835/151‑1 (MQ=255)
                                                                                                                         tATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCt                                                             <  2:2207000/104‑1 (MQ=11)
                                                                                                                         tATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCg              >  1:4760152/1‑151 (MQ=11)
                                                                                                                           tCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGa            <  1:342282/151‑1 (MQ=11)
                                                                                                                           tCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTCGCGGTGGCTTACACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCCACCACAGCCTGCAGAAGGGCGACGCACCggc            >  2:6737282/1‑150 (MQ=11)
                                                                                                                            cGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGaa           <  1:6543921/151‑1 (MQ=11)
                                                                                                                             gAGCTGGCCACTACGGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAAt          <  1:7055058/151‑1 (MQ=255)
                                                                                                                             gAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAAt          <  1:1496898/151‑1 (MQ=255)
                                                                                                                             gAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAAt          <  2:7236886/151‑1 (MQ=11)
                                                                                                                                cTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTGGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATaca       >  1:6802354/1‑151 (MQ=255)
                                                                                                                                  gACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACAcc     <  1:134950/151‑1 (MQ=255)
                                                                                                                                    ccactaCTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATACACCCt   <  1:3722923/151‑1 (MQ=11)
                                                                                                                                     cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTg  >  2:6688411/1‑151 (MQ=255)
                                                                                                                                     cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTg  <  2:1444454/151‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
TCACCGAGCAGGCCATTACGGAAAATCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5596985‑5597268

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: