Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,624,548 G→A 6.6% intergenic (+3141/+1423) PP_5654 → / ← PP_4091 phosphotriesterase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,624,5480GA 6.6% 218.7 / 3.9 75intergenic (+3141/+1423)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base A (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-02

GTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTCTATAAACAT  >  NC_002947/4624407‑4624694
                                                                                                                                              |                                                                                                                                                  
gTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACATCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTg                                                                                                                                             >  7:43015/1‑149 (MQ=255)
    tAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGttt                                                                                                                                          <  4:296255/148‑1 (MQ=255)
     aGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTc                                                                                                                                         >  5:137101/1‑148 (MQ=255)
      gATTATGTCCCCTTAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCaa                                                                                                                                       <  5:206964/149‑1 (MQ=255)
       ataaagTCCCCAGAGAAAAAAATAATCGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAg                                                                                                                                      <  6:306129/144‑1 (MQ=255)
       aTTATGTCCCCTGAGAAAAAAATAATTGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCaa                                                                                                                                       <  3:212643/148‑1 (MQ=255)
        ttATGTCCCCTGAGAAAAAAATAATTGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGc                                                                                                                                     <  4:303222/149‑1 (MQ=255)
              ccccTGAGAAAAAGATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaa                                                                                                                                >  8:228545/1‑148 (MQ=255)
               cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGGTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaac                                                                                                                               >  3:56296/1‑148 (MQ=255)
               cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaac                                                                                                                               >  6:51547/1‑148 (MQ=255)
               cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCAAAAACCTCATGACATTGAAATGCTGAATAGAGAAATTTACAAATGGCTTTGAAATAAAAAAAAAGAAAGCAGGAAAATAAAAATGAAATTTTTTAAAAGcttaaaa                                                                                                                               >  8:239692/1‑142 (MQ=255)
                ccTGAGAAAAAAATAATGGCTACCATTAGAAAACGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaacaa                                                                                                                             <  4:86043/149‑1 (MQ=255)
                 cTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCAGCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGTAATTTGTTTCAAGCGTTaaca                                                                                                                              >  1:281157/1‑148 (MQ=255)
                    agaAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACt                                                                                                                           <  2:186635/147‑1 (MQ=255)
                     gaAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCa                                                                                                                         <  5:62307/148‑1 (MQ=255)
                        aaaatTAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACAATGAAATGCTGCCTAGAGAAAGTTTAACTTGGCTTTCCAAAGACAATCTAGCTAGCAGGGAATTTGTAATGGATTTTGTTTCAAGCGTTAACAACTCAAAAt                                                                                                                     >  2:89545/1‑149 (MQ=255)
                         aaaaTAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAAtt                                                                                                                    >  8:245711/1‑149 (MQ=255)
                             taatGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTATCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCt                                                                                                                <  8:258443/149‑1 (MQ=255)
                                   cTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTTGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGt                                                                                                           >  3:358513/1‑148 (MQ=255)
                                      ccATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCAt                                                                                                        <  7:245711/148‑1 (MQ=255)
                                           aGAAACCGATACCGAAGACCGCATGACATTGAAATGCTGGCTCGCGAACTTTCCACATGGCTTTGCAATGACAAACTAGCTCGCAGGTAAATTCTAATGGAATTTGTTTCAAGCGTTAACTACTCAAAATTCCCTGACGTCTTTCTAc                                                                                                   <  4:269732/148‑1 (MQ=255)
                                             aaaaGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTaa                                                                                                >  7:285541/1‑149 (MQ=255)
                                             aaaaGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTaa                                                                                                >  4:240956/1‑149 (MQ=255)
                                             aaaaGGATACCGAAGACCTAATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCACTGACGTCATTCAACTaa                                                                                                >  8:49998/1‑149 (MQ=255)
                                              aaaGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAAc                                                                                               <  4:282192/149‑1 (MQ=255)
                                                  gATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGATATTTACTCATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCttt                                                                                           <  4:56296/149‑1 (MQ=255)
                                                   aTACCGAAGACCTAATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCttt                                                                                           <  7:239692/148‑1 (MQ=255)
                                                    tACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTg                                                                                          >  4:346372/1‑148 (MQ=255)
                                                    tACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                         <  8:285541/149‑1 (MQ=255)
                                                    tACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                         <  3:346372/149‑1 (MQ=255)
                                                    tACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                         <  5:51547/149‑1 (MQ=255)
                                                       cGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGtt                                                                                                                                  >  4:304843/1‑105 (MQ=255)
                                                       cGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGtt                                                                                                                                  <  3:304843/105‑1 (MQ=255)
                                                         aaGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACt                                                                                    <  2:80875/149‑1 (MQ=255)
                                                           gACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTg                                                                                  <  1:189454/149‑1 (MQ=255)
                                                              cTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAg                                                                               >  8:185119/1‑149 (MQ=255)
                                                                       aTTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATc                                                                       >  4:9106/1‑149 (MQ=255)
                                                                            aaTGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTc                                                                  <  7:49998/149‑1 (MQ=255)
                                                                             aTGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCa                                                                 <  4:279410/149‑1 (MQ=255)
                                                                                 tGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGa                                                              >  5:28435/1‑148 (MQ=255)
                                                                                 tGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGa                                                              >  5:34020/1‑148 (MQ=255)
                                                                                    cTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGTGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGc                                                          <  2:281157/149‑1 (MQ=255)
                                                                                         gAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGt                                                      <  1:241009/148‑1 (MQ=255)
                                                                                         gAAATTTACACATGGATTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCGCTGAAGTCATTCAACTAACCTTTGAATACTATAAGCTCCTATCCAGTCATGAAACCTAGt                                                      >  4:313360/1‑148 (MQ=255)
                                                                                          aaaTTTACACATGGCGTTGGAAAGACAAACTAGCTAGCAGGGAAATTGTTATGGAATTTGTTTGAAGCGGTTACAACTCAAAATTACCGGACGTCATTCAACTTACCTTTGAATACTTGAAGCGCCTTTCCACTACAGCAAGCGAGttt                                                    >  4:368466/1‑148 (MQ=255)
                                                                                            aTTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACt                                                  <  7:185119/149‑1 (MQ=255)
                                                                                              tgagggaTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACtat                                                <  1:79198/143‑1 (MQ=255)
                                                                                                    aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACTTCATTCAACTAAACTTTGAATACTTTAAGCGACTATCCACTCATGAAAGAGAGTTACTATACgat                                           >  5:389461/1‑147 (MQ=255)
                                                                                                    aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGAGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                           >  4:368529/1‑148 (MQ=255)
                                                                                                    aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTAATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTAATATACgag                                           >  5:195568/1‑148 (MQ=255)
                                                                                                    aTGGCTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                           >  6:69262/1‑147 (MQ=255)
                                                                                                     tGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                           >  6:370172/1‑147 (MQ=255)
                                                                                                       gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagga                                         >  6:290681/1‑147 (MQ=255)
                                                                                                       gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagga                                         <  5:290681/147‑1 (MQ=255)
                                                                                                       gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                                       <  1:202838/149‑1 (MQ=255)
                                                                                                       gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                                       <  5:370172/149‑1 (MQ=255)
                                                                                                        gcttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                                       <  5:69262/146‑1 (MQ=255)
                                                                                                         tttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACtt                                                                                   <  6:329486/103‑1 (MQ=255)
                                                                                                         tttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACtt                                                                                   >  5:329486/1‑103 (MQ=255)
                                                                                                                    aaaCTAGCTACCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAgg                          >  4:80045/1‑149 (MQ=255)
                                                                                                                          gatagcAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCTCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCTaaaa                                >  6:286685/3‑137 (MQ=255)
                                                                                                                             ttgaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt                 >  5:58609/8‑149 (MQ=255)
                                                                                                                              gcAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt                 <  7:228545/148‑1 (MQ=255)
                                                                                                                               gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                          <  4:241516/117‑1 (MQ=255)
                                                                                                                               gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                          <  4:236662/117‑1 (MQ=255)
                                                                                                                               gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                          >  3:236662/6‑122 (MQ=255)
                                                                                                                               gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                          >  3:241516/6‑122 (MQ=255)
                                                                                                                                 gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                                  <  8:171642/128‑1 (MQ=255)
                                                                                                                                 gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                                  >  7:171642/1‑128 (MQ=255)
                                                                                                                                  ggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGAtt            <  4:358513/149‑1 (MQ=255)
                                                                                                                                   gAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTc           <  6:137101/149‑1 (MQ=255)
                                                                                                                                           aaTGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt                 >  7:297863/1‑135 (MQ=255)
                                                                                                                                           aaTGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt                 <  8:297863/135‑1 (MQ=255)
                                                                                                                                            aTGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAATAAATAGGGGATTTATTTGATGCTATAAAcat  >  8:349408/1‑149 (MQ=255)
                                                                                                                                            aTGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTCTATaaa     <  3:240956/146‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTCTATAAACAT  >  NC_002947/4624407‑4624694

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: