Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,287 C→T *257* (TAG→TAA yqeC ← putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,2870CT98.9% 306.6 / NA 94*257* (TAG→TAAyqeCputative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (44/49);  minor base G (0/1);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACG  >  NC_000913/3014150‑3014418
                                                                                                                                         |                                                                                                                                   
gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                   <  1:682103/139‑1 (MQ=255)
gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                   <  2:1077422/139‑1 (MQ=255)
     tCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTca                                                                                                                              >  2:672400/1‑139 (MQ=255)
             cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                                                      >  2:659065/1‑139 (MQ=255)
             cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                                                      >  1:166061/1‑139 (MQ=255)
                 tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  1:1394538/139‑1 (MQ=255)
                 tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  2:1986242/139‑1 (MQ=255)
                   tCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGcgcgcg                                                                                                                >  1:1194606/1‑139 (MQ=255)
                        aaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGc                                                                                                           >  1:1976260/1‑139 (MQ=255)
                         aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                          <  2:166061/139‑1 (MQ=255)
                          tttGCGCCATAGCAAACGGACAATTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  >  2:194906/1‑129 (MQ=255)
                          tttGCGCCATAGCAAACGGACAATTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  1:194906/129‑1 (MQ=255)
                            tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  1:29770/139‑1 (MQ=255)
                             gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  >  2:1153212/1‑127 (MQ=255)
                             gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  1:1153212/127‑1 (MQ=255)
                              cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGTGCGCGATTGCAGGATg                                                                                                     <  2:551232/139‑1 (MQ=255)
                              cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  2:1418955/139‑1 (MQ=255)
                              cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  2:1316899/139‑1 (MQ=255)
                              cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  1:672400/139‑1 (MQ=255)
                                  aTAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCt                                                                                                 <  2:1194606/139‑1 (MQ=255)
                                    aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTg                                                                                               >  1:1239990/1‑139 (MQ=255)
                                    aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTg                                                                                               >  1:1104252/1‑139 (MQ=255)
                                       aaaCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATTCCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            <  1:1374900/139‑1 (MQ=255)
                                           ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  <  2:1050526/113‑1 (MQ=255)
                                           ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                  >  1:1050526/1‑113 (MQ=255)
                                           ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                                        <  2:1976260/139‑1 (MQ=255)
                                             aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGcc                                                                                      >  2:552109/1‑139 (MQ=255)
                                             aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGcc                                                                                      >  2:695485/1‑139 (MQ=255)
                                               aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcgc                                                                                                                 >  1:144691/1‑110 (MQ=255)
                                               aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcgc                                                                                                                 <  2:144691/110‑1 (MQ=255)
                                               aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCcca                                                                                    <  1:552109/139‑1 (MQ=255)
                                                     gCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCCTACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                                                                              <  2:1239990/139‑1 (MQ=255)
                                                     gCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                                                                              <  1:1245918/139‑1 (MQ=255)
                                                      ctctTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                                        <  2:1952523/128‑1 (MQ=255)
                                                      ctctTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                                        >  1:1952523/1‑128 (MQ=255)
                                                           cggtttttGCACTCTCTGTGAAAAGATTCACACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGTCCGTTAATTCACAAAGCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                                        <  1:1227984/135‑1 (MQ=255)
                                                           cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                                        >  2:312217/1‑139 (MQ=255)
                                                                  tGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAAt                                                                                                                                 <  2:1819551/75‑1 (MQ=255)
                                                                  tGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAAt                                                                                                                                 >  1:1819551/1‑75 (MQ=255)
                                                                           tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:10882/92‑1 (MQ=255)
                                                                           tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       >  1:10882/1‑92 (MQ=255)
                                                                            gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                          >  2:1508554/1‑88 (MQ=255)
                                                                            gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                          <  1:1508554/88‑1 (MQ=255)
                                                                                 aaTATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                      >  2:1179361/1‑135 (MQ=255)
                                                                                 aaTATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                      <  1:1179361/135‑1 (MQ=255)
                                                                                  atatTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTACAGGATg                                                                                                     <  2:247037/87‑1 (MQ=255)
                                                                                  atatTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTACAGGATg                                                                                                     >  1:247037/1‑87 (MQ=255)
                                                                                     ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                   >  1:274330/1‑54 (MQ=255)
                                                                                     ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:1620069/82‑1 (MQ=255)
                                                                                     ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       >  1:1620069/1‑82 (MQ=255)
                                                                                     ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                      <  2:1864137/83‑1 (MQ=255)
                                                                                     ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                      >  1:1864137/1‑83 (MQ=255)
                                                                                      tCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGaa                                             <  1:1382818/139‑1 (MQ=255)
                                                                                          aCCAACCTCAATGTATATATTTCAGTAACCCGGAATACCCGGGCCTGTAATTCAGAAATCTGCGTGCGATTGCAGGATGTTCGTGGATATCGCCGAGCCAGATTGCCTCTACTTCGTGTTGTTGCAGCGGTTGAAGGAg                                         <  1:632376/139‑1 (MQ=255)
                                                                                             aaCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagccag                                                                               <  2:57081/98‑1 (MQ=255)
                                                                                             aaCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagccag                                                                               >  1:57081/1‑98 (MQ=255)
                                                                                              aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                     >  1:1433520/1‑139 (MQ=255)
                                                                                              aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                     >  2:1576920/1‑139 (MQ=255)
                                                                                               ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGc                                    >  2:1944366/1‑139 (MQ=255)
                                                                                                cTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                   <  1:312217/139‑1 (MQ=255)
                                                                                                     tGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgtg                              <  2:1433520/139‑1 (MQ=255)
                                                                                                       gATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCa                            <  1:1944366/139‑1 (MQ=255)
                                                                                                         atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            >  1:188842/5‑74 (MQ=255)
                                                                                                         atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                            <  2:188842/70‑1 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:1595864/59‑1 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       >  1:1595864/3‑61 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     >  2:857024/3‑63 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                     <  1:857024/61‑1 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGaa                                             <  2:594524/117‑1 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGaa                                             >  1:594524/3‑119 (MQ=255)
                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                         <  1:659065/137‑1 (MQ=255)
                                                                                                              tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                      <  2:1104252/139‑1 (MQ=255)
                                                                                                                  aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagcc                                                                                 <  2:161244/76‑1 (MQ=255)
                                                                                                                  aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagcc                                                                                 >  1:161244/1‑76 (MQ=255)
                                                                                                                  aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCa                            <  2:944239/129‑1 (MQ=255)
                                                                                                                  aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCa                            >  1:944239/1‑129 (MQ=255)
                                                                                                                      aCCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    <  2:1736153/37‑1 (MQ=255)
                                                                                                                      aCCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    >  1:1736153/1‑37 (MQ=38)
                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgc                                                                                                                   >  1:1733264/1‑37 (MQ=255)
                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTgcgc                                                                                                                   <  2:1733264/37‑1 (MQ=255)
                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    >  1:958527/1‑36 (MQ=255)
                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                    <  2:958527/36‑1 (MQ=255)
                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACAt             >  1:96101/1‑139 (MQ=255)
                                                                                                                                cccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                        <  2:21612/119‑1 (MQ=255)
                                                                                                                                cccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                        >  1:21612/1‑119 (MQ=255)
                                                                                                                                cccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaaaaa    >  2:1270068/1‑139 (MQ=255)
                                                                                                                                 ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                        >  2:515410/1‑102 (MQ=255)
                                                                                                                                 ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                        <  1:515410/102‑1 (MQ=255)
                                                                                                                                 ccGGGCCGTTAATCCACAAATCTGCGCGCGATTGCGGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTGCGTCGTGTTTTTACAGCGTTTGTAGAAGATCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAAc   >  1:1054634/1‑139 (MQ=255)
                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       >  1:924179/1‑38 (MQ=255)
                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                       <  2:924179/38‑1 (MQ=255)
                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACg  >  2:1800882/1‑139 (MQ=255)
                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACg  >  2:648951/1‑139 (MQ=255)
                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACg  <  1:648951/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACG  >  NC_000913/3014150‑3014418

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