Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,312,862 T→A 32.0% R289R (CGT→CGA argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,312,8620TA32.0% 494.3 / inf 334R289R (CGT→CGAargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (118/109);  new base A (69/38);  total (187/147)
Fisher's exact test for biased strand distribution p-value = 3.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 336 total aligned reads displayed.

GCTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGTTTGACTC  >  CP009273/3312714‑3313011
                                                                                                                                                    |                                                                                                                                                     
 cTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTgg                                                                                                                                                     >  2:1518272/1‑150 (MQ=255)
  tttGAACAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGt                                                                                                                                                    >  1:3108492/1‑150 (MQ=255)
  tttGAACAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAAACGTATCATCGAAGCGAAAAGCCGTGGt                                                                                                                                                    >  1:4680951/1‑150 (MQ=255)
   ttGAACAAGGTCATCCGGTGGCGCTGACCGGTAAACCCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTa                                                                                                                                                   <  1:2798210/150‑1 (MQ=255)
       aCAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTa                                                                                                                                               <  2:2301602/150‑1 (MQ=255)
       aCAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTa                                                                                                                                               >  2:3761538/1‑150 (MQ=255)
       aCAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTa                                                                                                                                               >  1:4610919/1‑150 (MQ=255)
         aaGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACg                                                                                                                                             <  2:2427767/150‑1 (MQ=255)
         aaGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACg                                                                                                                                             <  2:1333540/150‑1 (MQ=255)
                tCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTcc                                                                                                                                      <  2:561146/150‑1 (MQ=255)
                tCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTcc                                                                                                                                      >  2:322021/1‑150 (MQ=255)
                 ccGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCg                                                                                                                                     <  1:2377740/150‑1 (MQ=255)
                  cGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCgg                                                                                                                                    >  1:2835180/1‑150 (MQ=255)
                   ggTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCggg                                                                                                                                   >  2:1408447/1‑150 (MQ=255)
                   ggTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCggg                                                                                                                                   <  2:1305845/150‑1 (MQ=255)
                   ggTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCggg                                                                                                                                   >  1:336703/1‑150 (MQ=255)
                   ggTGGCGCTGAACGGAAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCggg                                                                                                                                   <  1:1630128/150‑1 (MQ=255)
                     tGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGa                                                                                                                                 >  2:1934496/1‑150 (MQ=255)
                       gcgcTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATg                                                                                                                               <  2:2362410/150‑1 (MQ=255)
                         gcTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGc                                                                                                                             >  2:1200331/1‑150 (MQ=255)
                          cTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCa                                                                                                                            <  2:2232591/150‑1 (MQ=255)
                            gAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCAct                                                                                                                          >  1:3930760/1‑150 (MQ=255)
                            gAACGGAAAAACCTTGAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCAct                                                                                                                          <  1:840590/150‑1 (MQ=255)
                               cGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCActgct                                                                                                                       <  2:252769/150‑1 (MQ=255)
                                ggTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCActgctg                                                                                                                      <  1:3185626/150‑1 (MQ=255)
                                 gTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACtgctgc                                                                                                                     >  2:100920/1‑150 (MQ=255)
                                  tAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGca                                                                                                                    >  2:1634679/1‑150 (MQ=255)
                                   aaaaCCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGcac                                                                                                                   >  1:604697/1‑150 (MQ=255)
                                    aaaCCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGcaca                                                                                                                  >  2:2207943/1‑150 (MQ=255)
                                     aaCCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACAt                                                                                                                 <  1:3546871/150‑1 (MQ=255)
                                       ccTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTg                                                                                                               <  1:4220677/150‑1 (MQ=255)
                                       catttAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTg                                                                                                               <  2:4048093/148‑1 (MQ=255)
                                          ttAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGt                                                                                                            <  1:2925953/150‑1 (MQ=255)
                                              cgacgaCGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGa                                                                                                        >  2:4201507/1‑150 (MQ=255)
                                               gacgacGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCATGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGaa                                                                                                       >  2:2870291/1‑150 (MQ=255)
                                                acgacgTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAAc                                                                                                      >  1:4050746/1‑150 (MQ=255)
                                                 cgacgTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACg                                                                                                     >  1:133756/1‑150 (MQ=255)
                                                     gTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTg                                                                                                 <  2:1060302/150‑1 (MQ=255)
                                                     gTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTg                                                                                                 <  2:1095778/150‑1 (MQ=255)
                                                       aGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGtt                                                                                               <  1:1172587/150‑1 (MQ=255)
                                                          aaTGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGAc                                                                                            <  2:3036468/150‑1 (MQ=255)
                                                           atgatgCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGAcc                                                                                           >  1:2871577/1‑150 (MQ=255)
                                                             gatgCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCgg                                                                                         >  2:3567167/1‑150 (MQ=255)
                                                                 cTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTAtt                                                                                     >  1:1269409/1‑150 (MQ=255)
                                                                 cTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTAtt                                                                                     >  1:3718603/1‑150 (MQ=255)
                                                                   ggAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTca                                                                                   >  1:2143296/1‑150 (MQ=255)
                                                                     aaGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTcact                                                                                 <  2:353677/150‑2 (MQ=255)
                                                                      aGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACaa                                                                                >  1:4602732/1‑150 (MQ=255)
                                                                         tAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGa                                                                             >  2:2126748/1‑150 (MQ=255)
                                                                         tAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGa                                                                             >  2:4680193/1‑150 (MQ=255)
                                                                           aCCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAg                                                                           >  2:2916433/1‑150 (MQ=255)
                                                                           aCCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAg                                                                           <  1:310463/150‑1 (MQ=255)
                                                                           aCCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAg                                                                           >  2:3874894/1‑150 (MQ=255)
                                                                           aCCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAg                                                                           >  2:4549036/1‑150 (MQ=255)
                                                                            ccGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGa                                                                          >  2:4506619/1‑150 (MQ=255)
                                                                            ccGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGa                                                                          >  2:2159483/1‑150 (MQ=255)
                                                                              gCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGaca                                                                        >  1:4141197/1‑150 (MQ=255)
                                                                               cATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGacac                                                                       >  1:4497284/1‑150 (MQ=255)
                                                                                    gcggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTg                                                                  <  1:1093566/150‑1 (MQ=255)
                                                                                     cggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGa                                                                 <  1:1875582/150‑1 (MQ=255)
                                                                                      ggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAg                                                                <  2:691383/150‑1 (MQ=255)
                                                                                      ggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAg                                                                >  2:1176420/1‑150 (MQ=255)
                                                                                      ggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAg                                                                >  1:2988565/1‑150 (MQ=255)
                                                                                      ggTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAg                                                                <  1:3446118/150‑1 (MQ=255)
                                                                                       gtcgtcACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGc                                                               <  2:4376395/150‑1 (MQ=255)
                                                                                           tcACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTa                                                           <  1:2506259/150‑1 (MQ=255)
                                                                                            cACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTAt                                                          >  1:43931/1‑150 (MQ=255)
                                                                                              cGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCa                                                        >  2:2782275/1‑150 (MQ=255)
                                                                                                gCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCAcg                                                      <  1:1511994/150‑1 (MQ=255)
                                                                                                   tGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACgcgc                                                   >  2:657709/1‑150 (MQ=255)
                                                                                                      gCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATg                                                <  2:3942150/150‑1 (MQ=255)
                                                                                                       cATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATgg                                               >  1:240083/1‑150 (MQ=255)
                                                                                                        aTGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgt                                              <  2:1540261/150‑1 (MQ=255)
                                                                                                        aTGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgt                                              <  2:336703/150‑1 (MQ=255)
                                                                                                          gAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgtcg                                            <  1:3781178/150‑1 (MQ=255)
                                                                                                          gAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgtcg                                            >  2:2016550/1‑150 (MQ=255)
                                                                                                          gAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgtcg                                            >  1:1753102/1‑150 (MQ=255)
                                                                                                           aGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGgtcgt                                           >  2:4121691/1‑150 (MQ=255)
                                                                                                               aCCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGt                                       >  1:3304544/1‑150 (MQ=255)
                                                                                                                  aGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTgg                                    <  2:1156486/150‑1 (MQ=255)
                                                                                                                    aTTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGc                                  >  1:429680/1‑150 (MQ=255)
                                                                                                                    aTTGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGc                                  >  1:2277453/1‑150 (MQ=255)
                                                                                                                     ttGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGcc                                 >  2:2602074/1‑150 (MQ=255)
                                                                                                                     ttGAAAACCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGcc                                 <  2:514644/150‑1 (MQ=255)
                                                                                                                        aaaaCCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTc                              >  1:4321060/1‑150 (MQ=255)
                                                                                                                        aaaaCCGTATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTc                              >  2:1777547/1‑150 (MQ=255)
                                                                                                                         aaaCCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTct                             >  1:3841744/1‑150 (MQ=255)
                                                                                                                          aaCCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTctg                            <  1:913156/150‑1 (MQ=255)
                                                                                                                              gTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGt                        <  2:4272267/150‑1 (MQ=255)
                                                                                                                               tATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTa                       >  1:3237058/1‑150 (MQ=255)
                                                                                                                               tATCATCGAAGCGAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTa                       <  1:3062906/150‑1 (MQ=255)
                                                                                                                                atcatcGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTAc                      >  1:1036658/1‑150 (MQ=255)
                                                                                                                                           gAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTg           >  2:4781659/1‑150 (MQ=255)
                                                                                                                                           gAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTg           <  2:1268083/150‑1 (MQ=255)
                                                                                                                                           gAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTg           >  2:777622/1‑150 (MQ=255)
                                                                                                                                           gAAAAGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTg           >  2:443975/1‑150 (MQ=255)
                                                                                                                                             aaaGCCGAGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGt         <  2:3474312/150‑1 (MQ=255)
                                                                                                                                                gCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGTTTg      <  1:1408447/150‑1 (MQ=255)
                                                                                                                                                gCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGTTTg      >  2:929720/1‑150 (MQ=255)
                                                                                                                                                  cGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGTTTGAc    >  2:2513806/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTAAAACCTTTAGCGACGACGTAGAAATGATGCTGGAAGCTAACCGCATCGGCGGTCGTCACGGCCTGGGCATGAGCGACCAGATTGAAAACCGTATCATCGAAGCGAAAAGCCGTGGTATTTACGAAGCTCCGGGGATGGCACTGCTGCACATTGCGTATGAACGCCTGTTGACCGGTATTCACAACGAAGACACCATTGAGCAGTATCACGCGCATGGTCGTCAGTTGGGCCGTCTGCTGTACCAGGGGCGTTGGTTTGACTC  >  CP009273/3312714‑3313011

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: