Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,999,692 G→T 12.8% D157Y (GAT→TAT)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,999,6920GT12.8% 802.2 / 63.6 234D157Y (GAT→TAT) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (96/108);  new base T (16/14);  total (112/122)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 236 total aligned reads displayed.

TGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGC  >  CP009273/3999543‑3999841
                                                                                                                                                     |                                                                                                                                                     
tgatgaCAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACg                                                                                                                                                       >  2:2567064/1‑150 (MQ=255)
  atgaCAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGAt                                                                                                                                                     >  2:1050775/1‑150 (MQ=255)
  atgaCAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGAt                                                                                                                                                     <  2:2350644/150‑1 (MQ=255)
    gaCAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACTAttt                                                                                                                                                   >  1:3693123/1‑150 (MQ=255)
     aCAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTc                                                                                                                                                  <  1:3203325/150‑1 (MQ=255)
       aGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGAGTTCCg                                                                                                                                                >  1:34150/1‑150 (MQ=255)
          cTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGccc                                                                                                                                             <  1:1522458/150‑1 (MQ=255)
           tGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCg                                                                                                                                            >  2:2126901/1‑150 (MQ=255)
           tGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCg                                                                                                                                            >  2:306994/1‑150 (MQ=255)
             ccGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGa                                                                                                                                          <  2:4319651/150‑1 (MQ=255)
              cGCACTGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGaa                                                                                                                                         >  2:3374383/1‑150 (MQ=255)
                  ccGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATg                                                                                                                                     <  1:385211/150‑1 (MQ=255)
                   cGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGc                                                                                                                                    <  1:2281691/150‑1 (MQ=255)
                      gCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCcgc                                                                                                                                 <  1:3061899/150‑1 (MQ=255)
                       cAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCcgcg                                                                                                                                >  1:4179643/1‑150 (MQ=255)
                        aaaTTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCgcgc                                                                                                                               >  2:581009/1‑150 (MQ=255)
                         aaTTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCt                                                                                                                              >  2:704380/1‑150 (MQ=255)
                          aTTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTc                                                                                                                             >  2:4613506/1‑150 (MQ=255)
                          aTTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTc                                                                                                                             <  2:4390192/150‑1 (MQ=255)
                           ttCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCg                                                                                                                            <  2:2633213/150‑1 (MQ=255)
                            tCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCgg                                                                                                                           <  2:1519956/150‑1 (MQ=255)
                              gTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGTATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCAACGATTTCCGCCCGGAATATGCCGCGCTCGGTc                                                                                                                         >  2:2236659/1‑150 (MQ=255)
                                cTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAg                                                                                                                       >  1:2697237/1‑150 (MQ=255)
                                cTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAg                                                                                                                       <  1:1440220/150‑1 (MQ=255)
                                  gCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGtt                                                                                                                     >  2:3743090/1‑150 (MQ=255)
                                    tttATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTgc                                                                                                                   >  2:1424414/1‑150 (MQ=255)
                                     ttATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTgcg                                                                                                                  <  2:2673497/150‑1 (MQ=255)
                                     ttATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTgcg                                                                                                                  <  1:1578523/150‑1 (MQ=255)
                                       atatCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCACTGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGcc                                                                                                                <  2:3631564/150‑1 (MQ=255)
                                            gCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCgg                                                                                                           >  1:1014/1‑150 (MQ=255)
                                              cccGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGtt                                                                                                         >  2:1315946/1‑150 (MQ=255)
                                              cccGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGtt                                                                                                         >  2:1868568/1‑150 (MQ=255)
                                              cccGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGtt                                                                                                         >  1:2450824/1‑150 (MQ=255)
                                               ccGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTc                                                                                                        <  2:2487097/150‑1 (MQ=255)
                                               ccGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTc                                                                                                        <  1:3833286/150‑1 (MQ=255)
                                                cGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTcc                                                                                                       >  2:3841598/1‑150 (MQ=255)
                                                 ggAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTccc                                                                                                      <  2:3314917/150‑1 (MQ=255)
                                                  gAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCg                                                                                                     >  2:2963229/1‑150 (MQ=255)
                                                    aCGCCTGTTTCTGGATAACTTTCTTGAGCATCTGGCGCCCTGGAATCCGGTGTTATTAGCCGTTGATGAAGCTCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGTATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGAc                                                                                                   <  1:4415043/150‑1 (MQ=255)
                                                    aCGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGAc                                                                                                   <  2:3483133/150‑1 (MQ=255)
                                                    aCGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGAc                                                                                                   <  2:1165690/150‑1 (MQ=255)
                                                      gCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGc                                                                                                 >  2:468783/1‑150 (MQ=255)
                                                       ccTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCt                                                                                                <  2:1689122/150‑1 (MQ=255)
                                                        cTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTg                                                                                               <  2:4543025/150‑1 (MQ=255)
                                                           aTGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCg                                                                                            >  1:4218445/1‑150 (MQ=255)
                                                            tGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGt                                                                                           >  2:4468280/1‑150 (MQ=255)
                                                                ggATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTAt                                                                                       >  2:3697377/1‑150 (MQ=255)
                                                                 gATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATg                                                                                      >  1:2346437/1‑150 (MQ=255)
                                                                  aTAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATgg                                                                                     <  2:3672504/150‑1 (MQ=255)
                                                                   tAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGgc                                                                                    >  2:2290226/1‑150 (MQ=255)
                                                                    aaCTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAATCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGgcg                                                                                   >  1:2768829/1‑150 (MQ=255)
                                                                    aaCTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGgcg                                                                                   <  2:395990/150‑1 (MQ=255)
                                                                     aCTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGgcgc                                                                                  >  1:2571980/1‑150 (MQ=255)
                                                                     aCTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGgcgc                                                                                  >  2:157893/1‑150 (MQ=255)
                                                                          cTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGAcc                                                                             <  2:2484097/150‑1 (MQ=255)
                                                                             gAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGcc                                                                          >  1:2514379/1‑150 (MQ=255)
                                                                             gAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGcc                                                                          >  2:1816745/1‑150 (MQ=255)
                                                                                cATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCcaca                                                                       <  1:4313515/150‑1 (MQ=255)
                                                                                cATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCcaca                                                                       <  1:570658/150‑1 (MQ=255)
                                                                                   cTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGcc                                                                    >  1:4597338/1‑150 (MQ=255)
                                                                                    tGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCcg                                                                   >  2:1627049/1‑150 (MQ=255)
                                                                                     ggCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCcga                                                                  >  1:1836734/1‑150 (MQ=255)
                                                                                       cgcACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCcgacg                                                                <  2:3801304/150‑1 (MQ=255)
                                                                                         cACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCgacgac                                                              <  1:746613/150‑1 (MQ=255)
                                                                                            tGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGAcacc                                                           >  2:839576/1‑150 (MQ=255)
                                                                                               aaTCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCAcg                                                        <  2:420337/150‑1 (MQ=255)
                                                                                                    ggTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCa                                                   <  2:67823/150‑1 (MQ=255)
                                                                                                     gtgtTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTNACCGCCACAGCCGACGACACCACGCGCCAg                                                  <  1:4281912/150‑1 (MQ=255)
                                                                                                       gtTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGa                                                >  2:2825649/1‑150 (MQ=255)
                                                                                                       gtTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGa                                                >  2:2836796/1‑150 (MQ=255)
                                                                                                         tattaGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGata                                              >  2:1537804/1‑150 (MQ=255)
                                                                                                          attaGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGatat                                             <  2:240623/150‑1 (MQ=255)
                                                                                                           ttaGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATc                                            <  2:523146/150‑1 (MQ=255)
                                                                                                             aGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGt                                          >  2:2858487/1‑150 (MQ=255)
                                                                                                                 gTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTgcgc                                      >  1:2988168/1‑150 (MQ=255)
                                                                                                                    gatgaAGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCctg                                   <  1:366682/150‑1 (MQ=255)
                                                                                                                    gatgaAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCctg                                   <  1:1299556/150‑1 (MQ=255)
                                                                                                                     atgaAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCctgc                                  >  1:4074370/1‑150 (MQ=255)
                                                                                                                        aaGCGCACTGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTgg                               <  1:4389169/150‑1 (MQ=255)
                                                                                                                        aaGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTgg                               >  2:3169988/1‑150 (MQ=255)
                                                                                                                            gcACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCt                           <  1:4470033/150‑1 (MQ=255)
                                                                                                                            gcACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCt                           <  1:3557935/150‑1 (MQ=255)
                                                                                                                             cACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTg                          <  2:1505875/150‑1 (MQ=255)
                                                                                                                              aCTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGa                         <  1:4150474/150‑1 (MQ=255)
                                                                                                                                tGTATCTCCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAAc                       <  2:1081694/150‑1 (MQ=255)
                                                                                                                                tGTATCTCCCAATGGGGCCACTATGTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAAc                       <  2:1081631/150‑1 (MQ=255)
                                                                                                                                 gTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACg                      <  2:4226644/150‑1 (MQ=255)
                                                                                                                                  tATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGa                     >  2:2035612/1‑150 (MQ=255)
                                                                                                                                     ctcCCAATGGGGCCACTATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATcc                  <  1:2030533/150‑1 (MQ=255)
                                                                                                                                          aaTGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGa             >  1:3745815/1‑150 (MQ=255)
                                                                                                                                           aTGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGAt            >  2:1894204/1‑150 (MQ=255)
                                                                                                                                           aTGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGAt            <  2:2211485/150‑1 (MQ=255)
                                                                                                                                            tGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGAtt           <  1:3241457/150‑1 (MQ=255)
                                                                                                                                             ggggCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCATGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTc          <  2:3574846/150‑1 (MQ=255)
                                                                                                                                                gCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCaaa       <  2:902928/150‑1 (MQ=255)
                                                                                                                                                  cACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATc     >  1:1141495/1‑150 (MQ=255)
                                                                                                                                                    cGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATcag   <  2:3473426/150‑1 (MQ=255)
                                                                                                                                                    cGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATcag   >  2:1203487/1‑150 (MQ=255)
                                                                                                                                                    cGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATcag   >  1:4355934/1‑150 (MQ=255)
                                                                                                                                                     gATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATcagc  <  2:415539/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATGACAGGCTGCCGCACCGGGCAAATTCGTCTGCTTTATATCGCCCCGGAACGCCTGATGCTGGATAACTTTCTTGAGCATCTGGCGCACTGGAATCCGGTGTTATTAGCCGTTGATGAAGCGCACTGTATCTCCCAATGGGGCCACGATTTCCGCCCGGAATATGCCGCGCTCGGTCAGTTGCGCCAGCGGTTCCCGACGCTGCCGTTTATGGCGCTGACCGCCACAGCCGACGACACCACGCGCCAGGATATCGTGCGCCTGCTGGGGCTGAACGATCCGCTGATTCAAATCAGC  >  CP009273/3999543‑3999841

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: