Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,541,504 C→T L110L (CTG→CTA)  RPPX_20280 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,541,5040CT99.6% 795.6 / NA 234L110L (CTG→CTA) RPPX_20280hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (105/128);  minor base G (0/1);  total (105/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 234 total aligned reads displayed.

CGCATGTGCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGG  >  CP009974/4541358‑4541651
                                                                                                                                                  |                                                                                                                                                   
cgcATGTGCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGt                                                                                                                                                 <  2:2252190/151‑1 (MQ=38)
 gcATGTGCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTa                                                                                                                                                <  1:1323431/151‑1 (MQ=38)
       gCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTc                                                                                                                                          >  1:4293336/1‑151 (MQ=38)
       gCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTc                                                                                                                                          <  1:3004277/151‑1 (MQ=255)
          gTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCAcg                                                                                                                                       <  2:3768603/151‑1 (MQ=38)
             tGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCAcggcg                                                                                                                                    <  2:3613684/151‑1 (MQ=255)
                gcAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCggtg                                                                                                                                 <  2:4353578/151‑1 (MQ=37)
                  aGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTTGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGTCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCggtggt                                                                                                                               <  1:2921480/151‑1 (MQ=255)
                    ggCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAg                                                                                                                             >  2:866565/1‑151 (MQ=37)
                          cGAACCGGTGATACGGGCGGTGGTGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGtcctcc                                                                                                                       <  2:4649816/151‑1 (MQ=255)
                          cGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGtcctcc                                                                                                                       <  1:5189254/151‑1 (MQ=255)
                          cGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGtcctcc                                                                                                                       <  1:4931257/151‑1 (MQ=255)
                           gAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGtcctccc                                                                                                                      >  1:3624502/1‑150 (MQ=255)
                           gAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCa                                                                                                                      <  1:5340096/151‑1 (MQ=255)
                            aaCCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGGTCTGGCCGGATAGGTCGCGCTCACGGCGGGGGTCGTCCCCCAg                                                                                                                     >  1:5106886/1‑151 (MQ=11)
                             aCCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCTGATAGGTAGCGCTCACGGCGGTGGGAGTCCTCCAGc                                                                                                                    >  2:4783036/1‑151 (MQ=255)
                             aCCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGc                                                                                                                    <  1:3168629/151‑1 (MQ=255)
                              ccGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTCGTCCTCCAGcc                                                                                                                   >  2:4714691/1‑151 (MQ=255)
                              ccGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGcc                                                                                                                   >  2:1356496/1‑151 (MQ=38)
                               cGGTGATACGGGCGGTGATGGCGTTGGTCCACAGTGGCGGGCTGTCGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGccc                                                                                                                  <  2:2810669/151‑1 (MQ=255)
                                ggTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCt                                                                                                                 >  1:1550804/1‑151 (MQ=255)
                                         ggCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTagctc                                                                                                        <  1:1009600/151‑1 (MQ=255)
                                         ggCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTagctc                                                                                                        <  1:1745644/151‑1 (MQ=255)
                                            ggTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTagctcagc                                                                                                     <  1:5139614/151‑1 (MQ=255)
                                               gATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGctct                                                                                                  >  2:2301628/1‑151 (MQ=255)
                                                gtGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCtt                                                                                                 <  2:3950398/150‑1 (MQ=255)
                                                 tGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCttt                                                                                                <  1:3256162/151‑1 (MQ=255)
                                                  ggCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTg                                                                                               <  2:3654584/151‑1 (MQ=255)
                                                  ggCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCACACACACAGCCGGTGCTCTGGCCGGCTAGGTAGCGCTCACGGCGGTGGTAGTCCTTCAGCCGTTGGTAGCTCAGCTCattt                                                                                               >  1:123677/1‑150 (MQ=14)
                                                   gCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGt                                                                                              <  1:1648484/151‑1 (MQ=255)
                                                     gTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTag                                                                                            >  2:1306565/1‑151 (MQ=255)
                                                      ttGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGTCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTaga                                                                                           >  1:5164576/1‑151 (MQ=255)
                                                        ggTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAgg                                                                                         <  1:623722/151‑1 (MQ=255)
                                                           ccACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGccc                                                                                      >  2:1913464/1‑151 (MQ=255)
                                                           ccACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGccc                                                                                      <  1:1141056/151‑1 (MQ=255)
                                                            cacaGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCcg                                                                                     >  1:1982870/1‑151 (MQ=255)
                                                              caGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCcgcg                                                                                   >  1:2986525/1‑151 (MQ=255)
                                                               aGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCgcgc                                                                                  <  1:2013810/151‑1 (MQ=255)
                                                                ggTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGcc                                                                                 >  1:4625090/1‑151 (MQ=255)
                                                                 gTGCGGGCTGTAGCCGATCCACACTTCCTTGCCCGTGCCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGTCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCg                                                                                <  2:4251245/151‑1 (MQ=255)
                                                                   gCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGAc                                                                              <  2:4961985/151‑1 (MQ=255)
                                                                     gggCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGt                                                                            <  2:3181993/151‑1 (MQ=255)
                                                                      ggCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTg                                                                           <  2:5430194/151‑1 (MQ=255)
                                                                       gCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGc                                                                          >  1:473353/1‑151 (MQ=255)
                                                                         tGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTc                                                                        <  2:5118971/151‑1 (MQ=255)
                                                                          gTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCa                                                                       >  1:352734/1‑151 (MQ=255)
                                                                            aGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatg                                                                     <  1:859883/151‑1 (MQ=255)
                                                                              ccGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaa                                                                   >  1:2333231/1‑151 (MQ=255)
                                                                              ccGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaa                                                                   <  2:2333231/151‑1 (MQ=255)
                                                                               cGATCGACACGTCCTTGGCCGTGTGCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaat                                                                  <  1:1896102/151‑1 (MQ=255)
                                                                                gATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaatg                                                                 <  1:3977869/151‑1 (MQ=255)
                                                                                gATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaatg                                                                 <  2:2955251/151‑1 (MQ=255)
                                                                                gATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCatgaatg                                                                 >  2:798871/1‑151 (MQ=255)
                                                                                  tCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCCCATGAATGAc                                                               >  1:1632605/1‑151 (MQ=255)
                                                                                   ccACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGAcc                                                              <  1:263303/151‑1 (MQ=255)
                                                                                    cacaCTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCa                                                             <  2:71146/151‑1 (MQ=255)
                                                                                    cacaCTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCa                                                             <  1:4416736/151‑1 (MQ=255)
                                                                                    cacaCTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCACGAATGACCa                                                             >  2:591732/1‑151 (MQ=255)
                                                                                        cTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCCTGAATGACCAGCtt                                                         >  2:4973375/1‑151 (MQ=255)
                                                                                        cTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCtt                                                         >  2:3914433/1‑151 (MQ=255)
                                                                                          tCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCtttt                                                       <  1:1503947/151‑1 (MQ=255)
                                                                                            cggggCCGTGTCCTGGCAGGTGTCGGCCAACACCCAGCCGGTGGTCGGGCCGGAGAGGTAGCGCTCACGGCGGGGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCg                                                     <  2:4063055/148‑1 (MQ=14)
                                                                                            cTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCg                                                     <  1:2776355/151‑1 (MQ=255)
                                                                                             ttGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGc                                                    >  2:1875040/1‑151 (MQ=255)
                                                                                                 ccGTGTCCTGGCAGGTGTCGTCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTg                                                <  2:5164576/151‑1 (MQ=255)
                                                                                                  cGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGCTAGCGCTCACGGCGGAGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGa                                               >  2:2755367/1‑151 (MQ=255)
                                                                                                  cGTGTCCTGGCAGGTGTCGGCCCACAACCAGCCGGTGCTATGGCCGGATAGGTAGCGCTCACGGTGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCACCGACGTGCTCAGGACTGACCAGCTTTTCGCCCTGa                                               >  2:3380966/1‑151 (MQ=11)
                                                                                                     gtCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAggg                                            >  1:3415444/1‑151 (MQ=18)
                                                                                                       ccTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGcc                                          >  1:2926405/1‑151 (MQ=25)
                                                                                                        cTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGa                                               >  2:4489313/1‑145 (MQ=25)
                                                                                                        cTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGa                                               <  1:4489313/145‑1 (MQ=25)
                                                                                                        cTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCt                                         <  2:1683370/151‑1 (MQ=21)
                                                                                                        cTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCt                                         <  2:4408520/151‑1 (MQ=25)
                                                                                                         tGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCtt                                        >  1:3203127/1‑151 (MQ=25)
                                                                                                              ggTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATg                                   >  1:2860881/1‑151 (MQ=25)
                                                                                                              ggTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATg                                   <  2:2477875/151‑1 (MQ=25)
                                                                                                               gtgtCGGCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGa                                  <  1:4031888/151‑1 (MQ=18)
                                                                                                                    ggCCCACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACg                             >  2:4537991/1‑151 (MQ=25)
                                                                                                                      cccACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTa                           <  1:3314500/151‑1 (MQ=25)
                                                                                                                       ccACACCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAg                          >  1:1838695/1‑151 (MQ=25)
                                                                                                                           acCCAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACAt                      >  2:1426608/1‑151 (MQ=25)
                                                                                                                            cccAGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATa                     >  1:2132551/1‑151 (MQ=21)
                                                                                                                               aGCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGcc                  <  2:20989/151‑1 (MQ=21)
                                                                                                                                gCCGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGcct                 <  2:24545/151‑2 (MQ=21)
                                                                                                                                 ccGGTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGcctg                <  1:3644923/151‑3 (MQ=21)
                                                                                                                                   ggTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTggg              <  1:3694627/151‑1 (MQ=25)
                                                                                                                                   ggTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTggg              <  2:1514485/151‑1 (MQ=25)
                                                                                                                                   ggTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTggg              <  1:1887810/151‑1 (MQ=25)
                                                                                                                                    gTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgc             <  2:3126717/151‑1 (MQ=21)
                                                                                                                                        tctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaa         <  1:3924466/151‑1 (MQ=19)
                                                                                                                                        tctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaa         <  1:995977/151‑1 (MQ=19)
                                                                                                                                        tctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaa         >  2:1132236/1‑151 (MQ=255)
                                                                                                                                         ctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaaa        >  2:2178798/1‑151 (MQ=19)
                                                                                                                                          tGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAg       <  1:5219326/151‑1 (MQ=18)
                                                                                                                                          tGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAg       <  1:3160892/151‑1 (MQ=25)
                                                                                                                                           ggCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGc      <  2:4661592/151‑1 (MQ=25)
                                                                                                                                            gCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcg     >  1:732311/1‑150 (MQ=255)
                                                                                                                                             ccGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgc    >  2:4444418/1‑149 (MQ=255)
                                                                                                                                             ccGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgc    >  2:1401336/1‑149 (MQ=255)
                                                                                                                                              cGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgcg   <  2:5244509/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCATGTGCCGTTTGCGCAGGGCTGCCGAACCGGTGATACGGGCGGTGATGGCGTTGGTCCACAGGTGCGGGCTGTAGCCGATCCACACTTCCTTGGCCGTGTCCTGGCAGGTGTCGGCCCACACCCAGCCGGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGG  >  CP009974/4541358‑4541651

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: ‑