Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,001,304 (G)7→8 coding (960/2259 nt) xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,001,2971.G95.7% 238.7 / 2.3 69coding (953/2259 nt)xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (34/32);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001161‑3001427
                                                                                                                                         |                                                                                                                                  
gCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg                                                                                                                                  <  2:485149/139‑1 (MQ=255)
   aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggg                                                                                                                               <  2:1284170/139‑1 (MQ=255)
   aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggg                                                                                                                               <  2:592844/139‑1 (MQ=255)
          cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggggaata                                                                                                                        >  2:400580/1‑134 (MQ=255)
          cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa                                                                                                                        <  1:1663995/139‑1 (MQ=255)
              ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGt                                                                                                                     >  1:402348/1‑139 (MQ=255)
              ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa                                                                                                                        <  2:961236/135‑1 (MQ=255)
              ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa                                                                                                                        >  1:961236/1‑135 (MQ=255)
                  gCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 <  2:432140/139‑1 (MQ=255)
                            ttGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGTTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATc                                                                                                       <  2:1011198/139‑1 (MQ=255)
                             tGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATcc                                                                                                      >  2:1360208/1‑139 (MQ=255)
                              gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTGTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATaa                                                                                                                        <  2:539670/120‑1 (MQ=255)
                              gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATaa                                                                                                                        >  1:539670/1‑120 (MQ=255)
                                   gCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtg                                                                                                >  1:472193/1‑139 (MQ=255)
                                    cAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgt                                                                                               >  1:311357/1‑139 (MQ=255)
                                     aaaTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              <  2:402348/139‑1 (MQ=255)
                                                 aCCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATc                                                                                                       <  2:408368/118‑1 (MQ=255)
                                                 aCCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATc                                                                                                       >  1:408368/1‑118 (MQ=255)
                                                     cgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTc                                                                              >  1:218852/1‑139 (MQ=255)
                                                              aacatccagTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTg                                                                     >  1:878372/9‑139 (MQ=255)
                                                                  aaaGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCTTCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCtat                                                                 >  2:1226981/1‑139 (MQ=255)
                                                                    aGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGc                                                                                                                  <  1:19844/88‑1 (MQ=255)
                                                                    aGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGc                                                                                                                  >  2:19844/1‑88 (MQ=255)
                                                                     ggTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            <  1:699214/109‑1 (MQ=255)
                                                                     ggTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            >  2:699214/1‑109 (MQ=255)
                                                                       ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATc                                                                                                       <  2:343458/96‑1 (MQ=255)
                                                                       ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATc                                                                                                       >  1:343458/1‑96 (MQ=255)
                                                                       ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGAc                                                                        <  2:1249964/127‑1 (MQ=255)
                                                                       ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGAc                                                                        >  1:1249964/1‑127 (MQ=255)
                                                                        tataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                           <  1:400580/139‑1 (MQ=255)
                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTccc                                                      >  1:654887/1‑136 (MQ=255)
                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTccc                                                      <  2:654887/136‑1 (MQ=255)
                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTc                                                   >  2:369588/1‑139 (MQ=255)
                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggggaata                                                                                                                        >  1:60943/1‑64 (MQ=255)
                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggggaata                                                                                                                        <  2:60943/69‑6 (MQ=255)
                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCt                                                                                                            >  1:475542/1‑81 (MQ=255)
                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCt                                                                                                            <  2:475542/81‑1 (MQ=255)
                                                                                    gTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCt                                               <  1:1360208/139‑1 (MQ=255)
                                                                                       cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATaa                                                                                                                        <  1:557084/63‑1 (MQ=255)
                                                                                       cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATaa                                                                                                                        >  2:557084/1‑63 (MQ=255)
                                                                                       cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                >  1:572198/1‑135 (MQ=255)
                                                                                       cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                <  2:572198/135‑1 (MQ=255)
                                                                                              acacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGc                                     <  1:1226981/139‑1 (MQ=255)
                                                                                               cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  >  1:1020484/1‑77 (MQ=255)
                                                                                               cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAATCGCTTACCTTTATCCTCGt                                                                                                  <  2:1020484/77‑1 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATcc                                                                                                      <  1:1294473/68‑1 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATcc                                                                                                      >  2:1294473/1‑68 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         <  2:1470087/129‑1 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         >  1:1470087/1‑129 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      >  1:402342/1‑132 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  1:1254512/132‑1 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      >  2:1254512/1‑132 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  2:402342/132‑1 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                               >  1:584679/1‑139 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                               >  2:481282/1‑139 (MQ=255)
                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                               >  2:209943/1‑139 (MQ=255)
                                                                                                         tATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCCGGCTGGTGCGATGCGTGGTTATGgcg                         <  2:218852/139‑1 (MQ=255)
                                                                                                               tctcACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTAcc                                                                                                             >  2:1356725/1‑50 (MQ=255)
                                                                                                               tctcACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTAcc                                                                                                             <  1:1356725/50‑1 (MQ=255)
                                                                                                                  cACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACaa                 <  1:209943/139‑1 (MQ=255)
                                                                                                                     gggCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATgg                           >  2:444224/1‑126 (MQ=255)
                                                                                                                     gggCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATgg                           <  1:444224/126‑1 (MQ=255)
                                                                                                                        cACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTa           >  2:514456/1‑139 (MQ=255)
                                                                                                                         aCTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAt          <  1:481282/139‑1 (MQ=255)
                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 <  1:368026/35‑1 (MQ=255)
                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 >  2:368026/1‑35 (MQ=255)
                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCg                                    <  2:339570/112‑1 (MQ=255)
                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCg                                    >  1:339570/1‑112 (MQ=255)
                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAtt         >  2:1564708/1‑139 (MQ=255)
                                                                                                                             cATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGc      <  2:472193/139‑1 (MQ=255)
                                                                                                                               tCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCg    >  1:864445/1‑139 (MQ=255)
                                                                                                                                acggtcggCGGGGGGGACTAAGGTCGCTTAGCTTTATCGTCGTGGTCCCTACGCTTACGGTTCAAAGACCTGCTAAACCAACTGCCCCTCGGCTGGTGCGAGGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  2:407702/132‑1 (MQ=255)
                                                                                                                                 gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      >  1:427667/1‑103 (MQ=255)
                                                                                                                                 gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  2:427667/103‑1 (MQ=255)
                                                                                                                                 gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  2:311357/139‑1 (MQ=255)
                                                                                                                                 gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  2:688035/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001161‑3001427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: