Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,561,969 A→T 6.1% intergenic (‑146/‑44) glpE ← / → glpD thiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,561,9690AT 6.1% 707.0 / 27.4 197intergenic (‑146/‑44)glpE/glpDthiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (89/96);  new base T (6/6);  total (95/102)
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.01e-01

Only 100 of 203 total aligned reads displayed.

ATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCA  >  NC_000913/3561822‑3562116
                                                                                                                                                   |                                                                                                                                                   
   cGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                                                                                                                                                 <  1:6238/149‑1 (MQ=255)
     tctctcTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGa                                                                                                                                                     >  4:1018866/1‑143 (MQ=255)
     tctctcTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                                                                                                                               >  4:918112/1‑149 (MQ=255)
       tctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                              <  1:1814578/148‑1 (MQ=255)
          ctttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAta                                                                                                                                          >  2:471156/1‑148 (MQ=255)
           tttctttACAATCAAGTGGGCAAATTTACCACACTGTTTACGTCGAAGCGGCAGATCAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAGTGTTACCAAAAGCGCGATTCTTTGCTAATATCTTCGATAACGAACAttt                                                                                                                                         <  4:1171529/149‑1 (MQ=255)
           tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAtat                                                                                                                                         <  1:1864304/149‑3 (MQ=255)
           tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAtat                                                                                                                                         >  2:1864304/1‑147 (MQ=255)
           tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATaa                                                                                                                                                  <  2:102215/140‑1 (MQ=255)
           tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                         >  4:457360/1‑149 (MQ=255)
               tttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                                                                                                                                     <  3:186073/149‑1 (MQ=255)
                 tACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGc                                                                                                                                   >  4:474012/1‑149 (MQ=255)
                     aacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                                                                                                                               >  4:1102483/1‑149 (MQ=255)
                      acaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTaa                                                                                                                              <  3:250948/149‑1 (MQ=255)
                      acaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTaa                                                                                                                              <  2:577857/149‑1 (MQ=255)
                        aaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACg                                                                                                                            >  4:368466/1‑149 (MQ=255)
                           tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaaa                                                                                                                         <  4:1635671/149‑1 (MQ=255)
                              gCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTaa                                                                                                                              >  1:425219/1‑141 (MQ=255)
                                 aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAtatta                                                                                                                                       >  3:599020/1‑125 (MQ=255)
                                 aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAtatta                                                                                                                                       <  4:599020/129‑5 (MQ=255)
                                   tttACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGa                                                                                                                 <  3:474012/149‑1 (MQ=255)
                                     tACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAgg                                                                                                               <  4:105496/149‑1 (MQ=255)
                                      aCCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt                                                                                                                                       <  3:1032207/124‑1 (MQ=255)
                                                 tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                                   <  3:457360/149‑1 (MQ=255)
                                                    gTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGc                                                                                                                                   <  1:1586364/114‑1 (MQ=255)
                                                     tCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                               >  4:883921/1‑149 (MQ=255)
                                                      cGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaaa                                                                                              <  3:368466/149‑1 (MQ=255)
                                                      cGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaaa                                                                                              <  3:918112/149‑1 (MQ=255)
                                                         aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                           <  1:1415304/149‑1 (MQ=255)
                                                             gCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGa                                                                                                                                                     <  1:724106/87‑1 (MQ=255)
                                                             gCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGa                                                                                                                                                     >  2:724106/1‑87 (MQ=255)
                                                             gCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                                                                                                                                     >  4:1647636/1‑103 (MQ=255)
                                                               aGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAtt                                                                                     >  3:282131/1‑149 (MQ=255)
                                                               aGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAtt                                                                                     <  3:1154077/149‑1 (MQ=255)
                                                                 aTAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt                                                                                                                                >  2:1442028/1‑104 (MQ=255)
                                                                   aaaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGa                                                                                 <  2:155506/149‑1 (MQ=255)
                                                                    aaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaaa                                                                                                                         <  2:280942/108‑1 (MQ=255)
                                                                    aaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgc                                                                                                             <  2:45960/120‑1 (MQ=255)
                                                                    aaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAAGGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaaa                                                                                                                         >  1:280942/1‑108 (MQ=255)
                                                                     aCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATa                                                                               >  1:247546/1‑149 (MQ=255)
                                                                       gCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgg                                                                             <  1:1509891/149‑1 (MQ=255)
                                                                         cATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACa                                                                                                                                            >  1:364254/1‑84 (MQ=255)
                                                                         cATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                                        >  3:1643079/1‑104 (MQ=255)
                                                                         cATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAt                                                                                      >  1:111717/1‑138 (MQ=255)
                                                                         cATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAt                                                                                      <  3:1486354/138‑1 (MQ=255)
                                                                         cATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAt                                                                                      >  4:1486354/1‑138 (MQ=255)
                                                                          aTAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcag                                                                                                           <  2:1141257/116‑1 (MQ=255)
                                                                                ttATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCAt                                                                    <  2:247546/149‑1 (MQ=255)
                                                                                  ataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                           <  2:1785587/124‑1 (MQ=255)
                                                                                    aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAAt                                                                <  4:282131/149‑1 (MQ=255)
                                                                                        atCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaaa                                                                                              >  2:60015/1‑115 (MQ=255)
                                                                                         tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGt                                                             >  4:1865714/1‑147 (MQ=255)
                                                                                         tCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGc                                                           <  3:883921/149‑1 (MQ=255)
                                                                                           aCTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAAt                                                                                                                   <  2:41052/91‑1 (MQ=255)
                                                                                             tctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt                                                                                                                                       >  3:1869383/1‑69 (MQ=255)
                                                                                               tAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAGTGAGGGCAGCATGGAAACCAAAGATCTGa                                                                                       <  4:751983/115‑1 (MQ=255)
                                                                                               tAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGa                                                                                       >  3:751983/1‑115 (MQ=255)
                                                                                               tAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTAt                                                     >  4:107090/1‑149 (MQ=255)
                                                                                                  aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                              >  1:1275935/1‑57 (MQ=255)
                                                                                                   aTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGa                                                                                       <  2:1365055/111‑1 (MQ=255)
                                                                                                   aTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGCGAATGAGGGCAGCATGGAAACCAAAGATCTGa                                                                                       >  1:1365055/1‑111 (MQ=255)
                                                                                                         tttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCaa                                                                 <  2:1046107/127‑1 (MQ=255)
                                                                                                          ttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcg                                          >  4:797346/1‑149 (MQ=255)
                                                                                                               tGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCa                                                                     >  2:806025/1‑117 (MQ=255)
                                                                                                               tGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCa                                                                     <  1:806025/117‑1 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCGATGGTGCTGGTATCGCGGCAGAcgc                                         >  4:1232827/1‑139 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              <  3:681934/118‑1 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              >  4:681934/1‑118 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgc                                         <  3:1232827/139‑1 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgc                               >  1:496970/1‑149 (MQ=255)
                                                                                                                     cTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgc                               >  3:353694/1‑149 (MQ=255)
                                                                                                                           cgcgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgccgc                                      >  3:1112987/1‑136 (MQ=255)
                                                                                                                            gcgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTAt                        <  2:529673/149‑1 (MQ=255)
                                                                                                                             cgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGa                                                                                                     <  1:76753/71‑1 (MQ=255)
                                                                                                                             cgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGa                                                                                                     >  2:76753/1‑71 (MQ=255)
                                                                                                                                ttCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaaa                                                                                              >  1:579738/1‑75 (MQ=255)
                                                                                                                                ttCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGt                    <  3:797346/149‑1 (MQ=255)
                                                                                                                                ttCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGt                    >  2:97086/1‑149 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              >  1:884885/1‑104 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              <  3:1098126/104‑1 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              <  2:884885/104‑1 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTgg                                                        <  2:1862050/110‑1 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTgg                                                        >  1:1862050/1‑110 (MQ=255)
                                                                                                                                   tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGt                    <  1:1018665/146‑1 (MQ=255)
                                                                                                                                     tGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGa               <  3:1439102/149‑1 (MQ=255)
                                                                                                                                        tAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                              <  1:465798/99‑1 (MQ=255)
                                                                                                                                            atGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGa        <  3:1066978/149‑1 (MQ=255)
                                                                                                                                              gTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAc                                           >  2:704805/1‑112 (MQ=255)
                                                                                                                                              gTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAc                                           <  1:704805/112‑1 (MQ=255)
                                                                                                                                              gTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATc                       >  1:560984/1‑132 (MQ=255)
                                                                                                                                              gTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAgg      >  4:1449771/1‑149 (MQ=255)
                                                                                                                                               ttCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgccgc                                      <  4:1341271/116‑1 (MQ=255)
                                                                                                                                               ttCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgccgc                                      >  3:1341271/1‑116 (MQ=255)
                                                                                                                                                 cGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgt                                                                                   <  4:197734/69‑1 (MQ=255)
                                                                                                                                                 cGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAg       >  4:1150243/1‑145 (MQ=255)
                                                                                                                                                 cGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGgcgc   >  3:594512/1‑149 (MQ=255)
                                                                                                                                                 cGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGgcgc   <  1:97086/149‑1 (MQ=255)
                                                                                                                                                  gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCa  >  1:1476998/1‑149 (MQ=255)
                                                                                                                                                  gATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCa  >  3:1669708/1‑149 (MQ=255)
                                                                                                                                                   aTAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTg                >  1:1037761/1‑134 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCA  >  NC_000913/3561822‑3562116

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: