Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,001,030 G→T 33.5% 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,0300GT33.5% 441.4 / 216.8 244V603F (GTT→TTT) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (68/94);  new base T (40/42);  total (108/136)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

CGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCC  >  CP009273/4000881‑4001179
                                                                                                                                                     |                                                                                                                                                     
cgccgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATg                                                                                                                                                       >  1:4359722/1‑150 (MQ=255)
   cgTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTg                                                                                                                                                    <  1:1719129/150‑1 (MQ=255)
    gTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGa                                                                                                                                                   <  2:4544220/150‑1 (MQ=255)
     tACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGAt                                                                                                                                                  <  2:703601/150‑1 (MQ=255)
         tggtgTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCg                                                                                                                                              >  2:3675338/1‑150 (MQ=255)
           gtgtTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGAt                                                                                                                                            <  2:1028977/150‑1 (MQ=255)
           gtgtTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGAt                                                                                                                                            >  1:1820590/1‑150 (MQ=255)
              tttAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGAc                                                                                                                                         <  2:3330684/150‑1 (MQ=255)
              tttAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGAc                                                                                                                                         <  2:1515357/150‑1 (MQ=255)
                tAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGa                                                                                                                                       <  1:4517338/150‑1 (MQ=255)
                 aaCGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGaa                                                                                                                                      >  1:1284718/1‑150 (MQ=255)
                 aaCGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGaa                                                                                                                                      <  1:3484109/150‑1 (MQ=255)
                    gacgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAagag                                                                                                                                   >  2:1701572/1‑150 (MQ=255)
                     acgCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGagt                                                                                                                                  >  1:3437479/1‑150 (MQ=255)
                        cAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGagtagt                                                                                                                               >  2:3260444/1‑150 (MQ=255)
                        cAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGagtagt                                                                                                                               <  2:1142433/150‑1 (MQ=255)
                              tgattgaGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGGTGGCGATGACGAAGAGTAGTcagcag                                                                                                                         >  2:3320333/1‑150 (MQ=255)
                               gattgaGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCagcagc                                                                                                                        <  2:2486454/150‑1 (MQ=255)
                                attgaGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAgcagca                                                                                                                       <  2:3283208/150‑1 (MQ=255)
                                  tgaGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATa                                                                                                                     >  2:3311317/1‑150 (MQ=255)
                                      aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATaaaaa                                                                                                                 <  2:1894117/150‑1 (MQ=255)
                                         gCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGt                                                                                                              <  1:881008/150‑1 (MQ=255)
                                            gAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGcc                                                                                                           >  2:774598/1‑150 (MQ=255)
                                              aCAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAg                                                                                                         <  2:1851806/150‑1 (MQ=255)
                                              aCAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAACAGTGCCAg                                                                                                         >  2:4063269/1‑150 (MQ=255)
                                              aCAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAg                                                                                                         <  2:2532366/150‑1 (MQ=255)
                                                aGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTa                                                                                                       >  2:1314990/1‑150 (MQ=255)
                                                 gATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTAt                                                                                                      <  1:2115014/150‑1 (MQ=255)
                                                     ccGATCACCGCCAGCGAAATTCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAg                                                                                                  <  2:937433/150‑1 (MQ=255)
                                                      cGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 <  1:3359083/150‑1 (MQ=255)
                                                      cGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 <  1:3357461/150‑1 (MQ=255)
                                                       gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGAc                                                                                                <  2:956528/150‑1 (MQ=255)
                                                        aTCACCGCCAGCGATATGCTCAGAGTTAACGGCGTTGGGATGCTCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               <  1:1085272/150‑1 (MQ=255)
                                                        aTCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               <  2:3059107/150‑1 (MQ=255)
                                                        aTCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               >  1:310953/1‑150 (MQ=255)
                                                         tCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTc                                                                                              >  1:2309554/1‑150 (MQ=255)
                                                          cACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTcc                                                                                             <  1:741409/150‑1 (MQ=255)
                                                             cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                          <  2:704275/150‑1 (MQ=255)
                                                              gCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTaa                                                                                         >  2:500433/1‑150 (MQ=255)
                                                                    gAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGttt                                                                                   >  2:2143625/1‑150 (MQ=255)
                                                                       aTGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTccc                                                                                >  2:2037115/1‑150 (MQ=255)
                                                                        tGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTcccc                                                                               >  2:1454097/1‑150 (MQ=255)
                                                                          cTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                             <  2:4523003/150‑1 (MQ=255)
                                                                           tCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgc                                                                            >  1:4112696/1‑150 (MQ=255)
                                                                            cAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgcg                                                                           >  1:2501840/1‑150 (MQ=255)
                                                                             aGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGa                                                                          <  2:435349/150‑1 (MQ=255)
                                                                              gCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAg                                                                         <  2:1486138/150‑1 (MQ=255)
                                                                               cGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                        <  1:3410719/150‑1 (MQ=255)
                                                                                gTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTc                                                                       >  1:245902/1‑150 (MQ=255)
                                                                                 ttAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                      <  1:1546479/150‑1 (MQ=255)
                                                                                    aCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAat                                                                   <  1:1405567/150‑1 (MQ=255)
                                                                                    aCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAat                                                                   <  2:3920429/150‑1 (MQ=255)
                                                                                     cGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatg                                                                  <  1:439352/150‑1 (MQ=255)
                                                                                     cGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatg                                                                  <  1:2241593/150‑1 (MQ=255)
                                                                                      ggCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                 >  1:461660/1‑150 (MQ=255)
                                                                                         gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                              <  1:2852316/150‑1 (MQ=255)
                                                                                           tGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                            >  1:4004845/1‑150 (MQ=255)
                                                                                                 gcgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                      <  2:3177409/150‑1 (MQ=255)
                                                                                                  cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTa                                                     >  1:941409/1‑150 (MQ=255)
                                                                                                   gcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAt                                                    <  2:889650/150‑1 (MQ=255)
                                                                                                    cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAtt                                                   <  2:2753956/150‑1 (MQ=255)
                                                                                                    cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAtt                                                   <  1:638163/150‑1 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  >  1:2295710/1‑150 (MQ=255)
                                                                                                     aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTAttt                                                  >  1:823854/1‑150 (MQ=255)
                                                                                                      aGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                                 <  2:4175615/150‑1 (MQ=255)
                                                                                                       gCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtct                                                <  1:2313767/150‑1 (MQ=255)
                                                                                                        cTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtctc                                               <  1:1118162/150‑1 (MQ=255)
                                                                                                         tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCa                                              <  1:2734743/150‑1 (MQ=255)
                                                                                                         tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCa                                              >  2:4556420/1‑150 (MQ=255)
                                                                                                          ggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCAc                                             <  1:4513336/150‑1 (MQ=255)
                                                                                                              cGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                         <  1:3235011/150‑1 (MQ=255)
                                                                                                              cGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                         >  2:1568363/1‑150 (MQ=255)
                                                                                                              cGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                         >  1:1957091/1‑150 (MQ=255)
                                                                                                                   tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCAt                                    <  1:4576503/150‑1 (MQ=255)
                                                                                                                     gCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGGTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATgg                                  <  1:3405767/150‑1 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 <  2:588163/150‑1 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 <  2:4305778/150‑1 (MQ=255)
                                                                                                                      cAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGt                                 >  2:2442770/1‑150 (MQ=255)
                                                                                                                        aaCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGc                               <  2:3733695/150‑1 (MQ=255)
                                                                                                                         aCCGTTTATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                              <  2:282973/150‑1 (MQ=255)
                                                                                                                              ttATGGCGCTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                         <  1:1765950/150‑1 (MQ=255)
                                                                                                                              ttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                         <  2:33211/150‑1 (MQ=255)
                                                                                                                                aTGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtgtg                       >  2:4147619/1‑150 (MQ=255)
                                                                                                                                 tGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTgg                      >  1:1442830/1‑150 (MQ=255)
                                                                                                                                  ggCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgc                     >  2:4466936/1‑150 (MQ=255)
                                                                                                                                   gcgcTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                    >  1:3601438/1‑150 (MQ=255)
                                                                                                                                   gcgcTGATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                    >  1:3601275/1‑150 (MQ=255)
                                                                                                                                        gATTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTAt               >  2:1056844/1‑150 (MQ=255)
                                                                                                                                         aTTCGTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg              <  1:4092424/150‑1 (MQ=255)
                                                                                                                                            cGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGc           >  1:3861923/1‑150 (MQ=255)
                                                                                                                                             gTGCGCATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGcc          >  2:2377295/1‑150 (MQ=255)
                                                                                                                                              tGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGccc         <  1:2817501/150‑1 (MQ=255)
                                                                                                                                                cgcATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCg       >  2:3801335/1‑150 (MQ=255)
                                                                                                                                                cgcATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCg       >  2:4196691/1‑150 (MQ=255)
                                                                                                                                                cgcATGTTGATGGAGATGACGAAGAGTAGTCAGCAGCATAAAAATGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTGGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCg       <  2:2998476/150‑1 (MQ=255)
                                                                                                                                                  cATTTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt     >  2:1233076/1‑150 (MQ=255)
                                                                                                                                                  cATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt     >  1:557380/1‑150 (MQ=255)
                                                                                                                                                   atgatgatgGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTTACAGTATGAAGACTCCGTAAACGATTCCCCCGCGATTCAAAAGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGTCGCTTATGAGCCCCGGta    >  2:3349099/1‑149 (MQ=255)
                                                                                                                                                   aTGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTc    <  1:1356456/150‑1 (MQ=255)
                                                                                                                                                     gTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTccc  <  2:1674422/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: