Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,624,546 A→T 7.0% intergenic (+3139/+1425) 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,5460AT 7.0% 199.0 / 5.2 71intergenic (+3139/+1425)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base T (3/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

AAATTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTCTATAAACAT  >  NC_002947/4624398‑4624694
                                                                                                                                                     |                                                                                                                                                    
aaaTTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTaa                                                                                                                                                      >  2:202838/1‑149 (MQ=255)
aaaTTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTaa                                                                                                                                                      >  2:241009/1‑149 (MQ=255)
         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                                                                                         <  5:51547/149‑1 (MQ=255)
                                                             tACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                         <  3:346372/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:241516/6‑122 (MQ=255)
                                                                                                                                        gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                          >  3:236662/6‑122 (MQ=255)
                                                                                                                                          gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                                  >  7:171642/1‑128 (MQ=255)
                                                                                                                                          gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                                  <  8:171642/128‑1 (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)
                                                                                                                                                     |                                                                                                                                                    
AAATTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCTCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTCTATAAACAT  >  NC_002947/4624398‑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: