Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,621,021 Δ1 bp 100% coding (31/417 nt) PP_5654 → phosphotriesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,621,0210A.89.1% 194.4 / 17.7 55coding (31/417 nt)PP_5654phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base . (24/25);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAA  >  NC_002947/4620878‑4621163
                                                                                                                                               |                                                                                                                                               
tctctcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaa                                                                                                                                             >  8:97747/1‑148 (MQ=255)
tctctcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaa                                                                                                                                             >  2:122088/1‑148 (MQ=255)
    tcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                                         >  6:290414/1‑148 (MQ=255)
      tcccTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                                         <  4:212600/144‑1 (MQ=255)
      gACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATACGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaa                                                                                                                                      <  5:290414/149‑4 (MQ=255)
        ccTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaaa                                                                                                                                     <  2:338520/148‑5 (MQ=255)
         cTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCACAAGCAAAACTCACACAAACTACCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGATAACATCTGCAGATGGTCCaaaaaaaaacaacg                                                                                                                                    >  4:240275/1‑143 (MQ=255)
          ttCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATATGGAGATGGTCCATAAAAAAACAAAGa                                                                                                                                  >  6:30580/1‑148 (MQ=255)
               agaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaa                                                                                                                                            >  3:89597/1‑134 (MQ=255)
               agaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaa                                                                                                                                            <  4:89597/134‑1 (MQ=255)
               agaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACAt                                                                                                                            >  2:172257/1‑149 (MQ=255)
                gaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACAt                                                                                                                            >  2:20827/1‑148 (MQ=255)
                       aGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAAc                                                                                                                    <  4:161270/149‑1 (MQ=255)
                          tAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAGCATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACAt                                                                                                                 <  6:144522/148‑1 (MQ=255)
                            gagCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATaaa                                                                                                              <  7:102480/149‑1 (MQ=255)
                            gagCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATaa                                                                                                               >  8:102480/1‑148 (MQ=255)
                             agCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATaaaa                                                                                                             <  8:199376/149‑1 (MQ=255)
                                                gCAAACAATAACTCCCAAGCAAAACTCACTCAATCTCTCGTAAAAATAATAAGTGTAACATGCAAATCAACAAAGCTAACATCTTGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCt                                                                                          <  8:292309/149‑1 (MQ=255)
                                                  aaaCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCtt                                                                                         >  3:190752/1‑148 (MQ=255)
                                                  aaaCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCtt                                                                                         >  8:87612/1‑148 (MQ=255)
                                                    aCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGAAGGAGATATTTTATGGTTCAc                                                                                      >  4:188478/1‑149 (MQ=255)
                                                     cAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACACTAAAAGCAGGAGATATTTTCTGCTTCAc                                                                                      <  1:20827/149‑1 (MQ=255)
                                                        tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAg                                                                                                         <  2:185506/126‑1 (MQ=255)
                                                        tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAg                                                                                                         >  1:185506/1‑126 (MQ=255)
                                                        tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAAAGAACAAACCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACAAAAAATAACATAAAACCATCACATATTTTCTGCTTCAACCAg                                                                                  >  3:292748/1‑149 (MQ=255)
                                                              ccAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGAc                                                                            >  5:60575/1‑149 (MQ=255)
                                                               cAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGCACAAAGCTAACATCTGGAGATTTGCTAAACAAAAACAAATACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGTAAGACg                                                                           <  5:242586/149‑1 (MQ=255)
                                                               cAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACg                                                                           >  5:59116/1‑149 (MQ=255)
                                                                  gCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGa                                                                         >  4:105137/1‑148 (MQ=255)
                                                                      aaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACata                                                                    <  7:97747/149‑1 (MQ=255)
                                                                           acacaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATattatt                                                               >  1:31261/1‑149 (MQ=255)
                                                                             acaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTAtttt                                                             <  3:227695/149‑1 (MQ=255)
                                                                              aaaaCTCCCGTAAAAATCATAATTTTAACATTCAAATTAAAAAATCAAACATCTGGAGATGTTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTg                                                            <  3:188478/148‑1 (MQ=255)
                                                                               aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAATGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                           <  3:105137/149‑1 (MQ=255)
                                                                               aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                           <  7:87612/149‑1 (MQ=255)
                                                                               aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                           <  1:172257/149‑1 (MQ=255)
                                                                                   tCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGc                                                         >  5:293770/1‑148 (MQ=255)
                                                                                     ccGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCGGGAAGACGGAACATATTATTTTGGGCAaat                                                     <  4:190752/149‑1 (MQ=255)
                                                                                          aaaCTCAGAAGTTGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACAt                                                                                                                 <  5:77941/84‑1 (MQ=255)
                                                                                          aaaCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAatat                                                <  3:102685/149‑1 (MQ=255)
                                                                                          aaaCTCAGAAGTGGAACATGCAAAAGAACAAAGCTAACATCTGGAGATGGACCAAAAAAAAAAAAACAATACATTAAAGAAAAt                                                                                                                 >  6:77941/1‑84 (MQ=255)
                                                                                              tCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaa                                             >  4:253102/1‑148 (MQ=255)
                                                                                               cAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCa                                              >  8:258730/1‑146 (MQ=255)
                                                                                               cAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCa                                              <  7:258730/146‑1 (MQ=255)
                                                                                                          cATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTGACCCAGGAAGACGGAAAATATTTTTTTTGGCAAATAATAACAAAAATAATAACttg                                >  5:118829/1‑147 (MQ=255)
                                                                                                           aTGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGc                               >  7:26513/1‑149 (MQ=255)
                                                                                                            tGCAAATGAACAAAGCTAACATCTGGAGATGGTCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGc                               >  1:139135/1‑148 (MQ=255)
                                                                                                             gCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCa                             >  4:271496/1‑149 (MQ=255)
                                                                                                               aaaTGAACAAAGCTAACATCTGGATATGGTCCAAAAAAAAACAAAGAATACATTAAAGAACATAAAAGCAGTAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                           <  3:240275/149‑1 (MQ=255)
                                                                                                                    aaCAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGCAACATATTATTTTGGCAAAATAATATCAAAAATAATAAATGGCCATTTaca                       >  4:123312/1‑146 (MQ=255)
                                                                                                                     aCAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCt                      <  6:293770/149‑1 (MQ=255)
                                                                                                                      cAAAGCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGACATGTAGCTg                     >  5:52476/1‑148 (MQ=255)
                                                                                                                        aaGCTAACATCTCTAGATGGTCCAAAAAAAAACAAAGTATACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACTGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGaaa                  <  4:292748/149‑1 (MQ=255)
                                                                                                                          gCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAata                >  5:255258/1‑149 (MQ=255)
                                                                                                                          gCTAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAACCACGACATATTTTCAGCTTCACCCACGAAGACGGAACATATAATTTTGGGCAAATAATATCAAAAATAATAAAAGGCCATGTAGCTGAAat                 >  5:259888/1‑148 (MQ=255)
                                                                                                                            tAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAatat               >  2:328186/1‑148 (MQ=255)
                                                                                                                            tAACATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATAtt              <  3:123312/149‑1 (MQ=255)
                                                                                                                              aCATCTGGAGTTGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATAtt              >  4:259025/1‑147 (MQ=255)
                                                                                                                              aCATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCa            <  1:328186/149‑1 (MQ=255)
                                                                                                                               cATCTTGTGATGTTCCAAAAAAAAACAAAGACTACATTTAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAg           <  6:259888/149‑1 (MQ=255)
                                                                                                                               cATCTGGAGATGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAg           <  6:255258/149‑1 (MQ=255)
                                                                                                                                      agaTGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATgg      >  6:162553/1‑148 (MQ=255)
                                                                                                                                        aTGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa  <  7:240630/149‑1 (MQ=255)
                                                                                                                                         tGGTCCAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa  <  3:253102/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAA  >  NC_002947/4620878‑4621163

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: