Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,001,030 G→T 5.6% V603F (GTT→TTT)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,001,0300GT 5.6% 961.0 / 21.4 250V603F (GTT→TTT) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (106/130);  new base T (6/8);  total (112/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 253 total aligned reads displayed.

CGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCC  >  CP009273/4000881‑4001179
                                                                                                                                                     |                                                                                                                                                     
cgccgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATg                                                                                                                                                       >  1:1086061/1‑150 (MQ=255)
cgccgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATg                                                                                                                                                       >  1:4134488/1‑150 (MQ=255)
  ccgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGtt                                                                                                                                                     <  2:3711101/150‑1 (MQ=255)
   cgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTg                                                                                                                                                    >  1:149791/1‑150 (MQ=255)
      aCGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATg                                                                                                                                                 >  2:2683517/1‑150 (MQ=255)
        gtggtgTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGc                                                                                                                                               >  2:3952863/1‑150 (MQ=255)
           gtgtTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGAt                                                                                                                                            <  2:807536/150‑1 (MQ=255)
             gtTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGa                                                                                                                                          <  1:4110806/150‑1 (MQ=255)
              tttAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGAc                                                                                                                                         <  1:4266788/150‑1 (MQ=255)
                 aaCGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGaa                                                                                                                                      <  1:3795262/150‑1 (MQ=255)
                 aaCGACGCAACCGTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGaa                                                                                                                                      >  2:1098208/1‑150 (MQ=255)
                  angacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAag                                                                                                                                     >  2:3563045/1‑150 (MQ=255)
                  acgacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAag                                                                                                                                     >  2:1932356/1‑150 (MQ=255)
                  acgacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAag                                                                                                                                     <  2:4191400/150‑1 (MQ=255)
                  acgacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAag                                                                                                                                     <  1:3958618/150‑1 (MQ=255)
                   cgacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAaga                                                                                                                                    <  2:677968/150‑1 (MQ=255)
                     acgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGagt                                                                                                                                  >  1:625429/1‑150 (MQ=255)
                        cAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGagtagt                                                                                                                               >  2:4062263/1‑150 (MQ=255)
                          aCCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTca                                                                                                                             <  2:1160744/150‑1 (MQ=255)
                          aCCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTca                                                                                                                             <  2:3115581/150‑1 (MQ=255)
                            cTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTcagc                                                                                                                           >  1:1699962/1‑150 (MQ=255)
                            cTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTcagc                                                                                                                           >  1:4266720/1‑150 (MQ=255)
                                  tgaGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATa                                                                                                                     <  1:4489659/150‑1 (MQ=255)
                                   gagaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATaa                                                                                                                    >  1:1805586/1‑150 (MQ=255)
                                    agaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATaaa                                                                                                                   <  2:4587678/150‑1 (MQ=255)
                                      aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATaaaaa                                                                                                                 <  1:3782827/150‑1 (MQ=255)
                                        ggCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAg                                                                                                               >  1:3664260/1‑150 (MQ=255)
                                         gCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGt                                                                                                              <  1:2381180/150‑1 (MQ=255)
                                         gCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGt                                                                                                              <  2:2502919/150‑1 (MQ=255)
                                                 gATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTAt                                                                                                      <  1:3330368/150‑1 (MQ=255)
                                                  aTGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATg                                                                                                     >  2:3837481/1‑150 (MQ=255)
                                                    gCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                   >  1:709186/1‑150 (MQ=255)
                                                        aTCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               >  1:2498270/1‑150 (MQ=255)
                                                         tCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTc                                                                                              >  1:1042130/1‑150 (MQ=255)
                                                           aCCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCg                                                                                            <  1:3904078/150‑1 (MQ=255)
                                                           aCCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCg                                                                                            <  2:3655428/150‑1 (MQ=255)
                                                             cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                          <  1:1906087/150‑1 (MQ=255)
                                                                cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCTAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAAc                                                                                       <  2:2921371/150‑1 (MQ=255)
                                                                cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAAc                                                                                       <  2:3137295/150‑1 (MQ=255)
                                                                cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAAc                                                                                       <  1:2358146/150‑1 (MQ=255)
                                                                 aGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACg                                                                                      <  2:1175831/150‑1 (MQ=255)
                                                                   cGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGtt                                                                                    <  2:4362798/150‑1 (MQ=255)
                                                                     aaaTGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTc                                                                                  <  1:2817523/150‑1 (MQ=255)
                                                                        tGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTcccc                                                                               >  2:1702563/1‑150 (MQ=255)
                                                                        tGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTcccc                                                                               >  1:3034942/1‑150 (MQ=255)
                                                                         gCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTccccc                                                                              <  1:873843/150‑1 (MQ=255)
                                                                           tCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgc                                                                            >  1:3286306/1‑150 (MQ=255)
                                                                            cAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgcg                                                                           <  1:1635716/150‑1 (MQ=255)
                                                                                 ttAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                      <  1:2705075/150‑1 (MQ=255)
                                                                                   aaCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCaaa                                                                    >  1:198235/1‑150 (MQ=255)
                                                                                   aaCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCaaa                                                                    <  2:2440079/150‑1 (MQ=255)
                                                                                    aCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCCGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAat                                                                   <  2:3937033/150‑1 (MQ=255)
                                                                                    aCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAat                                                                   <  2:4536360/150‑1 (MQ=255)
                                                                                      ggCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                 <  2:603430/150‑1 (MQ=255)
                                                                                       gCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgta                                                                <  1:1208615/150‑1 (MQ=255)
                                                                                        cGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtat                                                               <  1:3565789/150‑1 (MQ=255)
                                                                                         gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                              <  2:1426298/150‑1 (MQ=255)
                                                                                         gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                              >  2:3255003/1‑150 (MQ=255)
                                                                                           tGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                            <  1:2820665/150‑1 (MQ=255)
                                                                                            gggATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTg                                                           >  1:3402405/1‑150 (MQ=255)
                                                                                               aTGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                        <  1:1523270/150‑1 (MQ=255)
                                                                                                 gcgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                      <  2:1284402/150‑1 (MQ=255)
                                                                                                 gcgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                      <  2:1284466/150‑1 (MQ=255)
                                                                                                  cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTa                                                     <  1:2620800/150‑1 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  >  1:1191737/1‑150 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  >  2:3408260/1‑150 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  <  2:2170996/150‑1 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  <  2:1041901/150‑1 (MQ=255)
                                                                                                      aGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                                 <  2:428224/150‑1 (MQ=255)
                                                                                                      aGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                                 <  2:1837460/150‑1 (MQ=255)
                                                                                                       gCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtct                                                <  1:2983994/150‑1 (MQ=255)
                                                                                                        cTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtctc                                               <  1:2740404/150‑1 (MQ=255)
                                                                                                        cTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtctc                                               >  1:3082403/1‑150 (MQ=255)
                                                                                                         tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCa                                              >  1:1564743/1‑150 (MQ=255)
                                                                                                         tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCa                                              <  1:4053295/150‑1 (MQ=255)
                                                                                                               gCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        <  2:757268/150‑1 (MQ=255)
                                                                                                               gCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        <  2:2074765/150‑1 (MQ=255)
                                                                                                               gCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        >  2:4493705/1‑150 (MQ=255)
                                                                                                               gCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        <  1:2456949/150‑1 (MQ=255)
                                                                                                                   tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCAt                                    <  1:2258919/150‑1 (MQ=255)
                                                                                                                   tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGAGGTTATTTCTCACCGTCGCCAt                                    >  2:3472468/1‑150 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTTTGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 <  2:858523/150‑1 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 <  2:2900053/150‑1 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 <  2:1232196/150‑1 (MQ=255)
                                                                                                                       aaaCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTg                                >  1:3718505/1‑150 (MQ=255)
                                                                                                                        aaCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGc                               >  1:1315186/1‑150 (MQ=255)
                                                                                                                             tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACAtt                          <  1:1387297/150‑1 (MQ=255)
                                                                                                                                aTGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtgtg                       >  2:2223037/1‑150 (MQ=255)
                                                                                                                                  ggCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgc                     <  2:2722868/150‑1 (MQ=255)
                                                                                                                                   gcgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTAAGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATATCTCACCGTCGCCATGGTGCACATTGTGgcg                    >  1:4059305/1‑150 (MQ=255)
                                                                                                                                    cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcgc                   >  2:1190858/1‑150 (MQ=255)
                                                                                                                                       tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTa                <  1:4358701/150‑1 (MQ=255)
                                                                                                                                         aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg              >  2:3792685/1‑150 (MQ=255)
                                                                                                                                         aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg              >  2:2159867/1‑150 (MQ=255)
                                                                                                                                         aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg              >  2:4504961/1‑150 (MQ=255)
                                                                                                                                            cGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGcc          >  2:4108239/1‑150 (MQ=255)
                                                                                                                                             gTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGcc          >  1:1885056/1‑150 (MQ=255)
                                                                                                                                             gTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGcc          >  1:1024491/1‑150 (MQ=255)
                                                                                                                                                  cATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt     >  1:1357362/1‑150 (MQ=255)
                                                                                                                                                     gTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTccc  <  2:2054709/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCC  >  CP009273/4000881‑4001179

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: