Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,813,962 T→C 7.8% Q190R (CAA→CGA)  ligB ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,813,9620TC 7.8% 1096.2 / 61.1 331Q190R (CAA→CGA) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (148/157);  new base C (20/6);  total (168/163)
Fisher's exact test for biased strand distribution p-value = 7.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 350 total aligned reads displayed.

CCAGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTTG  >  CP009273/3813813‑3814111
                                                                                                                                                     |                                                                                                                                                      
 cAGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATttgttgttg                                                                                                                                                       <  2:4019146/150‑1 (MQ=255)
 cAGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATttgttgttg                                                                                                                                                       >  2:4531854/1‑150 (MQ=255)
     tctTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGata                                                                                                                                                   >  1:4003124/1‑150 (MQ=255)
      ctTTTAAACGATCAGACATTAACTGCGGTCCATCCTGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGatat                                                                                                                                                  <  2:1534990/150‑1 (MQ=255)
          tAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGccc                                                                                                                                              <  1:4899540/150‑1 (MQ=255)
           aaaCGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGcccc                                                                                                                                             <  2:4222226/150‑1 (MQ=255)
            aaCGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCAGTTATTCCCCCCATTGGTTGTTGGATATGCCCCt                                                                                                                                            >  2:4252468/1‑150 (MQ=255)
              cGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCAGTTGTTGTTGGATATGCCCCTcg                                                                                                                                          >  1:4766398/1‑150 (MQ=255)
                aacAGACATGACCTGCGGTCCATCCGGCCATGCCCAGACAAACACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTcgcg                                                                                                                                        <  1:2718024/148‑1 (MQ=255)
                   aGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTg                                                                                                                                     >  2:4967890/1‑150 (MQ=255)
                   aGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTg                                                                                                                                     >  1:4692723/1‑150 (MQ=255)
                   aGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTg                                                                                                                                     >  2:4351488/1‑150 (MQ=255)
                    gACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTgg                                                                                                                                    <  2:1061064/150‑1 (MQ=255)
                     aCATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGga                                                                                                                                   <  2:1644911/150‑1 (MQ=255)
                      cATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGgag                                                                                                                                  <  2:4488025/150‑1 (MQ=255)
                      cATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGgag                                                                                                                                  >  2:1212393/1‑150 (MQ=255)
                        ttAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGaa                                                                                                                                <  1:1125999/150‑1 (MQ=255)
                          aaCTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCACCTTTTGCGCGGGCATTTATTCCCCCCAGTTGTTGTTGGATATGCCCCTCGCGCTGGAGaact                                                                                                                              >  1:2882700/1‑148 (MQ=255)
                          aaCTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAat                                                                                                                              >  2:3435454/1‑150 (MQ=255)
                          aaCTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAat                                                                                                                              <  1:2994944/150‑1 (MQ=255)
                           aCTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAata                                                                                                                             <  2:3774009/150‑1 (MQ=255)
                           aCTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAata                                                                                                                             >  1:3289402/1‑150 (MQ=255)
                            cTGCGGTCCCTCCGGCCAGGCCCAGACAACAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCGGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTAGTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAatat                                                                                                                            <  2:556718/150‑1 (MQ=255)
                             tGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATAtt                                                                                                                           <  2:845479/150‑1 (MQ=255)
                               cGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTc                                                                                                                         <  2:4572487/150‑1 (MQ=255)
                               cGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTc                                                                                                                         <  1:3127/150‑1 (MQ=255)
                                ggTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTcc                                                                                                                        >  1:1907579/1‑150 (MQ=255)
                                 gTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTccc                                                                                                                       <  1:4560948/150‑1 (MQ=255)
                                  tCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTcccc                                                                                                                      <  2:3068078/150‑1 (MQ=255)
                                    cATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCt                                                                                                                    >  1:23155/1‑150 (MQ=255)
                                        cGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTAGTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAg                                                                                                                <  2:1993805/150‑1 (MQ=255)
                                         ggCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGc                                                                                                               >  1:4304227/1‑150 (MQ=255)
                                             aTGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTAc                                                                                                           <  1:1439259/150‑1 (MQ=255)
                                              tGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACt                                                                                                          <  1:1533979/150‑1 (MQ=255)
                                               gCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTg                                                                                                         >  1:3356223/1‑150 (MQ=255)
                                               gCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTg                                                                                                         >  2:4638851/1‑150 (MQ=255)
                                                cccAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGt                                                                                                        <  2:2322565/150‑1 (MQ=255)
                                                     aCAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCt                                                                                                   >  1:2541927/1‑150 (MQ=255)
                                                          aaCGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAgg                                                                                              >  1:2681107/1‑150 (MQ=255)
                                                            cGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGcc                                                                                            >  2:2641782/1‑150 (MQ=255)
                                                            cGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGcc                                                                                            >  1:4649968/1‑150 (MQ=255)
                                                                cAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCt                                                                                        <  1:2801184/150‑1 (MQ=255)
                                                                cAGAGAATTCAGCGTGTCGCTATCGACCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCATCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCt                                                                                        >  2:1101540/1‑150 (MQ=255)
                                                                 agagAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTa                                                                                       <  2:378473/150‑1 (MQ=255)
                                                                 agagAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTa                                                                                       >  2:1950383/1‑150 (MQ=255)
                                                                 agagAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTa                                                                                       <  2:376692/150‑1 (MQ=255)
                                                                   agaATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAAc                                                                                     <  2:2102500/150‑1 (MQ=255)
                                                                   agaATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAAc                                                                                     <  1:1065383/150‑1 (MQ=255)
                                                                    gaATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACg                                                                                    <  1:480144/150‑1 (MQ=255)
                                                                       ttCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGtt                                                                                 <  1:5060723/150‑1 (MQ=255)
                                                                          aGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGc                                                                              >  1:2153134/1‑150 (MQ=255)
                                                                               gtCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGcac                                                                         >  1:163737/1‑150 (MQ=255)
                                                                                tCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGcaca                                                                        >  2:1814164/1‑150 (MQ=255)
                                                                                  gCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGc                                                                      <  1:3801044/150‑1 (MQ=255)
                                                                                   cTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCg                                                                    <  1:447215/150‑1 (MQ=255)
                                                                                    tATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCgg                                                                   >  1:2071668/1‑150 (MQ=255)
                                                                                    tATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCgg                                                                   >  1:846617/1‑150 (MQ=255)
                                                                                    tATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCgg                                                                   >  2:2357677/1‑150 (MQ=255)
                                                                                         tCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAAtt                                                              <  2:1008422/150‑1 (MQ=255)
                                                                                          ccTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTa                                                             >  2:3954395/1‑150 (MQ=255)
                                                                                              gcgcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCt                                                         >  2:2631686/1‑150 (MQ=255)
                                                                                              gcgcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCt                                                         <  2:2430400/150‑1 (MQ=255)
                                                                                               cgcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTa                                                        <  2:361677/150‑1 (MQ=255)
                                                                                                gcATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTaa                                                       <  1:3043391/150‑1 (MQ=255)
                                                                                                  atcatcAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACt                                                     >  2:328020/1‑150 (MQ=255)
                                                                                                   tcatcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACtt                                                    <  2:714014/150‑1 (MQ=255)
                                                                                                     atcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTc                                                  <  2:4860856/150‑1 (MQ=255)
                                                                                                     atcaAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTc                                                  <  2:1474191/150‑1 (MQ=255)
                                                                                                       caAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTg                                                <  1:406455/150‑1 (MQ=255)
                                                                                                       caAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTg                                                >  1:3031719/1‑150 (MQ=255)
                                                                                                        aaGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGc                                               >  1:3661860/1‑150 (MQ=255)
                                                                                                        aaGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCACTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGc                                               >  2:2263140/1‑150 (MQ=255)
                                                                                                        aaGCCAGCAACTTTTGCGCGGGCATTTAGTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGc                                               <  2:4114719/150‑1 (MQ=255)
                                                                                                         aGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCg                                              >  1:1371142/1‑150 (MQ=255)
                                                                                                          gCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGt                                             <  2:959808/150‑1 (MQ=255)
                                                                                                          gCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGt                                             >  1:3041754/1‑150 (MQ=255)
                                                                                                            cAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTcc                                           >  2:2045978/1‑150 (MQ=255)
                                                                                                               cAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGt                                        >  1:5200968/1‑150 (MQ=255)
                                                                                                                   ttttGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTc                                    <  2:1337812/150‑1 (MQ=255)
                                                                                                                            ggCATTTATTCCCCCCCTTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCAGAAGCGTAAGGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCa                            <  2:2849070/150‑1 (MQ=255)
                                                                                                                            ggCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAg                           <  1:3392376/150‑1 (MQ=255)
                                                                                                                             gCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAgg                          <  1:4285911/150‑1 (MQ=255)
                                                                                                                              cATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGc                         >  1:3505415/1‑150 (MQ=255)
                                                                                                                              cATTTATTCCCCCCATTTGTTGTCGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGc                         >  1:1856227/1‑150 (MQ=255)
                                                                                                                               aTTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGcc                        <  1:3804746/150‑1 (MQ=255)
                                                                                                                                  tATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCAtt                     <  1:5019545/150‑1 (MQ=255)
                                                                                                                                   aTTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTa                    <  1:4921250/150‑1 (MQ=255)
                                                                                                                                   aTTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTa                    >  1:1018591/1‑150 (MQ=255)
                                                                                                                                     tCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTAcc                  <  1:1679330/150‑1 (MQ=255)
                                                                                                                                     tCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGACACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTAcc                  <  1:88844/150‑1 (MQ=255)
                                                                                                                                       cccccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACcgc                >  2:4089586/1‑150 (MQ=255)
                                                                                                                                          ccATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTATGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGAc             >  2:2726264/1‑150 (MQ=255)
                                                                                                                                           cATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACt            >  2:1609282/1‑150 (MQ=255)
                                                                                                                                           cATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACt            >  1:1190977/1‑150 (MQ=255)
                                                                                                                                            aTTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTg           >  2:524250/1‑150 (MQ=255)
                                                                                                                                                  tgttgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGACGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGAGGAATTTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCt     <  2:634554/150‑1 (MQ=255)
                                                                                                                                                   gttgGATATGCCCCTCGCGCTGGTGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCtt    >  2:1738650/1‑150 (MQ=255)
                                                                                                                                                   gttgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTACAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCtt    <  2:4787844/150‑1 (MQ=255)
                                                                                                                                                     tgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTTg  >  1:2136727/1‑150 (MQ=255)
                                                                                                                                                     tgGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTTg  <  2:4003124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCAGCTCTTTTAAACGATCAGACATTAACTGCGGTCCATCCGGCCATGCCCAGACAAAAACGCCCAGAGAATTCAGCGTGTCGCTATCGTCCTGGCGCATCATCAAGCCAGCAACTTTTGCGCGGGCATTTATTCCCCCCATTTGTTGTTGGATATGCCCCTCGCGCTGGAGAAATATTTCCCCCTGAAGCGTACTGTTGGCTAAAGGCCCGCTAACGGTTTGCGGCACAGCGGAAATTAAGCTAACTTTCTGCGTCCAGTCCTCGCCTTTCAGGCCATTACCGCGACTGATTGCTTTG  >  CP009273/3813813‑3814111

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: