Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,756,936 T→C D390G (GAT→GGT)  RPPX_16780 ← phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,756,9360TC99.3% 993.3 / ‑2.9 279D390G (GAT→GGT) RPPX_16780phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (144/133);  total (145/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 279 total aligned reads displayed.

AATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCACCTTGCGC  >  CP009974/3756788‑3757082
                                                                                                                                                    |                                                                                                                                                  
aaTCAGGAATGCCGACGCCAGCGACTGGGCGGCATTCATCCGCGGATCACGATGGGTGTGATACCGGTCAGACAACGGATCCTCGGGGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCg                                                                                                                                                  <  2:5828640/151‑1 (MQ=255)
aaTCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCg                                                                                                                                                  <  1:6174921/151‑1 (MQ=255)
aaTCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCg                                                                                                                                                  <  1:2110508/151‑1 (MQ=255)
         tGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATgg                                                                                                                                         <  2:6306257/151‑1 (MQ=255)
          gCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGc                                                                                                                                        >  2:544354/1‑151 (MQ=255)
          gCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGc                                                                                                                                        >  2:5352076/1‑151 (MQ=255)
           ccAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCt                                                                                                                                       <  2:1595185/151‑1 (MQ=255)
            cAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTg                                                                                                                                      >  1:4537268/1‑151 (MQ=255)
            cAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTg                                                                                                                                      >  1:4573978/1‑151 (MQ=255)
             aaCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGc                                                                                                                                     <  1:5063527/151‑1 (MQ=255)
              aCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcc                                                                                                                                    >  2:1418198/1‑151 (MQ=255)
              aCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcc                                                                                                                                    >  1:723681/1‑151 (MQ=255)
               cTCGAGCGACTGATCGGGATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGccc                                                                                                                                   <  2:5907966/151‑1 (MQ=255)
                cccAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCt                                                                                                                                  <  2:885535/150‑1 (MQ=255)
                      gACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCt                                                                                                                            >  1:4745612/1‑151 (MQ=255)
                        cTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTgg                                                                                                                          >  1:5256315/1‑151 (MQ=255)
                          gATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtg                                                                                                                        <  2:2260180/151‑1 (MQ=255)
                          gATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtg                                                                                                                        >  2:5214579/1‑151 (MQ=255)
                           aTCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                                                       >  2:2575273/1‑151 (MQ=255)
                           aTCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                                                       >  2:1545379/1‑151 (MQ=255)
                               gCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACtt                                                                                                                   <  2:1012689/151‑1 (MQ=255)
                                 aTTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCCTGATATGGATACCGCCGGCATGGCTGCCCTCGGCCCGGTGCACTTgg                                                                                                                 >  2:4670875/1‑151 (MQ=255)
                                   tCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaa                                                                                                               >  2:3549902/1‑151 (MQ=255)
                                    cATCCGCGTATCATCATTGGTTTGATACCGGTCAGACAACGCATCCGCGGTGATCGGCCGGCCACCGCCGATGCATGCGGTGACGTTCCGCCCAGTCACCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaag                                                                                                              <  2:729038/151‑1 (MQ=255)
                                    cATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaag                                                                                                              <  2:4903876/151‑1 (MQ=255)
                                     aTCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaaga                                                                                                             >  1:701062/1‑151 (MQ=255)
                                      tCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaagaa                                                                                                            <  2:1127781/151‑1 (MQ=255)
                                        cgcgGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACt                                                                                                          <  2:5588880/151‑1 (MQ=255)
                                         gcgGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTg                                                                                                         >  1:1887500/1‑151 (MQ=255)
                                           ggATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACGTGGAAGACCTGCt                                                                                                       >  2:6128266/1‑151 (MQ=255)
                                             aTCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTc                                                                                                     >  1:3752115/1‑151 (MQ=255)
                                              tCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCa                                                                                                    >  2:3156607/1‑151 (MQ=255)
                                                acaATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGTGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCAcc                                                                                                  <  1:179012/151‑1 (MQ=255)
                                                 caATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCt                                                                                                 >  2:1782414/1‑151 (MQ=255)
                                                 caATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAACAACTGCTTCAccc                                                                                                 >  2:3508558/1‑150 (MQ=255)
                                                   aTGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCa                                                                                               <  2:4745612/151‑1 (MQ=255)
                                                     gggTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATc                                                                                             >  1:1242482/1‑151 (MQ=255)
                                                       gtgtGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCa                                                                                           <  2:343716/151‑1 (MQ=255)
                                                        ggtgATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAg                                                                                          <  2:5122429/150‑1 (MQ=255)
                                                             tACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCt                                                                                     >  2:3774466/1‑151 (MQ=255)
                                                             tACCGGTCAGACAACGCATCCTCGGCGATCGGCCGGGCACCGTCGATACATCCCGTGACGTTCTGCCCAGTCATCTCGATATAGATACCGCCGACATGGCTGTCCTCCGGCTCGAGTACTTGGTGGAACTGCTTTACATATTCCAGAATCt                                                                                     >  2:4077215/1‑151 (MQ=255)
                                                               ccGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgc                                                                                   >  1:1031965/1‑151 (MQ=255)
                                                                 ggTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgcgc                                                                                 <  1:2054652/151‑1 (MQ=255)
                                                                 ggTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgcgc                                                                                 <  1:4062206/151‑1 (MQ=255)
                                                                  gTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCa                                                                                >  2:3327836/1‑151 (MQ=255)
                                                                  gTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCa                                                                                >  1:5927586/1‑151 (MQ=255)
                                                                   tCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCaa                                                                               >  2:1921302/1‑151 (MQ=255)
                                                                     aGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCaaaa                                                                             >  1:3633709/1‑151 (MQ=255)
                                                                      gACAACGCATCCTGGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAAt                                                                            <  2:5699640/151‑1 (MQ=255)
                                                                       aCAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATc                                                                           <  1:2711537/151‑1 (MQ=255)
                                                                             cATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGt                                                                     >  2:2037621/1‑151 (MQ=255)
                                                                                ccTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCtt                                                                  >  2:6330579/1‑151 (MQ=255)
                                                                                ccTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCtt                                                                  >  1:3093838/1‑151 (MQ=255)
                                                                                  tCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGt                                                                >  1:1466647/1‑151 (MQ=255)
                                                                                    gggggtCGGCAGGGCACCGCAGATGCATTCGGTGACGTTCTGCGCAGTGATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATGGCGGGTCTGGTaa                                                              <  2:6037145/146‑1 (MQ=255)
                                                                                     gTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAAc                                                             >  1:6006807/1‑151 (MQ=255)
                                                                                            gacGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCgctgct                                                      <  1:6343151/149‑1 (MQ=255)
                                                                                            gCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCgctgct                                                      <  2:2100932/151‑1 (MQ=255)
                                                                                                ggCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGcc                                                  <  2:1241533/151‑1 (MQ=255)
                                                                                                  cACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCtt                                                >  1:3888277/1‑151 (MQ=255)
                                                                                                   aCCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTg                                               >  1:4233374/1‑151 (MQ=255)
                                                                                                    ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                              >  1:736997/1‑151 (MQ=255)
                                                                                                    ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                              >  1:2103366/1‑151 (MQ=255)
                                                                                                    ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                              <  1:6286800/151‑1 (MQ=255)
                                                                                                    ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                              <  2:642220/151‑1 (MQ=255)
                                                                                                     cgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGAt                                             >  1:2415909/1‑151 (MQ=255)
                                                                                                       ccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTtggtgg                                           >  2:2839422/1‑151 (MQ=255)
                                                                                                        cgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgt                                          >  1:1751625/1‑151 (MQ=255)
                                                                                                         gATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgtg                                         <  2:6241046/151‑1 (MQ=255)
                                                                                                         gATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgtg                                         <  2:4623989/151‑1 (MQ=255)
                                                                                                           tGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGtt                                       <  2:3752514/151‑1 (MQ=255)
                                                                                                               ttCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCg                                   <  2:249092/151‑1 (MQ=255)
                                                                                                                tCGGTCACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGt                                  <  2:4472872/151‑1 (MQ=255)
                                                                                                                 cGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTg                                 <  2:3633709/151‑1 (MQ=255)
                                                                                                                 cGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTg                                 >  2:4937244/1‑151 (MQ=255)
                                                                                                                  ggTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGc                                >  1:2647043/1‑151 (MQ=255)
                                                                                                                   gTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCa                               >  1:5356601/1‑151 (MQ=255)
                                                                                                                       cgTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatccg                           >  1:4973944/2‑149 (MQ=255)
                                                                                                                       cgTTCTGCCCAGTCATCTCGAGATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcgg                           <  2:5082193/150‑1 (MQ=255)
                                                                                                                        gTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcgga                          <  2:2725956/151‑1 (MQ=255)
                                                                                                                        gTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcgga                          >  1:1957112/1‑151 (MQ=255)
                                                                                                                         ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                         >  1:6436030/1‑151 (MQ=255)
                                                                                                                         ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                         >  1:3787886/1‑151 (MQ=255)
                                                                                                                         ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                         >  1:5531494/1‑151 (MQ=255)
                                                                                                                         ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                         >  1:1659723/1‑151 (MQ=255)
                                                                                                                            tGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggatcgg                      >  2:43167/1‑151 (MQ=255)
                                                                                                                            tGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggatcgg                      >  2:129860/1‑151 (MQ=255)
                                                                                                                             gCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGCGCAccggatcgga                     >  1:582051/1‑151 (MQ=255)
                                                                                                                              cccagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGaa                    <  1:2152980/147‑1 (MQ=255)
                                                                                                                              cccAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGaa                    >  1:2286498/1‑151 (MQ=255)
                                                                                                                                cagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACt                  <  2:2757763/149‑1 (MQ=255)
                                                                                                                                cagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACt                  <  2:1217330/149‑1 (MQ=255)
                                                                                                                                cagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACt                  <  2:3293963/149‑1 (MQ=255)
                                                                                                                                 agTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTc                 <  2:5299190/150‑1 (MQ=255)
                                                                                                                                  gTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTcc                >  2:2330080/1‑151 (MQ=255)
                                                                                                                                    cagctcgATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTACCCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCcac              <  1:3263906/145‑1 (MQ=255)
                                                                                                                                    cATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCcac              <  1:5422551/151‑1 (MQ=255)
                                                                                                                                       ctcgATATGGATACCGCCGGCGTGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGc           <  2:1836216/148‑1 (MQ=255)
                                                                                                                                              tggATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCACCTTgc    <  2:4537268/150‑1 (MQ=255)
                                                                                                                                               ggATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCACCTTgcg   <  1:4395423/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCACCTTGCGC  >  CP009974/3756788‑3757082

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: