Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,957,218 C→T D217D (GAC→GAT RPPX_21945 → ATPase AAA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,957,2180CT99.5% 434.4 / NA 187D217D (GAC→GATRPPX_21945ATPase AAA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (82/104);  minor base G (1/0);  total (83/104)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

ACTGGGGCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAG  >  CP009974/4957070‑4957365
                                                                                                                                                    |                                                                                                                                                   
acTGGGGCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGAtgc                                                                                                                                                   <  2:1520539/151‑1 (MQ=255)
acTGGGGCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGAtgc                                                                                                                                                   >  2:1723559/1‑151 (MQ=255)
  tGGGGCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGa                                                                                                                                                 >  1:4843927/1‑151 (MQ=14)
   ggggCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGgcgatgcgat                                                                                                                                                <  2:4165662/151‑1 (MQ=11)
     ggCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATcc                                                                                                                                              <  1:3274390/151‑1 (MQ=11)
       ccTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTc                                                                                                                                            <  1:1068378/151‑1 (MQ=11)
       ccTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTc                                                                                                                                            <  1:296711/151‑1 (MQ=11)
       ccTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTc                                                                                                                                            <  2:435635/151‑1 (MQ=255)
        cTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCg                                                                                                                                           <  1:3542785/151‑1 (MQ=255)
        cTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCg                                                                                                                                           <  2:3496920/151‑1 (MQ=255)
         tGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGAGGCGATCCGCGa                                                                                                                                          >  2:244351/1‑151 (MQ=11)
          ggCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGCTCCTCGAc                                                                                                                                         >  2:3958013/1‑151 (MQ=11)
           gCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGAcc                                                                                                                                        <  1:106139/151‑1 (MQ=11)
               cGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggtt                                                                                                                                    <  1:4015140/151‑4 (MQ=11)
               cGTTGACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggtt                                                                                                                                    <  1:1220320/151‑4 (MQ=11)
                gTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttg                                                                                                                                   >  1:4195319/1‑147 (MQ=11)
                gTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttg                                                                                                                                   >  2:2098476/1‑147 (MQ=11)
                 tttACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttgg                                                                                                                                  >  1:1260423/1‑146 (MQ=11)
                  ttACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttggt                                                                                                                                 >  1:1407733/1‑145 (MQ=11)
                  ttACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttggt                                                                                                                                 <  2:1703340/151‑7 (MQ=11)
                  ttACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCggttggt                                                                                                                                 >  2:773170/1‑145 (MQ=11)
                        cAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACaa                                                                                                                           <  1:412795/151‑1 (MQ=11)
                        cAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACaa                                                                                                                           >  2:2505272/1‑151 (MQ=11)
                         aGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGCGCACAAc                                                                                                                          >  2:2246199/1‑151 (MQ=11)
                            gcAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCt                                                                                                                       >  1:3976263/1‑151 (MQ=11)
                            gcAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCGTGGGCGATGCGATCCTCGACCGGGTGGTGCACAACGcc                                                                                                                       >  2:284707/1‑150 (MQ=11)
                             cAACGCCGCGACATGCTCGAACTGCTGGGCGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCtt                                                                                                                      <  2:3659868/151‑1 (MQ=11)
                              aaCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTa                                                                                                                     >  1:645968/1‑151 (MQ=11)
                              aaCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTa                                                                                                                     >  2:5103230/1‑151 (MQ=11)
                               aCGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTAt                                                                                                                    <  1:1178023/151‑1 (MQ=11)
                                cgccgcGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATc                                                                                                                   <  2:2856062/151‑1 (MQ=11)
                                  ccgcGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCgg                                                                                                                 <  1:3885777/151‑1 (MQ=11)
                                  ccgcGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCgg                                                                                                                 <  1:4533239/151‑1 (MQ=11)
                                                                cTACGGCGACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAg                                                                                   <  2:3206121/151‑1 (MQ=11)
                                                                 tACGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGc                                                                                  <  1:3001951/151‑1 (MQ=11)
                                                                  aCGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCg                                                                                 <  2:242747/151‑1 (MQ=11)
                                                                  aCGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCg                                                                                 >  1:906375/1‑151 (MQ=11)
                                                                   cGGCAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGt                                                                                <  2:4054014/151‑1 (MQ=11)
                                                                      cAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcg                                                                             <  2:2491419/151‑1 (MQ=11)
                                                                      cAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcg                                                                             <  1:2057087/151‑1 (MQ=11)
                                                                      cAACCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcg                                                                             <  2:3771218/151‑1 (MQ=11)
                                                                       aaCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcga                                                                            >  2:4601276/1‑151 (MQ=11)
                                                                        aCCGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgac                                                                           <  1:786759/151‑1 (MQ=11)
                                                                         ccGCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgacg                                                                          >  1:2001745/1‑151 (MQ=11)
                                                                         ccGCTCGACACTGGGGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGcgacg                                                                          <  2:897132/151‑1 (MQ=11)
                                                                           gCTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaa                                                                        <  2:4837517/151‑1 (MQ=14)
                                                                            cTCGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGaaa                                                                       >  1:1777780/1‑151 (MQ=11)
                                                                              cGACACTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAAtt                                                                     >  2:36597/1‑151 (MQ=11)
                                                                                 cacTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAc                                                                  <  2:4727819/151‑1 (MQ=14)
                                                                                 cacTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAc                                                                  <  1:5103230/151‑1 (MQ=14)
                                                                                 cacTGGTGACCAGCCAGAGGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAc                                                                  <  2:3310936/151‑1 (MQ=14)
                                                                                  acTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAcg                                                                 >  1:5231724/1‑151 (MQ=11)
                                                                                  acTGGTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAcg                                                                 <  1:4274282/151‑1 (MQ=14)
                                                                                     ggTGACCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGAcggcg                                                              >  2:1202712/1‑151 (MQ=14)
                                                                                       ggaCCAGCAAGATGCCGGTGGATAAATGGCATGCACTGCTCGGAGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACggcggc                                                            <  2:2281687/149‑1 (MQ=11)
                                                                                         aCCAGCCAGATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAg                                                          <  2:1003663/151‑1 (MQ=14)
                                                                                           cagccagATGCCGGTGGATACATGGCATGCGCTGCTCGGGGATCCGACCTTGGGCCATGCGACCCTCGACGGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAggg                                                        <  1:4398045/151‑1 (MQ=14)
                                                                                           cagccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAggg                                                        >  2:2569019/1‑151 (MQ=14)
                                                                                            agccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGa                                                       >  1:1961941/1‑151 (MQ=14)
                                                                                             gccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAc                                                      >  2:3849259/1‑151 (MQ=14)
                                                                                             gccagATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGAc                                                      <  1:4601276/151‑1 (MQ=14)
                                                                                                 gATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCa                                                  >  1:4295893/1‑151 (MQ=14)
                                                                                                 gATGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCa                                                  >  2:4353299/1‑151 (MQ=14)
                                                                                                  aTGCCGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAg                                                 >  1:2502838/1‑151 (MQ=14)
                                                                                                     ccGGTGGATAAATGGCATGCACGGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACt                                              <  2:2985898/151‑1 (MQ=14)
                                                                                                      cGGTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTa                                             >  1:3732107/1‑151 (MQ=14)
                                                                                                       ggTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTaa                                            >  1:4762994/1‑151 (MQ=14)
                                                                                                        gTGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAAc                                           <  2:5213341/151‑1 (MQ=14)
                                                                                                         tGGATAAATGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACa                                          <  2:984893/151‑1 (MQ=14)
                                                                                                            ataaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATg                                       >  1:1829733/3‑151 (MQ=14)
                                                                                                            ataaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATg                                       <  1:3638698/149‑1 (MQ=14)
                                                                                                             taaaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGc                                      <  1:1342597/150‑1 (MQ=14)
                                                                                                              aaaGGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCa                                     <  2:1879455/151‑1 (MQ=14)
                                                                                                               aaTGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCaa                                    >  1:4222470/1‑151 (MQ=14)
                                                                                                                 tGGCATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAAAGCTTATCGGATCGAACTGAGGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGCCTTCAGACTAACAATGCAAcc                                  <  2:1777504/151‑1 (MQ=14)
                                                                                                                    cATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTg                               <  1:2611512/151‑1 (MQ=255)
                                                                                                                    cATGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTg                               >  1:1981068/1‑151 (MQ=14)
                                                                                                                     aTGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGc                              <  2:645968/151‑1 (MQ=255)
                                                                                                                     aTGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGc                              >  1:2174842/1‑151 (MQ=255)
                                                                                                                      tGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCg                             <  2:3423356/151‑1 (MQ=255)
                                                                                                                      tGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCg                             <  2:806730/151‑1 (MQ=255)
                                                                                                                      tGCACTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCg                             >  2:3025120/1‑151 (MQ=255)
                                                                                                                          cTGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGc                         <  2:1077214/151‑1 (MQ=255)
                                                                                                                           tGATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCt                        <  1:3659727/151‑1 (MQ=255)
                                                                                                                            gATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCaa                                                                                    <  2:135452/90‑1 (MQ=11)
                                                                                                                            gATCGGCGATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTg                       >  1:2215371/1‑151 (MQ=255)
                                                                                                                                   gATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCg                <  1:4000810/151‑1 (MQ=255)
                                                                                                                                   gATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCg                >  2:1235631/1‑151 (MQ=255)
                                                                                                                                   gATCCGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCg                >  1:3160973/1‑151 (MQ=255)
                                                                                                                                      ccGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGTCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCACACTAACAATCTAACCCTGCGCCCCTGTGCTctcact             >  2:3902615/1‑146 (MQ=14)
                                                                                                                                      ccGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACt             <  1:4665204/151‑1 (MQ=255)
                                                                                                                                       cGACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTg            >  2:2759229/1‑151 (MQ=255)
                                                                                                                                        gACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGc           >  1:3175091/1‑151 (MQ=255)
                                                                                                                                        gACCTTGGGCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGc           <  1:4027637/151‑1 (MQ=255)
                                                                                                                                              gggCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTcc     <  1:4017787/151‑1 (MQ=255)
                                                                                                                                               ggCGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCa    >  1:1688821/1‑151 (MQ=255)
                                                                                                                                                gcgatgcgatCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCaa   <  1:1961526/151‑1 (MQ=255)
                                                                                                                                                gcgatgcgatCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCaa   >  1:4566445/1‑151 (MQ=255)
                                                                                                                                                 cGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAg  >  1:2868377/1‑151 (MQ=255)
                                                                                                                                                 cGATGCGATCCTCGACCGGTTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAg  <  2:2634875/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTGGGGCCTGGCACCGTTTACCGCAGCGCAACGCCGCGACATGCTCGAACTGCTGGACGACCGCTACGGCAACCGCTCGACACTGGTCACCAGCCAGATGCCGGTGGACAAATGGCATGCGCTGATCGGCGATCCGACCTTGGGCGACGCGATCCTCGACCGTCTGGTGCACAACGCTTATCGGATCGAACTGAAGGGCGAATCGATGCGCAAGCGTGCGACGAAATTGACGGCGGCAGGGACTTCAGACTAACAATGCAACCCTGCGTCGCTGCGCTCCGACTGCCTGTCCAAG  >  CP009974/4957070‑4957365

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