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% 1107.6 / NA 393T614T (ACC→ACTRPPX_24915hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (182/210);  total (182/211)
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 400 total aligned reads displayed.

CTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5597013‑5597268
                                                                                                            |                                                                                                                                                   
cgttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTcc                                                                                                           <  1:4888475/149‑1 (MQ=2)
cTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTcc                                                                                                           >  1:3254891/1‑151 (MQ=2)
cTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTcc                                                                                                           <  2:2894169/151‑1 (MQ=2)
 tttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcg                                                                                                          <  1:1494976/151‑1 (MQ=9)
 tttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcg                                                                                                          <  1:5426714/151‑1 (MQ=9)
 tttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcg                                                                                                          <  1:2290250/151‑1 (MQ=9)
 tttAGCCAGCAGGTGCAGCGGTAGTTGAAGAGCGGCAGGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcg                                                                                                          <  2:6116613/151‑1 (MQ=11)
  ttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcgg                                                                                                         <  1:3348594/151‑1 (MQ=14)
  ttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcgg                                                                                                         <  1:2031503/151‑1 (MQ=9)
   tAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggc                                                                                                        >  1:585603/1‑151 (MQ=14)
   tAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggc                                                                                                        >  1:1587003/1‑151 (MQ=14)
   tAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggc                                                                                                        >  1:2825940/1‑151 (MQ=12)
    aGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGCCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcg                                                                                                       <  2:4071495/151‑1 (MQ=9)
     gCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcgg                                                                                                      >  1:6356889/1‑151 (MQ=17)
      ccAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggc                                                                                                        <  2:1345205/148‑1 (MQ=17)
      ccAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggc                                                                                                        >  1:1345205/1‑148 (MQ=17)
      ccAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggt                                                                                                     >  1:6463636/1‑151 (MQ=12)
       cAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtg                                                                                                    >  1:2634807/1‑151 (MQ=12)
       cAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtg                                                                                                    >  1:3528239/1‑151 (MQ=12)
         ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggt                                                                                                  <  2:3046716/151‑1 (MQ=12)
         ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggt                                                                                                  >  2:4146867/1‑151 (MQ=12)
         ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggt                                                                                                  >  1:4538739/1‑151 (MQ=12)
         ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGCTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggg                                                                                                  >  1:1942157/1‑150 (MQ=9)
               gCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtggtggt                                                                                            <  1:861168/151‑1 (MQ=12)
                 aGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGttt                                                                                          <  2:1986459/151‑1 (MQ=12)
                   cGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaa                                                                                        <  2:3241129/151‑1 (MQ=12)
                   cGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaa                                                                                        <  2:1588514/151‑1 (MQ=12)
                    ggTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                       >  1:4945334/1‑151 (MQ=12)
                      tATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGGGGGGGTTTAAACg                                                                                     >  2:1165393/1‑151 (MQ=9)
                           gAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAggg                                                                                >  2:5487941/1‑151 (MQ=255)
                            aaGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGa                                                                               >  2:5216733/1‑151 (MQ=19)
                              gagCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtg                                                                             >  1:3788953/1‑151 (MQ=255)
                                gCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgtt                                                                           >  1:6233460/1‑151 (MQ=255)
                                gCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgct                                                                           >  2:623163/1‑149 (MQ=255)
                                 cGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGTGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTGGGGTTGCTCACCTCCGGCGGTGGTGGGGGTTTAAACGCAGGGAtgttg                                                                          >  2:3297564/1‑151 (MQ=11)
                                  ggCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgctgt                                                                         >  2:4529032/1‑151 (MQ=9)
                                    cATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTCCTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTgg                                                                       <  1:2987483/151‑1 (MQ=255)
                                    cATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTgg                                                                       >  2:2305653/1‑151 (MQ=9)
                                    cATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTgg                                                                       >  2:6035027/1‑151 (MQ=255)
                                      tGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGtt                                                                     >  2:1715254/1‑151 (MQ=255)
                                      tGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGtt                                                                     <  2:394225/151‑1 (MQ=255)
                                       ggCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTc                                                                    >  1:606248/1‑151 (MQ=255)
                                         cacaGGTTGAAGAGATGAAAGCGTTGCCCATCTCGGCTAGCCGCACGGAGGCTATCGAGATGACCTCTACTGCGCGCGTAAAACCCAACAACGTAGGACAGCTCACCTCCCACCGTGGTGGTGGAGTAAACGCAGAGATGCTATGGTTCaa                                                                  >  2:6468442/1‑151 (MQ=2)
                                           caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACa                                                                >  2:1476166/1‑151 (MQ=255)
                                           caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACa                                                                >  1:6152252/1‑151 (MQ=255)
                                           caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACa                                                                <  2:5084536/151‑1 (MQ=12)
                                            aGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACAt                                                               >  1:970754/1‑151 (MQ=255)
                                            aGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACAt                                                               <  2:4871457/151‑1 (MQ=255)
                                              gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCt                                                             <  1:1049992/151‑1 (MQ=255)
                                                tGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTa                                                           >  2:4828042/1‑151 (MQ=255)
                                                 gaacggaTGAAAGGGTGGCGGATCTCGGGTAGCCGCACGGTGTCTATCGGGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCGCCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTAt                                                          <  2:5151583/146‑1 (MQ=14)
                                                 gAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTAt                                                          <  2:4028687/151‑1 (MQ=255)
                                                 gAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTAt                                                          >  2:2482345/1‑151 (MQ=255)
                                                  gagagATGAAAGCGTTGCGGATGTCGGGTAGCCGCACGGTGGCTATCGAGCTGACCACTACTTCGCGCGAGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATc                                                         <  1:2392533/150‑1 (MQ=14)
                                                   agagATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCg                                                        >  2:3577073/1‑151 (MQ=255)
                                                    gagaTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGc                                                       <  2:310202/151‑1 (MQ=255)
                                                    gagaTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGc                                                       <  2:1294267/151‑1 (MQ=255)
                                                       aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCa                                                    >  1:484206/1‑151 (MQ=255)
                                                       aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTGCAACATCTTATCGCTCa                                                    >  2:5366507/1‑151 (MQ=255)
                                                        tGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaa                                                   <  1:660478/151‑1 (MQ=255)
                                                         gAAAGCGTTTCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaac                                                  >  1:3897827/1‑150 (MQ=11)
                                                         gAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaac                                                  >  2:3002969/1‑150 (MQ=255)
                                                          aaaGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaaca                                                 <  2:1863621/151‑3 (MQ=255)
                                                           aaGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGAc                                                <  2:585603/151‑1 (MQ=255)
                                                             gCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACgg                                              >  1:4303045/1‑151 (MQ=255)
                                                              cGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTCGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGc                                             <  2:2158741/151‑1 (MQ=255)
                                                                 tGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTa                                          <  1:2120271/151‑1 (MQ=255)
                                                                   cgcATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTaca                                        <  1:6328509/151‑1 (MQ=11)
                                                                    gcATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaa                                       >  2:5288605/1‑151 (MQ=255)
                                                                       tctcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAg                                    >  1:6015999/1‑151 (MQ=255)
                                                                        ctcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                   <  1:1700942/151‑1 (MQ=255)
                                                                         tcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGcc                                  <  1:4003218/151‑1 (MQ=255)
                                                                           gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCGCCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTg                                <  1:4437826/151‑1 (MQ=11)
                                                                           gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTg                                <  1:3166608/151‑1 (MQ=255)
                                                                           gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTg                                <  2:2712890/151‑1 (MQ=255)
                                                                           gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTg                                >  2:1097729/1‑151 (MQ=14)
                                                                           gggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTg                                <  1:5503202/151‑1 (MQ=14)
                                                                            ggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGc                               >  1:2439423/1‑151 (MQ=255)
                                                                             gTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCa                              <  1:4097020/151‑1 (MQ=255)
                                                                              tAGCCGCGCGGTGTCTATCGAGCTGACCGCTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCag                             <  2:3833855/151‑1 (MQ=14)
                                                                                gCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaa                           >  2:747705/1‑151 (MQ=255)
                                                                                  cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                         >  1:4604710/1‑151 (MQ=255)
                                                                                  cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCagaaga                         >  1:3049210/1‑151 (MQ=255)
                                                                                     aCGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcg                      >  2:4217673/1‑149 (MQ=11)
                                                                                       ggTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGac                    <  1:1318374/151‑2 (MQ=11)
                                                                                          gtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCa                 <  2:4768245/151‑1 (MQ=255)
                                                                                            cTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCAcc               <  1:5045550/151‑1 (MQ=255)
                                                                                             tATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCTGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACTGCCTACAACAGCCTGCAGAAGAGCGACGCACCg              >  1:2878505/1‑151 (MQ=11)
                                                                                              aTCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCgg             <  2:6406467/151‑1 (MQ=255)
                                                                                                 gAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAAt          <  1:1236536/151‑1 (MQ=255)
                                                                                                 gAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAAt          >  1:4219934/1‑151 (MQ=255)
                                                                                                  aGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATa         <  1:5284993/151‑1 (MQ=255)
                                                                                                  aGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATa         <  2:5450293/151‑1 (MQ=11)
                                                                                                   gCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATac        <  1:1866995/151‑1 (MQ=255)
                                                                                                    cTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATaca       >  1:5662270/1‑151 (MQ=11)
                                                                                                      gACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCa                                                    >  2:2604090/1‑104 (MQ=14)
                                                                                                      gACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGGGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGGGAGGCACCGGAATACAcc     >  2:818033/1‑151 (MQ=14)
                                                                                                       accactACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATACAccc    <  2:5721306/151‑1 (MQ=11)
                                                                                                        ccactaCTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCt   <  1:5391510/151‑1 (MQ=255)
                                                                                                        ccactaCTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCt   >  1:216721/1‑151 (MQ=255)
                                                                                                            |                                                                                                                                                   
CTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5597013‑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: