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.7% 249.2 / ‑5.7 78*257* (TAG→TAAyqeCputative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (41/36);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGG  >  NC_000913/3014157‑3014419
                                                                                                                                  |                                                                                                                                    
gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                                             >  1:1202519/1‑139 (MQ=255)
      cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                                                       >  1:1069391/1‑139 (MQ=255)
            tCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGcgcgcg                                                                                                                 <  2:1069391/139‑1 (MQ=255)
                 aaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                    >  1:329084/1‑115 (MQ=255)
                  aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                           <  1:343675/139‑1 (MQ=255)
                  aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                           <  1:747087/139‑1 (MQ=255)
                     tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        <  2:1397341/139‑1 (MQ=255)
                     tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        <  1:100534/139‑1 (MQ=255)
                     tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGCGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        <  2:1573777/139‑1 (MQ=255)
                      gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                       >  2:1202888/1‑139 (MQ=255)
                      gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                       >  1:665385/1‑139 (MQ=255)
                      gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                       >  2:656665/1‑139 (MQ=255)
                       cgacATAGCAAACAGACAATTTTGCTCTTCGAGGTTTGCAGCCTCTGTGAAAATATTCATACCAACGTCAATGGATATCTTTCAGTAACCCGGAAAACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                      <  1:189966/136‑1 (MQ=255)
                            tAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                 >  1:364365/1‑139 (MQ=255)
                                    ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAta                                                                                                      >  1:1387744/1‑125 (MQ=255)
                                    ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAta                                                                                                      <  2:1387744/126‑2 (MQ=255)
                                      aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                                             <  1:208805/117‑1 (MQ=255)
                                      aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                                             >  2:208805/1‑117 (MQ=255)
                                             tGCTCTTCGAGTTTTGCAGCCTCTGTTAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCctagccag                                                                                <  2:106308/139‑1 (MQ=255)
                                              gCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGGATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                                                                               >  1:860788/1‑139 (MQ=255)
                                                    cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                                         <  2:1028363/139‑1 (MQ=255)
                                                          ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgt                                                                   <  1:1483620/139‑1 (MQ=255)
                                                            gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                                 <  2:607506/139‑1 (MQ=255)
                                                                    tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                         >  1:1408285/1‑91 (MQ=255)
                                                                    tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                         <  2:1408285/91‑1 (MQ=255)
                                                                    tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                         >  2:259801/1‑139 (MQ=255)
                                                                      tgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                       <  2:460260/139‑1 (MQ=255)
                                                                      cgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTTTTGTTGCa                                                       <  2:1386273/138‑1 (MQ=255)
                                                                              ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                   <  2:256233/71‑1 (MQ=255)
                                                                              ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                                   >  1:256233/1‑71 (MQ=255)
                                                                              ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                                                   <  1:445829/103‑1 (MQ=255)
                                                                              ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                                                   >  2:445829/1‑103 (MQ=255)
                                                                                cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                >  1:178290/1‑136 (MQ=255)
                                                                                cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                <  2:178290/136‑1 (MQ=255)
                                                                                cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                             >  1:1432976/1‑139 (MQ=255)
                                                                                cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGAAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                             >  1:478851/1‑139 (MQ=255)
                                                                                  tACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGa                                           >  1:1270668/1‑139 (MQ=255)
                                                                                   aCCAACCTCAATGGATATCTTTCAGTAACCCGGGATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                          <  2:860788/139‑1 (MQ=255)
                                                                                    ccAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        >  2:1651744/1‑76 (MQ=255)
                                                                                    ccAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        <  1:1651744/76‑1 (MQ=255)
                                                                                       aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                         >  1:1226953/1‑136 (MQ=255)
                                                                                       aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                         <  2:1226953/136‑1 (MQ=255)
                                                                                       aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                      >  2:1655605/1‑139 (MQ=255)
                                                                                        ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                          <  2:269576/134‑1 (MQ=255)
                                                                                        ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                          >  1:269576/1‑134 (MQ=255)
                                                                                         cTCGATGGATATCTTTCAGTGACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCTGTTGAAGGAGCTCGCt                                    >  1:1123431/1‑139 (MQ=255)
                                                                                         cTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                    >  2:1621688/1‑43 (MQ=255)
                                                                                             aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                     >  1:597827/1‑54 (MQ=255)
                                                                                             aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                     <  2:597827/54‑1 (MQ=255)
                                                                                                  atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCaa                            >  1:1099995/5‑139 (MQ=255)
                                                                                                    atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                          >  2:997175/3‑139 (MQ=255)
                                                                                                    atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                          >  1:38023/3‑139 (MQ=255)
                                                                                                    atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                          >  1:434949/3‑139 (MQ=255)
                                                                                                     tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                         <  2:1202519/138‑1 (MQ=255)
                                                                                                     tcTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                         <  1:1121373/138‑1 (MQ=255)
                                                                                                       tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        <  2:942249/58‑1 (MQ=255)
                                                                                                       tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                        >  1:942249/1‑58 (MQ=255)
                                                                                                       tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                       <  2:665385/139‑1 (MQ=255)
                                                                                                         tCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCAtt                     >  2:1278957/1‑139 (MQ=255)
                                                                                                           aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct                   <  1:656665/139‑1 (MQ=255)
                                                                                                                cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                       <  1:525679/98‑1 (MQ=255)
                                                                                                                cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                       >  2:525679/1‑98 (MQ=255)
                                                                                                                cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                      <  2:838705/115‑1 (MQ=255)
                                                                                                                cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                      >  1:838705/1‑115 (MQ=255)
                                                                                                                   ggAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                                    >  1:464327/1‑82 (MQ=255)
                                                                                                                   ggAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCACAGCCAGATTGCCTCTAcg                                                                    <  2:464327/82‑1 (MQ=255)
                                                                                                                       tACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                           >  1:449410/1‑39 (MQ=38)
                                                                                                                       tACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                                           <  2:449410/39‑1 (MQ=255)
                                                                                                                         cccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGccc                                                                                      <  2:789731/58‑1 (MQ=255)
                                                                                                                         cccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGccc                                                                                      >  1:789731/1‑58 (MQ=255)
                                                                                                                          ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCa                             >  2:518209/1‑114 (MQ=255)
                                                                                                                          ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCa                             <  1:518209/114‑1 (MQ=255)
                                                                                                                            gggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                                    >  1:460726/1‑73 (MQ=255)
                                                                                                                            gggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                                    <  2:460726/73‑1 (MQ=255)
                                                                                                                            gggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACgg  >  1:808841/1‑139 (MQ=255)
                                                                                                                               ccGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                      <  2:781453/36‑1 (MQ=255)
                                                                                                                               ccGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                      >  1:781453/1‑36 (MQ=255)
                                                                                                                               ccGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                                              >  1:633298/1‑60 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGG  >  NC_000913/3014157‑3014419

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: