Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,579,422 C→A 5.9% E254D (GAG→GAT ggt ← gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,579,4220CA 5.9% 1298.6 / 42.0 336E254D (GAG→GATggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (143/173);  new base A (6/14);  total (149/187)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 350 total aligned reads displayed.

GGATTTGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTAC  >  CP009273/3579274‑3579570
                                                                                                                                                    |                                                                                                                                                    
ggATTTGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCAtc                                                                                                                                                      <  2:952502/149‑1 (MQ=255)
ggATTTGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCAtc                                                                                                                                                      <  2:446659/149‑1 (MQ=255)
 gATTTGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCAtct                                                                                                                                                     <  1:2822833/149‑1 (MQ=255)
   tttGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTcc                                                                                                                                                   <  1:5074565/149‑1 (MQ=255)
   tttGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTcc                                                                                                                                                   >  2:2271619/1‑149 (MQ=255)
   tttGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTcc                                                                                                                                                   >  2:3344683/1‑149 (MQ=255)
   tttGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATATcc                                                                                                                                                   <  2:2104734/149‑1 (MQ=255)
    ttGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCt                                                                                                                                                  <  2:4618836/149‑1 (MQ=255)
     tGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTg                                                                                                                                                 >  1:4370931/1‑149 (MQ=255)
      gTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTgg                                                                                                                                                <  2:2921503/149‑1 (MQ=255)
      gTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTgg                                                                                                                                                <  2:1203727/149‑1 (MQ=255)
      gTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTgg                                                                                                                                                >  1:2151999/1‑149 (MQ=255)
      gTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTgg                                                                                                                                                <  1:1165382/149‑1 (MQ=255)
          gATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGa                                                                                                                                            >  2:4369457/1‑149 (MQ=255)
           atatGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGAt                                                                                                                                           <  2:4033101/149‑1 (MQ=255)
              tGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGATTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTATTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTg                                                                                                                                        >  1:74620/1‑149 (MQ=255)
              tGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTg                                                                                                                                        <  1:3975965/149‑1 (MQ=255)
                 aTCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTc                                                                                                                                     >  2:30067/1‑149 (MQ=255)
                 aTCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTc                                                                                                                                     >  1:4133388/1‑149 (MQ=255)
                  tCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTcc                                                                                                                                    >  2:2578608/1‑149 (MQ=255)
                      gccgGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCaa                                                                                                                                <  1:157941/149‑1 (MQ=255)
                         ggaggaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCg                                                                                                                             >  1:4802444/1‑149 (MQ=255)
                          gaggaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGt                                                                                                                            <  1:5021831/149‑1 (MQ=255)
                            ggaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGc                                                                                                                          >  2:1035156/1‑149 (MQ=255)
                            ggaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGc                                                                                                                          <  1:2705489/149‑1 (MQ=255)
                             gaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGcc                                                                                                                         <  2:3644585/149‑1 (MQ=255)
                             gaTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATATCCTGGGCGATCTGTTCCGCAATCGTGcc                                                                                                                         <  1:574510/149‑1 (MQ=255)
                              aTGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCt                                                                                                                        >  1:4684389/1‑149 (MQ=255)
                               tGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCtt                                                                                                                       <  1:3765666/149‑1 (MQ=255)
                                 gCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTta                                                                                                                     <  1:3461722/149‑1 (MQ=255)
                                 gCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTta                                                                                                                     <  1:3668547/149‑1 (MQ=255)
                                   ggTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTtata                                                                                                                   <  1:319164/149‑1 (MQ=255)
                                   ggTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTtata                                                                                                                   <  2:3109352/149‑1 (MQ=255)
                                        cATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGACtt                                                                                                              >  2:4707330/1‑149 (MQ=255)
                                        cATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAAtt                                                                                                              <  1:1356088/149‑1 (MQ=255)
                                        cATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAAtt                                                                                                              >  2:4953581/1‑149 (MQ=255)
                                         aTGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATATCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTc                                                                                                             <  1:781705/149‑1 (MQ=255)
                                          tGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCg                                                                                                            >  1:1255156/1‑149 (MQ=255)
                                            gAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTc                                                                                                          <  1:1796764/149‑1 (MQ=255)
                                                  aCCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGcc                                                                                                    >  1:129032/1‑149 (MQ=255)
                                                   ccTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCg                                                                                                   >  1:3819580/1‑149 (MQ=255)
                                                    cTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGt                                                                                                  <  2:712850/149‑1 (MQ=255)
                                                    cTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGt                                                                                                  >  1:3090030/1‑149 (MQ=255)
                                                         aCCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCa                                                                                             >  2:1752710/1‑149 (MQ=255)
                                                         aCCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCa                                                                                             >  1:2488755/1‑149 (MQ=255)
                                                                 tAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAatcat                                                                                     <  1:96039/149‑1 (MQ=255)
                                                                  aaTCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAatcatc                                                                                    <  1:2054448/149‑1 (MQ=255)
                                                                    tCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTTCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTAAGCAatcatcac                                                                                  >  1:1677309/1‑147 (MQ=255)
                                                                     cgccgcTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTcc                                                                                 <  2:4523655/149‑1 (MQ=255)
                                                                      gccgcTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCa                                                                                <  2:560221/149‑1 (MQ=255)
                                                                       ccgcTTATCGGAGTGCGTTCGACCGCTTGATAGGCTGCTAAATCGTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAg                                                                               <  2:3258794/149‑1 (MQ=255)
                                                                           ttATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGctc                                                                           <  2:2871964/149‑1 (MQ=255)
                                                                           ttATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGctc                                                                           >  2:3809757/1‑149 (MQ=255)
                                                                           ttATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGctc                                                                           >  1:3723897/1‑149 (MQ=255)
                                                                           ttATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATATCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGctc                                                                           <  1:2151527/149‑1 (MQ=255)
                                                                            tATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGctct                                                                          >  1:670502/1‑149 (MQ=255)
                                                                             aTCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCtt                                                                         >  2:5116466/1‑149 (MQ=255)
                                                                              tCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAACTCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCttt                                                                        <  2:3842367/149‑1 (MQ=255)
                                                                                   gTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAg                                                                   <  1:1017557/149‑1 (MQ=255)
                                                                                     gCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCCCCAGGCTCTTTGCCAGGt                                                                 <  2:1431224/149‑1 (MQ=255)
                                                                                      cGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGtt                                                                <  1:2310768/149‑1 (MQ=255)
                                                                                       gTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTc                                                               <  1:4191745/149‑1 (MQ=255)
                                                                                           gACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTAAGCAATCATATCCAGGCTATTTGCCAGGTTCGCCt                                                           >  1:3542181/1‑149 (MQ=255)
                                                                                               gCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCAc                                                       >  2:4255383/1‑149 (MQ=255)
                                                                                               gCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCAc                                                       >  2:3714801/1‑149 (MQ=255)
                                                                                                         ctgctAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCg                                             >  1:2954928/1‑149 (MQ=255)
                                                                                                           gctAAATCTTCTGTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGATCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGcc                                           >  2:1542533/1‑149 (MQ=255)
                                                                                                            ctAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGccc                                          >  2:4694827/1‑149 (MQ=255)
                                                                                                             tAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCt                                         <  2:2060979/149‑1 (MQ=255)
                                                                                                             tAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCAAATCCAGGCTCTTTGCCAGGTTCGCCAGCACCAGCGTGTCGCCCt                                         >  1:1865424/1‑149 (MQ=255)
                                                                                                              aaaTCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCtt                                        <  2:2811129/149‑1 (MQ=255)
                                                                                                              aaaTCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCtt                                        >  2:1824571/1‑149 (MQ=255)
                                                                                                               aaTCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCACCGTGTCGCCCttt                                       >  1:2764006/1‑149 (MQ=255)
                                                                                                                 tcttctTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCAACTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCttttt                                     >  1:1361946/1‑149 (MQ=255)
                                                                                                                  cttcttTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTc                                    <  2:2376267/149‑1 (MQ=255)
                                                                                                                  cttcttTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTc                                    >  2:131465/1‑149 (MQ=255)
                                                                                                                     cttTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGc                                 >  1:2122205/1‑149 (MQ=255)
                                                                                                                       ttAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCgg                               >  2:508933/1‑149 (MQ=255)
                                                                                                                           tGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCg                           >  2:3037395/1‑149 (MQ=255)
                                                                                                                             aTCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGcc                         >  2:1654557/1‑149 (MQ=255)
                                                                                                                              tCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGccc                        >  1:4263683/1‑149 (MQ=255)
                                                                                                                              tCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGccc                        <  1:2649951/149‑1 (MQ=255)
                                                                                                                                 aGCCACCGTTTTTCTGCATATCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCctct                     <  2:5148478/149‑1 (MQ=255)
                                                                                                                                  gCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCtt                    >  2:4004946/1‑149 (MQ=255)
                                                                                                                                   ccACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCttt                   <  2:2968780/149‑1 (MQ=255)
                                                                                                                                   ccACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCttt                   <  1:4476568/149‑1 (MQ=255)
                                                                                                                                    cACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTc                  <  1:1157260/149‑1 (MQ=255)
                                                                                                                                     aCCGTTTTTCTGCATATCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTcc                 <  2:649165/149‑1 (MQ=255)
                                                                                                                                      ccGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCCGGTTCGCCTGCACCAGCGTGTCGCCCATTTTCAGCGGCTCGCCCTCTTTCCa                >  2:1379084/1‑149 (MQ=255)
                                                                                                                                      ccGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCa                <  1:5117495/149‑1 (MQ=255)
                                                                                                                                      ccGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCa                <  1:2571422/149‑1 (MQ=255)
                                                                                                                                         tttttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCagaa             <  2:2641887/149‑1 (MQ=255)
                                                                                                                                         tttttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCagaa             <  2:3538444/149‑1 (MQ=255)
                                                                                                                                            ttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGAt          <  2:4840177/149‑1 (MQ=255)
                                                                                                                                             tCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATa         <  1:1067924/149‑1 (MQ=255)
                                                                                                                                                gCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCt      <  2:1242484/149‑1 (MQ=255)
                                                                                                                                                 cATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCtt     <  2:1949599/149‑1 (MQ=255)
                                                                                                                                                 cATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATAAAAACCACGCTCTTTGCCAGGGTCGCCTGCACCACCGTGTCGCCCATTTTCAGCGGCGCGCCCGCTTTGCAGAAGATatcgt     >  1:2229660/1‑145 (MQ=255)
                                                                                                                                                  aTCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCttt    >  2:1145775/1‑149 (MQ=255)
                                                                                                                                                    ctcCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTAc  <  1:292594/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATTTGTACGATATGGATCCCGCCGGAGGATGGCGGTGGCATGGAGTAAACCTGATACCCGCGATAATCGCCGCTTATCGGAGTGCGTTCGACCGCTTTATAGGCTGCTAAATCTTCTTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTAC  >  CP009273/3579274‑3579570

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